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PREFACE 

This  manual  is  intended  to  fill  a  long  standing  gap  in  the 
documentation  of  SYMBOLANG.   Until  now  there  has  been  no  reasonably 
detailed  reference  for  users  of  the  package.   Past  reports  have 
assumed  much  expertise  on  the  part  of  the  reader.   We  are  probably 
guilty  of  the  same  sin,  but  hopefully  to  a  lesser  extent. 

"Old  hands"  at  SYMBOLANG  will  note  that  there  have  been  major 
changes  made.   Functions  of  an  arbitrary  number  of  variables  are 
now  allowed.   Expressions  may  appear  as  exponents,  instead  of 
exponents  being  restricted  to  constants.   List  space  is  conserved 
by  allowing  expressions  to  share  common  subexpressions.   Evaluation 
and  substitution  have  been  combined  and  generalized.  Most  of  the 
old  routines  are  declared  unfit  for  new  code,  and  a  bunch  of 
strangers  have  taken  their  place.   Yet,  despite  the  brief  flirtation 
with  the  name  SYMB66,  we  still  call  the  package  SYMBOLANG. 

Arnold  Lapidus  and  Max  Goldstein,  the  original  developers  of 
SYMBOLANG,  provided  an  algebraic  manipulation  package  with  at  least 
one  major  virtue  -  it  could  grow  freely.  The  internal  representation 
has  accommodated  itself  well  to  extensions  and  could  probably  accept 
still  more.   Being  a  collection  of  subroutines  for  use  by  FORTRAN 
programs,  the  package  can  always  absorb  new  features  by  simple 
addition  of  new  routines.   It  is  very  likely  that  SYMBOLANG  will 
continue  to  grow.   At  least  tl»e  author  hopes  that  to  be  the  case. 

The  package  itself  is  the  work  of  Arnold  Lapidus,  Max  Goldstein, 


Susan  S.  tloffberg,  and,  recently,  the  author.  Great  assistance  has 

been  provided  by  Samuel  Greenspan  and  Alfred  Magnus,  and  others 

who  have  been  at  the  AEC  Computing  and  Applied  Mathematics  Center 

of  New  York  University.   A  singular  debt  is  owed  to  Professor  Goldstein 

for  his  inspiration  and  continued  work  on  SYMBOLANG. 

The  author  also  wishes  to  thank  Professor  Nicholas  Findler  for 
providing  a  considerable  part  of  the  incentive  for  the  writing  of 
this  manual,  and  Frances  C.  Bernstein  for  many  hours  of  patient 
proofreading,  indexing,  and  wifely  cheer. 

Herbert  J.  Bernstein 

New  York ,  New  York 
February  1970 
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1.   INTRODUCTION 

This  is  a  user's  manual  for  the  algebraic  manipulation  package 
SYMBOI.ANG,  a  collection  of  FORTRAN  callable  subroutines  to  perform 
arithmetic  operations,  substitutions,  evaluations,  and  differentiations 
on  expressions  represented  as  SLIP  lists.  This  is  not  a  text  on 
algebraic  manipulation,  nor  even  on  algebraic  manipulation  with 
SYMBOLANG,  but  a  description  of  the  external  features  and  some 
internal  features  of  version  3.34  of  SYMBOLANG  as  implemented  on 
a  Control  Data  6600  operating  under  the  SCOPE  3.1.2  operating  system 
using  the  RUN  2.3  FORTRAN  compiler. 

The  reader  is  assumed  to  have  access  to  and  some  familiarity 
with  references  [1],[2]  and  [14]  which  describe  SCOPE  3.1.2,  RUN  2.3 
and  SLIP  respectively. 

1.1.   An  Overview 

The  manual  is  divided  into  three  major  sections:   INTRODUCTION 
which  discusses  background,  system  access,  general  structure,  etc., 
BASIC  ROUTINES  which  describes  the  most  commonly  used  SYMBOLANG 
calls,  and  BELLS,  WHISTLES,  AND  FRILLS  which  covers  features  which 
should  not  often  be  needed.   Each  section  is  further  subdivided  into 
subsections,  some  of  which  end  with  exercises.   It  is  the  author's 
intention  that  the  user  will  at  least  once  read  the  manual  linearly 
rather  than  randomly,  attempt  tlie  exercises,  and  ponder  the  examples. 
However,  as  a  concession  to  those  who  can  only  read  manuals  randomly, 
an  effort  has  been  made  to  provide  sufficient  redundant  information 
to  permit  considerable  skipping.  Thus  the  phrases  "list  or  expression 
lis"   and  "variable  or  array  element  var"   are  semi-inf initely  repeated 
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rather  than  rely  on  the  reader's  memory  of  notation. 

The  bulk  of  the  information  presented  consists  of  descriptions 
of  functions  and  subroutines.   The  user  who  feels  a  need  for  some 
sort  of  a  preprocessor  to  facilitate  his  coding  efforts  should 
consult  [9]  for  a  description  of  the  SYMBOLANG  preprocessor  QUEST, 
[11]  for  a  simple  genrral  macro  processor  which  might  be  useful 
for  this  purpose,  or  the  example  in  §3.10  which  happens  to  be  a 
preprocessor  for  SYMBOLANG  written  in  SYMBOLANG. 

The  attitude  is  assumed  throughout  that  a  function  is  a 
sul)routine  and  that  a  subroutine  of  at  least  one  argument  may  he 
treated  as  a  function.   I.ocv  [10,  pp.  52-53]  provides  a  nice 
desciption  of  the  results  of  such  an  attitude.   Those  reading  this 
manual  for  the  purpose  of  converting  SYMBOLANG  to  another  environment 
should  be  warned  ttiat  many  compilers  consider  this  attitude  subversive. 
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1.2.   Background 

SYMBOLANn  was  developed  by  Lapidus  and  Goldstein  [8]  in  1965 
for  use  on  an  IBM  7094,  extensively  revised  in  1967  by  Lapidus, 
Goldstein,  and  lloffberg  [9]  for  use  on  a  CPC  6600  under  the  Ciiippewa 
Operating  System,  and  with  major  revisions  and  extensions  brought 
to  its  current  form  in  1969  (see  [3]).   In  its  three  incarnations 
it  has  lx3en  a  pact-age  of  FORTRAN  callable  subroutines,  mostly  written 
in  FORTRAN  and  based  on  the  list  processing  package  SUP  developed 
by  Weizenbaum  [M],  which  allows  one  to  perform  such  operations 
as  multiplying  the  expression  f[l  +  x]]    by  tlie  expression  [[1  -  x]] 
to  obtain  the  expression  [[1  -  x^]]    rather  than  some  numerical  value. 
That  is,  SYMBOLANG  is  one  of  the  large  class  of  "formula  manipulation' 
or  "algebraic  manipulation"  systems  such  as  are  described  in  [5] ,  [6] , 
[12],  and  [13]. 

Wiien  dealing  with  numbers  it  is  reasonable  to  allocate  a  certain 
fixed  amount  of  space  to  represent  each  number,  and  just  throw  away 
any  digits  that  overflow  the  allowed  space  as  insignificant.   With 
expressions,  results  may  grow  and  shrink  over  a  very  wide  range, 
and  we  rarely  know  what  may  be  safely  discarded.   Rather  than 
adopt  the  rather  expensive  approach  of  allowing  the  maximum  possible 
space  for  each  expression,  it  is  common  in  formula  manipulation 
to  use  list  processing  techniques  which  allocate  to  an  expression 
just  the  space  it  needs  and  allow  one  to  return  unneeded  space  to 
a  list  of  available  space  for  future  use. 

The  list  processing  system  chosen  for  SYMBOLANG  is  SLIP,  itself 
a  package  of  FORTRAN  callable  subroutines  mostly  written  in  FORTRAN. 
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The  SLIP  routines  allow  one  to  create  lists  of  items,  which  may 
themselves  be  lists,  insert  and  delete  items,  and  scan  the  items 
on  a  list.   It  is  most  convenient  to  think  of  a  SLIP  list  as  a 
parenthesized  group  of  quantities: 

(  ) 

or   iquani  quan2    ■••   quan,) 

where  we  call  the  first  list  empty,  and  the  second  list  a  list  of 
k   elements,  the  leftmost  (top)  element  being  the  quantity  quan\ 
and  the  rightmost  (bottom)  clement  quan,.      In  certain  contexts  we 
shall  also  consider  the  open  and  close  parentheses  as  a  special 
element  of  a  list,  the  "head",   livery  list,  even  an  empty  list, 
has  a  head. 

A  SLIP  list  is  referenced  by  a  quantity  which  is  interpreted 
by  the  routines  of  the  package  as  a  pointer  to  the  head  of  the  list, 
which  itself  contains  pointers  to  the  left-  and  rightmost  elements 
of  tl\e  list,  which  themselves  contain  pointers  to  the  elements  to 
their  left  and  right,  etc.   We  use  the  symbol  lis,   possibly  with 
a  subscript,  to  denote  a  quantity  which  we  intend  to  be  treated  as 
a  list  within  I'ORTRAN  code.   As  we  liave  done  above,  we  use  the 
symbol  quan,    for  quantities  on  wliich  we  place  no  a  priori  restrictions, 
Such  a  quantity  might  be  a  floating  point  number,  or  an  integer,  or 
a  list,  or  a  Hollerith  constant,  etc. 

Mien  we  wish  a  more  restrictive  notation,  we  use  int   for  an 
integer  quantity,  hoi   for  a  Hollerith  quantity,  symb   for  a  Hollerith 
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quantity  without  embedded  blanks  (called  an  "expression  symbol"), 
holarray   for  a  Hollerith  quantity  from  which  more  than  ten  characters 
worth  of  information  may  be  obtained,  var   for  a  variable  or  array 
element,  intvav   for  an  integer  variable  or  array  element,  and  array 
for  an  array.  The  user  should  refresh  his  memory  of  FORTRAN 
sufficiently  to  recall  that  a  quantity  may  be  a  literal,  variable, 
array  element,  or  value  of  a  function  call,  and  that  values  of 
function  calls  may  be  used  in  places  where  variables  or  array  elements 
are  called  for  as  arguments  (causing  the  possible  loss  of  that  value). 

Since  a  SLIP  list  may  be  placed  on  a  SLIP  list,  certain 
precautions  are  needed  to  avoid  losing  track  of  lists.   First,  a 
list  must  never  be  placed  on  itself  (see  [15]).   Second,  quantities 

other  tlian  lists  to  be  placed  on  a  list  are  restricted  to 

-293 
floating  point  quantities  in  the  range  10 

322 
to  10    in  magnitude 


integer  quantities  in  the  range  -(2   -1)  through  2   -1 
Hollerith  quantities  (up  to  10  characters) 


We  represent  expressions  as  SLIP  lists  as  follows; 

<expression>  *-   (  {<tenn>}o  ) 
<term>  <-  (  quan   {<factor>  <power>}o  ) 
<factor>  *-   {<expression>}o  symb 
<power>  *■  quan    \    <expres5ion> 
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where  the  left-hand  side  of  each  line  represents  the  class  of 
objects  being  defined,  <  and  >  enclose  class  names,  {  and  }  indicate 
that  the  enclosed  objects  are  to  be  included  in  the  definition  for 
each  number  of  repetitions  in  the  range  indicated  by  the  numbers  on 
the  right  brace,  |  separates  alternative  definitions,  quan   is  treated 
as  the  class  of  floating  point  quantities  other  than  zero,  and 
S]imh   is  the  class  of  expression  symbols.   The  meaning  of  this 
syntax  is  simple:   an  expression  is  a  list  of  terms  and  represents 
the  sum  of  those  terms;  a  term  is  a  list  beginning  with  a  constant 
coefficient  followed  by  factors  with  powers  and  represents  the 
product  of  the  constant  by  the  factors  raised  to  the  powers;  a  factor 
is  a  list  fragment  consisting  of  expressions  followed  by  an  expression 
symbol  and  represents  the  expression  symbol  applied  as  a  function  to 
the  expressions  as  its  arguments  or  represents  the  expression  symbol 
itself  if  there  are  no  arguments;  and  a  power  is  either  an  expression 
or  a  constant. 

For  convenience  in  displaying  examples  of  expressions  as  lists, 
we  introduce  commas  between  list  elements,  and  use  double  brackets 
as  above  to  enclose  real  world  expressions. 

The  expression  [[0]]  may  be  thought  of  as  a  sum  of  no  terms. 
We  represent  it  as  the  list  (  )  which  is  empty. 

The  expression  [[36.5]]  is  another  example  of  a  "constant" 
expression;  i.e.  it  involves  only  numbers, not  symbols.   It  consists 
of  a  single  term  which  has  only  a  coefficient.   Ke  represent  it  as 
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the  list  (  (36.5)  )  which  has  one  element,  a  sublist.  The  sublist 
has  one  element,  the  number  36.5. 

The  expression  [[x]]    is  a  non-constant  expression  consisting 
of  one  term  with  coefficient  one,  a  single  factor  IHX,  and  the  power 
one.   Ke  represent  it  as  the  list  (  (1.,  IHX,  1.)  ). 

The  expression  [  [-18st)i(x)-  z^  ♦  1]]  is  a  moderately  complex 
expression,  which  we  represent  as  the  list  (  (1.),  (-18.,  IIIZ,  3., 
(  (1.,  IHX,  1.)  ),  3HSIN,  (  (1.,  IIIY,  1.)  )  )  ). 

The  expression  [[ui3er(0,  x)]]    involves  a  function  of  more  than 
one  variable.   We  represent  it  as  the  list  (  (1.,  (  ),  (  (1.,  IHX,  1.)  ) 
4HUSER,  1.)  ). 

In  general,  we  choose  a  Hollerith  constant  of  from  1  to  10 
non-blank  cliaracters  to  represent  each  symbol  in  the  expression, 
represent  all  function  arguments  as  lists,  break  up  the  expression 
into  a  sum  of  terms,  and  each  term  into  a  product  of  factors 
raised  to  powers,  and  introduce  the  trivial  coefficient  and  power 
one  where  necessary.   All  numbers  are  in  floating  point,  and  usually 
one  uses  HOKTRAiN  variable  naming  conventions  for  forming  expression 
symbols.   Rather  than  force  ourselves  to  use  the  distributive  law 
to  expand  parenthesized  subexpressions,  which  may  not  even  be  possible 
with  fractional  powers,  we  treat  parentheses  without  a  function  name 
as  an  appearance  of  a  special  function  name  which  we  arbitrarily 
translate  into  the  expression  symbol  3H1.*. 
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Tn   this  manner  we  may  represent  expressions   as  general   as 


m  m(l.O)      ^0,j(i,0)        m(l.l)  p^   j^^    ^) 


j'(i,i)=o      ^>n^,^) 


m(l,2)  Pp    .(.    ps 

J(l,2)=0      2,j(i,2)^^2,J(i,2),l'    ^2,j(l,2),2^ 


m(i,r^) 


jdlrJ^O    V,J(l,r.)^^r.,j(i,r^),l> 


1  ■ 


,    e         .,.  .  )       ^  ^ 
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as   the   following   list: 


m(l,0) 


°i^^O,j(i,0)    ^0,j(i,0)^j^.^Q)^Q 


•l,j(i,l)    ^l,j(l,l)    ^l,j(l,l) 


m(i,l) 
j(i,l)=0 


{e,     ,,,     ,,    3H1.*    p. 


} 


m'(i,l) 
j'(l,l)=( 


m(l,2) 
^^2,j(l,2),l    ^2,j(l,2),2    ^2,j(i,2)    P2 , j (1 , 2 ) ^ ^ ^ . ^ 2 ) =0 


{e 


r^,j(i,r^),l    ' 


m(i,r^) 
'i,J(i,r.),r.    V,j(i,r.)    ^^^i  ,  J  (i  ,r .  )  ^  .  ^ .  ^^_  ^^^ 


1  =  0 
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What  has  been  presented  so  far  does  not  force  any  degree 
of  uniqueness  in  the  representation  of  an  expression  as  a  list. 
That  matter  is  handled  in  §1.4.   For  most  purposes  the  user 
need  not  concern  himself  with  this  matter.   Very  complex  expressions 
may  be  formed  by  use  of  the  routine  INLIST  described  in  §2.2  and 
§3.2.   Reasonably  complex  expressions  may  be  built  "by  hand"  from 
simple  expressions  of  just  one  term  with  at  most  one  factor  by 
addition,  multiplication,  and  taking  powers  using  the  arithmetic 
routines  in  §2.6  and  the  list  construction  routines  in  §2.1. 
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1.3.   Residual  Notational  Matters 

We  are  about  to  look  at  some  actual  functions,  so  notation 
had  best  be  made  clear. 

An  expression  within  SYMBOLANG  code  is,  of  course,  just  a 
particular  form  of  a  list.   We  use  the  symbol  lis   to  represent  an 
expression  as  well  as  any  other  list.   However,  there  are  many 
routines  which  will  hang  or  do  other  interesting  things  if  presented 
with  an  arbitrary  list  where  they  expect  an  expression.   Further, 
there  are  many  list  modifications  which  are  perfectly  proper  for 
lists  other  than  expressions,  but  which  would  cause  too  much  dis- 
organization if  applied  to  expressions.   Nondestructive  operations 
specified  for  lists  may  always  be  applied  to  expressions.  Thus  the 
user  should  be  careful  to  distinguish  the  phrases  "list  lis", 
"expression  lis",    and  "list  or  expression  lis". 

In  most  cases  where  a  Hollerith  quantity,  hoi   or  symb,    is 
called  for,  only  the  first  ten  characters  will  be  used.   Also,  by 
Hollerith,  we  mean  exactly  that,  nnH,  not  nnL  or  nnR. 

Wliere  the  difference  is  important,  the  letter  0  will  be  distinguished 
from  the  digit  0  by  writing  the  letter  as  0.   Thus  we  use  SYMB0LANr)  or 
SYMBOLANG,  hut  not  SYMBOLANG. 

Some  functions  which  return  lists  or  expressions  liave  floating 
point  names.  Most  have  integer  names.  The  user  is  cautioned  to  avoid 
undesired  mode  conversion. 
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1.4.   Ordering  and  Simplification 

All  right  thinking  users  will  agree  that  A*B*C-C*B*A  is  best 
treated  as  zero.   Not  everyone  will  agree  that  X-(X)  is  best 
treated  as  zero.   SYMBOLANC  is  right  thinking  but  not  everyone. 
Identical  terms  differing  only  in  their  coefficients  are  combined, 
and  a  zero  coefficient  causes  a  term  to  be  dropped  from  an  expression. 
Identical  factors  differing  only  in  their  powers  are  combined,  and 
a  zero  power  causes  a  factor  to  be  dropped.   If  the  constant  expressions 
[[0)]  and  [[1]]  are  raised  to  any  power,  they  arc  unchanged.   In 
short,  when  operating  on  expressions,  SYMBOLANG  performs  complete 
identity  removal  in  the  sense  of  [13,  p.  557].   In  order  to  keep 
this  task  finite,  SYMBOI.ANC  imposes  an  ordering  on  all  expressions 
and  their  components.   All  operations  on  expressions  preserve  this 
ordering,  some  assume  it,  and  most  impose  it. 

'Hie  user  may  simplify  a  hand  created  expression  by  calling 
the  routine  SMPL. 

SMPL(Hs)  reorganizes  the  expression  lis    into  simplified 
form,  returning  lis   as  its  value.  This  is  an  actual  change  to  the 
expression,  and  one  of  the  few  cases  where  it  is  proper  to  alter 
an  existing  expression. 

This  operation  is  rather  expensive,   and  should  be  avoided  as 
unnecessary  in  most  cases,  since  simplification  spreads  in  SYMBOLANG 
with  lightning  speed. 
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The  ordering  imposed  is  built  up  from  orderings  on  numbers 
and  expression  symbols.   Numbers  are  ordered  arithmetically,  and 
expression  symbols  lexicographically  using  the  following  ordered 
alphabet: 

0123456789ABCDEFf; 
HIJKLMN0PQRSTUVWX 
YZ  +  -*/()$=,.   =   []:^ 

-►VAt4'<>£>_— j 

with  blank  coming  before  all  other  characters.  Thus  3M1.*  is 
less  than  IHA,  which  is  less  than  2HAA,  which  is  less  than  1H$. 
We  now  build  a  lexicographic  ordering  of  expressions  from 
the  following  rules: 

A.  All  numbers  are  less  than  all  expressions. 

B.  The  empty  expression  C  )  is  less  than  all  non-empty 
expressions. 

C.  An  expression  (terrni,   term2,    ...,  term,)    is  less  than 
the  expression  {term'\,    term'2,    ....  term',,)    if 

(1)  term\    is  less  than  tem'i;   or 

(2)  there  is  an  i   such  that  term,   exists  and  is  less 

than  term',   which  also  exists,  and  term.  -  term'  .   for 
^  ,13 

,7  =  1 ,  2,  .  .  .  ,  t-1;  or 

(3)  k   is  less  than  k'   and  term.  =  term'  .   for  ,?  =  1,  .  .  .  ,k. 

.7        ,7 

I).   A  term  with  just  a  coefficient  is  less  than  a  term  with 
factors,  and  also  less  than  any  other  term  with  just  a 
coefficient  where  that  coefficient  is  greater. 
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E.  A  term  (aoe ,  factor\,  power i ,    ....  factor,,  powerp) 
is  less  than  the  term  (ooe' ,   factor' i,  poiier'\,    ..., 
factor'^,,  power' ^,)    if 

(1)  factory    is  less  tlian  factor'\;   or 

(2)  there  is  a  leftmost  factor  or  power  in  the  first 
term  which  is  less  than  the  same  numbered  factor  or 
power  in  the  second  term,  and  where  all  factors  and 
powers  to  the  left  agree;  or 

(3)  k   is   less  than  k'   and  all  factors  and  powers 
through  the  kth   agree;  or 

(4)  the  terms  differ  only  in  their  coefficients,  and 
the  first  term  has  a  lesser  coefficient. 

F.  A  factor  without  arguments  is  less  than  any  factor 
with  arguments,  which  in  turn  is  less  than  any  factor  with 
a  still  greater  number  of  arguments. 

fi.  Two  factors  with  an  equal  number  of  arguments  are 
compared  first  by  their  expression  symbols.   If  those  are 
equal,  then  the  arguments  are  compared  from  left  to  right. 

Factors  within  terms  and  terms  within  expressions  are  reordered  to 
achieve  expressions  of  the  lowest  possible  ordering.   Constant 
expression  powers  are  replaced  by  the  constants,  terms  and  factors  are 
combined  where  possible,  and  dropped  when  trivial.   Parentheses 
are  not  handled  any  differently  from  any  other  function. 

The  result  is  to  produce  expressions  in  which  constant  terms 
lead,  symbols  without  arguments  come  next,  then  functions  of  one 
variable,  two,  etc. 

14 


1.15 


1.5.   System  Access 

The  user  gains  access  to  SYMBOLANG  by  running  ordinary  batch 
production  jobs.   It  may  also  be  possible  to  use  SYMBOLANG  in  a 
batch- like  time  sharing  environment  as  in  [7],  but  the  package  is 
not  interactive  in  the  sense  of  [6] . 

Assuming  normal  installation  procedure  has  been  followed,  some 
of  the  routines  will  have  been  added  to  the  system  library,  and  the 
rest  will  be  obtained  by  use  of  the  programs  SLIP  [3]  and  SELECT  [4] 
which  should  themselves  be  in  the  system  library.   In  such  an 
environment,  a  SYMBOLANG  job  would  have  the  structure 

Job  card 

RUN(S)        C0MPILP.  USER  PR0GRAM  T0  FILE  LG0 

SLIP(LIB)     PUT  SLIP  AND  SYMB0LANG  LIB  0N  LIB 

SELECTCN)     TRANSFER  NEEDED  R0UTINES  FROM  LIB  T0  LG0 

LG0.  L0AD  AND  EXECUTE  LG0 

end-of- record  aard 

user  program  and  auhroutines 
end-of-record  card 

user  data 
end-of-file  card 

The  job  card  should  allow  a  field  length  of  at  least  70000  octal 
words.   Time  depends  completely  on  the  task.   The  user  program  and 
subroutines  should  follow  ordinary  FORTRAN  coding  conventions,  subject 
to  the  following  constraints: 
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The  program  should  provide  the  standard  input  unit, 
SLINPIJT,  and  the  standard  output  unit,  6LnUTPlJT,  by  starting 
with  a  program  header  card,  say,  of  the  form 
PR0GRAM  TEST(1NPUT,0UTPIIT) 

A  list  of  available  space  must  be  initialized,  see  §2.11, 
say  by 

DIMENSI0N  SPACE (5000) 

CALL  I NITAS (SPACE, 5000) 
This  call  to  INITAS  also  provides  100  "public"  lists  in  blank 
C0MM0N  which  may  be  used  for  scratch  work.   They  are  accessed 

by 

C0MM0N  AVSL,X(100) 
where  AVSL  is  a  pointer  not  to  be  tampered  with,  and  X  is  the 
array  of  public  lists.   X(99)  and  X(IOO)  are  used  extensively 
within  SYMBOLANC  and  should  not  be  touched. 

C0MM0N  blocks  with  labels  consisting  of  six  characters 
and  beginning  with  LSQ  should  not  be  created  by  the  user,  nor 
should  he  create  one  with  the  label  SYMSYS. 

In  addition  to  entry  points  described  in  this  manual,  the 
user  should  not  create  any  routines  with  names  of  six  characters 
beginning  with  LSQ,  LTQ,  XSQ,  or  XTQ. 

Tlie  generality  of  the  last  constraint  is  to  allow  for  the  future 
expansion  of  the  package.  The  entry  points  currently  in  the  package 
are: 
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ADD 

ADVLEL 

ADVLER 

ADVLL 

ADVLNL 

AUVLNR 

ADVLR 

ADVLWL 

ADVLWR 

ADVSEL 

ADVSER 

ADVSL 

ADVSNL 

ADVSNR 

ADVSR 

ADVSWL 

ADVSWR 

AL0AD 

B0T 

BREAK 

CPYTRM 

DELETE 

DVSIIM 

DVTRNt 

EVALUE 

GETC0E 

IlITENT 

INITAS 

INITRD 

INLIST 

INLSTL 

INLSTR 

INSBST 

INSUBT 

INTARG 

INTENT 

IRALST 

IRARDR 

ITSVAL 

K0KE 

LCNTR 

LDATVL 

LDIF 

LEXIC0 

LIST 

LI STAY 

LISTMT 

LIST0N 

L0CARG 

L0CT 

L0FRDR 

L00KUP 

LPNTR 

LRDRCP 

I.R[)R0V 

LSQADD 

LSQCDR 

LSQCMP 

LSQCNM 

LSOCPY 

LSQCXP 

LSQDEF 

LSODES 

LSODSF 

LSQERR 

LSQGAR 

LSQGNF 

LSQIDR 

LSQINI 

LSOLC0 

LSQMEX 

LSQMNL 

LSQMTM 

LSQNAM 

LSQ0UT 

LSQPNT 

LSQRAZ 

LSQSBS 

LSQTRC 

F,SnTYP 

LSQUNM 

LSQVAL 

l.SQVVL 

LSSCPY 

LSTEQL 

LVLRVl 

LVLRVT 

MADATR 

MADLFT 

MADNBT 

MADNTP 

MAURC.T 

MAKEDL 

MANY 

M0N0 

MRKLSS 

MR K EST 

MTDLST 

MTLIST 

NAMEDL 

NAMTST 

NEWB0T 

NtWT0P 

NEWVAL 

N0ATVL 

NUCELL 

NULSTL 

NULSTR 

NUMPY 

NXTLFT 

NXTRGT 

P0PB0T 

P0PMIU 

P0PT0P 

P0WER 

PR I  PUT 

PRLSTS 

PUT EST 

RCELL 

REED 

SAME 

SBST 

SEQLL 

SEQLR 

SEORDR 

SEQSL 

SEQSR 

SMPL 

S0I,VE 

SRTRM 

SUB 

SUBSBT 

SUBST 

SUBSTP 

SIIBT 
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SUMPY     TERM      T0P 
TSTC0N    VALARG    VISIT 


TRCAL 


TRDNC 


The  following  entry  points  are  usually  obtained  from  the  system 
library  and  should  not  be  redefined: 


ABN0RML 

ACC0F.R 

ADVIN. 

AL0r, 

AL0niO 

AND 

ATAN 

BACC11K= 

BKSPRU. 

CI01. 

C0NT 

C0S 

DAT. 

END 

EQUAL 

EXIT 

LXP 

FIZBAK. 

FIZBA. 

GETBA 

K01)l:R 

ID 

IFENDF 

IFTIIEN 

INCHEK 

IND 

INHALT 

INIBIP 

INPIJTC 

INPUTS 

INTGER 

KRAKER 

LAN0RM 

LNKL 

LNKR 

MA1)0V 

MVWDS . 

0PF,N. 

0R 

0UTPTC 

0UTPTS 

0VWRT 

P0SFIL. 

P0SFI. 

Q8NTRY 

RBAIEX 

RBAREX 

RDPRII. 

REVIND 

SB ARCS 

SETDIR 

SETIND 

SHTRAY 

SHIFTR 

SHIN 

SIN 

SI0. 

SQIN 

SO0UT 

ST0P 

STRDIR 

STRIND 

SYSTEM 

SYSTEMC 

SYSTEMP 

TAN 

THEN IF 

ZER0 

This  list  is  a  mixture  of  SLIP,  SYMBOLANG,  and  miscellaneous 
routines,  not  all  of  which  will  be  described.   Unless  there 
is  some  pressing  reason  to  do  otherwise,  we  shall  refer  to 
all  as  SYMBOLANG  routines  on  the  theory  that  SYMBOLANG  is  an 
extension  of  SLIP  which  is  an  extension  of  FORTRAN. 
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1.6.   Elementary  Programming  Techniques 

The  most  commonly  repeated  pattern  of  code  in  SYMBOLANC 
programming  is 

create  operands 

perform  expression  operation 

erase  unneeded  operands 
For  example,  using  the  input  techniques  of  §2.2,  addition  as  per 
§2.6  and  erasure  as  per  §2.4,  we  may  input, add  and  erase  two 
expressions  by 

LA=INLIST(LA,5LINPUT) 

LB=INLIST(I,B,5LINPUT) 

LC=LSQAUD(LA,LB) 

CALL  LSQDF.S(LA) 

CALL  LSQUES(LB) 
The  erasure  of  unneeded  operands  is  very  important  to  avoid  running 
out  of  space  from  which  to  create  expressions. 

Since  mode  conversions  are  often  not  desired,  information 
may  be  transferred  without  conversion  by 

STRDlR^quan,    var)   which  stores  the  quantity  quan   into 
the  variable  or  array  element  var   without  mode  conversion. 

imG\iR{quan)   provides  the  quantity  quan   with  an  integer  name. 
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M0N0{(7Man)  also  provides  the  quantity  quan   with  an 
integer  name. 

REMiquan)   provides  the  quantity  quan   with  a  floating 
point  name. 

SM-\E(_quan)    also  provides  the  quantity  quan   with  a 
floating  point  name. 

The  final  programming  technique  to  be  covered  is  recursion. 
The  routines  TnW'i  and  VISIT  provide  the  ability  to  recursively 
execute  blocks  of  FORTRAN  code.   The  entry  to  the  code  is  provided 
by  executing  an  ASSIGN  statement  to  store  the  location  of  the 
first  line  of  code  to  be  executed  in  an  integer  variable.   Since 
the  fonn  of  an  ASSIGN  statement  is  usually 

ASSIGN  100  T0  L0C 
a  variable  established  this  way  will  be  designated  by  the  symbol 
too. 

VlSlTiloc)    saves  the  location  from  which  the  call  is  made 
and  transfers  control  to  the  FORTRAN  statement  at  the  location 
specified  by  loc.      VISIT  return  a  value  which  is  provided  as  the 
argument  to  the  call  to  TERM  which  returns  control. 

JV.RMiquan)    returns  control  to  the  location  from  which 
the  last  VISIT  was  paid,  discarding  that  location  from  the  list 
on  which  it  was  saved.  The  quantity  quan   is  returned  as  the  value 
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of  VISIT. 

For  example,  the  following  subroutine  will  compute  the  function 
N*(N-1)*(N-2)*...*2*1  recursively. 

FUNCTION  IFAC(N) 
ASSIGN  100  T0  L0C 
M=N 

CALL  STRDIR(VISIT(L0C).IFAC) 
RETURN 
100   IF(M.LE.1)CALL  TERM(l) 
M=M-  1 

CALL  STRDIR(VISIT(L0r),K) 
M=M-»- 1 

CALL  TERM(M*K) 
END 

Note  that  we  have  been  careful  to  restore  M  to  its  original  state 
on  each  pass  through  the  code.  More  complex  saving  and  restoring 
may  be  done  with  the  routines  MANY  and  LSQUNM  described  in  §2.1 
and  §2.4. 
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1.7.   Not  so  Elementary  Programming  Techniques 

SYMBOLAKG  routines  make  use  of  some  special  facilities  with 
which  the  user  may  wish  to  enhance  his  own  code. 

The  sliort  simple  routines  such  as  STRDIR,  which  fit  into  one 
word  of  machine  code,  actually  replace  the  code  used  to  call  them 
with  in-line  code.   This  is  done  when  the  routine  is  executed,  so 
that  a  loop  of  sucli  calls  will  execute  the  routine  once  and  the 
in-line  code  thereafter.   (This  execution  time  overwriting  is 
believed  to  be  unique  to  SYMBOLANG) .   The  effect  on  most  SYMBOLANG 
code  is  an  improvement  of  30%  over  the  timing  with  more  conventional 
coding.   If  the  user  has  a  similar  short  routine,  he  may  make  it 
overwrite  its  call  with  the  routine  0VWRT. 

f^VYIRT iguan)    replaces  the  call  to  the  routine  which  called 
0VWRT  with  the  word  of  machine  code  given  by  the  quantity  quan.      For 
example  0VWRT(46OOO46OOO46OOO46OOOB)  would  replace  the  call  with  a 
do-nothing  instruction. 

Some  compression  of  code  can  be  achieved  by  use  of  the  functions 
IFTHEN  and  TIIENIF  which  act  as  conditional  expressions. 

IFTUENiquani,   quan2 quan.)    for  k  =   1,  2,  3,  ....  60, 

considers  its  odd  numbered  arguments  (except  possibly  the  last)  to 
be  logical  expressions,  and  returns  the  leftmost  even  numbered 
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argument  following  an  odd  numbered  argument  which  is  TRUE.   If  the 
total  number  of  arguments  is  odd,  the  last  argument  is  returned 
if  all  the  other  odd  numbered  arguments  are  FALSE.   If  the  total 
number  of  arguments  is  even,  and  no  odd  numbered  argument  is  TRUE, 
zero  is  returned. 

TIIENIF  is  identical  to  IITIIEN,  except  for  the  name. 

For  example 

X=TIIENIF(J.LT.O,0. ,J.GT.5,5. ,FL0AT(J)) 

is  a  short  way  of  writing 

X  =  FL0AT(J) 

IF(J.LT.O)X=0. 

IF(J.GT.5)X=5. 

IFTHEN  and  THENIF  are  examples  of  SYMBOLANG  routines  which 
handle  variable  length  argument  lists.  The  user  may  write  routines 
which  handle  such  argument  lists  by  use  of  the  functions  SBARGS, 
L0CARG,  INTARG,  and  VALARG,  as  well  as  the  subroutine  AL0AD. 

SBARGS (tntvar)  stores  the  number  of  actual  arguments 
used  in  the  call  to  the  routine  which  called  SBARGS  in  the  integer 
variable  or  array  element  intvar.      The  number  of  arguments  is  also 
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returned  as  the  function  value. 

AL0AD(ar*ray ,  var\,   var>2,    ...,   y'^ir,)  for  k  =   1,  2,  ...,  59, 
stores  the  machine  addresses  of  the  variables  var .   into  the  array 
elements  arrayii)   where  array   is  an  array  of  dimension  k.    Such  a 
call  to  AL0AD  should  be  made  to  store  the  addresses  of  the  first 
six  arguments  of  any  routine  which  wishes  to  use  I,0CARG,  INTARG,  or 
VALARG. 

L0CARG(int)  returns  the  machine  address  of  argument  int 
of  the  routine  which  called  L0CARG.  The  intefjer  quantity  int   must 
lie  between  1  and  the  actual  number  of  aguments  used,  and  a  call 
to  AL0AD  for  the  first  six  arguments  must  have  been  made. 

INTARG (i>7t)  returns  the  value  of  argument  int   of  the 
routine  which  called  INTARG.   The  same  restrictions  apply  as  do  to 
I,0CARG. 

VALARG  is  INTARG  with  another  name. 

Thus  a  subroutine  which  wishes  to  index  through  its  argument  list 
might  begin 

FUNCTI0N  G0G0(A,B,C,D,E,F,G,H,I,J,K,L,M,N,0) 

DIMI-NSI0N  LAF^G(6) 

CALL  SBARGS(NARGS) 

CALL  AL0AD(LARG,A,B,C,D,F.,F) 
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100  JJ=1,NARGS 


.  .  =  VALARG(JJ)  .  , 


Since  machine  addresses  have  been  mentioned,  we  will  look  at 
a  few  routines  to  make  some  use  of  them.   In  the  rest  of  this  manual 
the  concept  of  machine  address  will  be  avoided  assiduously. 

STRIND(<7uan,  tnt)    stores  the  quantity  quan   into  the 
machine  address  specified  by  the  integer  quantity  int. 

MAI)0V(i»ar)  returns  the  machine  address  of  the  variable 
or  array  element  var.      Thus  STRIND(yar,  MAP0V(yar')  )  is  a  peculiar 
way  to  do  nothing. 

C0NT(int)  returns  the  contents  of  machine  address  int. 

INIlALTCtnt)  also  returns  the  contents  of  machine  address 
int,   but  with  an  integer  name. 

These  last  four  routines  are  examples  of  SLIP  "primitives".   The 
full  set  consists  of  C0NT,  ID,  INIiALT,  LNKL,  I.NKR,  MAD0V,  SF.TDIR, 
SETINI),  STRDIR,  STRIND,  AND,  HQUAL,  INTHKR,  LAN0RM,  0R,  SHIN, 
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SQIN,  and  SQ0UT.   See  [10]  and  [14]  for  details  on  these  routines, 
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2.   BASIC  ROUTINES 

The  following  are  the  more  important  calls  to  SYflEOLANG  routines. 
Since  the  full  package  consists  of  over  150  routines,  some  with 
several  alternative  calling  sequences,  the  decision  as  to  which  are 
"basic"  is  somewhat  arbitrary. 

2.1.   Creation  of  Lists  and  Expressions 

The  functions  LIST,  LSQMNL,  MANY,  and  LSSCPY  may  be  used  to 
construct  lists  and  expressions  from  scratch  or  from  existing  lists. 

LIST(9)  returns  a  new  empty  list,  (  ),  v;hich  is  the  correct 
representation  of  [[0]]. 

LSQMNL  (qt/on  1 ,  qv.an2,    •••,  quan,)  ,   k   =   1,  2,  ....  60,  returns 
a  new  list  containing  the  quantities  in  the  argument  list,  i.e.  {qiAarti. 
qucm2,    ...,  quan,).      For  example,  the  following  return  the  correct 
representations  of  the  indicated  expressions: 

[[3.6]]  LSQMNL (LSQMNL (3. 6)) 

[[!]]  LSQMNL (LSQMNL (1.)) 

[[-6.9E36]]  LSQMNL (LSQMNL (-6. 9F36)) 

[[x]]  LSQMNL (LSQMNL(1. ,1HX,1.)) 

[[5y"^]]  LSQMNL  (LSQMNL  (5.  ,1I1Y,3.)) 


[[-4  .Ssinf^alpka)  ]] 


LSQMNL (LSQMNL (-4. 5, 

LSQMNL(LSQMNL(1. .SHALPHA, 1 . ) ) ,3HSIN, 
LSQMNL (LSQMNL(1 .  ,4UBF.TA,  1 .  )) ) 
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[[useri  iO.S)]]  LSOMNL(LSnM.\'L(l.  ,LIST(9)  .LSnMNL(LSQMNL(S.)) , 

SHUSERl,!.)) 

MANY(Zts,  quani,  quan2,    •■-,  quan-,)  ,   k   =  1,  2,  ...,  59, 
appends  quarii ,  quan2,    ■■■,  quan.,    in  turn,  to  the  bottom  of  lis, 
returning  lis   as  its  value.  This  is  an  actual  change  to  the  list. 

LA  =  LSQMNLCLQ1,LQ2,LQ3,LQ4) 
is  equivalent  to 

LA  =  LIST(9) 

CALL  MANY(LA,LQ1,LQ2,LQ3,LQ4) 
and  to 

LA  =  LSQMNL(LQl) 

CALL  MANY(LA,LQ2,LQ3,L04) 
This  routine  is  especially  useful  for  stacking  items  to  be  saved 
during  recursion.   It  may  also  be  used  to  build  up  expressions.   For 
example,  the  representation  of  [[1+0(0,0:  ,7)]]  may  be  constructed  by: 

C   PUT  FIRST  TERM  IN  AN  EXPRESSIf)N 

LATl  =  LSQMNL(LSQMNL(1.)) 
C   BEGIN  THE  SEC0ND  TERM 

LATT2  =  LS0MNL(1.) 
C  F(?RM  THE  ARGUMENTS  FOR  USE  BY  THE  SECOND  TERM 

LARl  =  LIST(9) 

LAR2  =  LSQMNL(LSQMNL(1. ,1HX,3.)) 

LARS  =  LSQMNL(LSQMNL(7.)) 
C  C0MPLETE  THE  SECOND  TERM 


28 


2.5 

CALL  MANY(LATT2 , LARl , LAR2 , LARS , IHG , 1 . ) 
C  F0RM  THE  FINAL  EXPRESS I0N 
CALL  MANY (LATl , LATT2) 

(This  last  step  is  correct  only  because  the  ordering  of  terms  is 
clear  here.   In  general,  LSQADD  as  described  in  §2.6  should  be  used.) 

LSSCPY(Z-is)  returns  a  complete  copy  of  the  list  or  expression 
lis.      This  is  useful  wlien  a  routine  destroys  or  modifies  an  argument 
which  must  be  available  later. 
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Exercises 

2.1.1.  Write  code  to  generate  the  correct  representation  of  the 
following  expressions: 

[[01]  [[1/311        [[xV3ll       [[x(u)/^]] 

[[111  [[xll  [[-^Vsn      [\xiSij)]] 

[[-in         [[x/31]        [[x^/31]       [[-IxiZu)-'']] 

2.1.2.  Describe  the  list  LB  that  results  from  the  following  code: 

LA=LSQMNL(1.) 

LAA=LSQMNL(-2.5,I.SQMNL(LA)  ,I.SOMNL(L.SSCPY(I.A))  , 
*   IIIX.LA) 
LB=LIST(9) 

CALL,  MANY(LA,1IIY,1.) 
CALL  MANY  (LB,  LM) 

2.1. .3.   What  expression  does  the  list  LB  formed  above  represent? 

2.1.4.  The  list  { (1 . , 1I!X,-1 . , IHX, 1 . ))  is  not  a  proper  expression 

since  the  reduction  to  ((1.))  by  cancelling  the  powers  of  IIIX,  is  required, 

Use  tlie  expression  symbol  3111.*  as  the  identity  function  to  represent 

parenthesization,  and  write  code  for  the  expression  [[j-/(j:)]].   (Do 

not  expend  too  much  energy  worrying  about  proper  ordering.) 


30 


2.5 


2.2.   Input  of  Expressions 

The  function  INLIST  provides  a  means  of  translating  FORTRAN-like 
expressions  into  their  internal  representation  on  lists.   Input  may 
be  taken  from  a  file  or  from  line  images  saved  on  a  list. 

INLIST(yar,  quan)   returns  a  list  representing  the  expression 
found  on  the  logical  unit  or  list  quan.     The  variable  or  array 
element  var   will  also  contain  the  newly  created  expression.   Regardless 
of  whether  quan   is  a  unit  number  or  a  list,  input  is  handled  line  by 
line,  scanning  for  an  expression  in  columns  7  through  72  of  each  line 
until  either  a  $  is  encountered  or  a  line  is  encountered  in  which  any  of 
columns  1  through  5  are  non-blank.   The  end  of  the  input  stream  will 
also  terminate  the  scan.  The  following  syntax  is  applied  to  the 
characters  obtained  in  this  manner  after  removing  any  blanks. 

<expression>  ♦-  <term>  {<addition  operator><term>}o 

<term>  ■*-   <factor>  { **<factor>}o  {<multiplication  operator> 

<factor>  {**<factor>}o}o 
<factor>  -<-  {<addition  operator>}o  {<number>  |  <symbol>  | 

{<symbol>}o(  <expression>  {  ,  <expression>}o  )} 
<addition  operator>  ♦-  +  |  - 
<multiplication  operator>  -^  *  |  / 
<symbol>  ■<-  <letter>  {<letter>  |  <digit>  |  .}o 


<number>  -<-  <whole>  {<fraction>}o  {<exponent>}o  | 


,  1 
<fraction>  {<exponent>}o 


<whole>  ■*■   {<digit>}i 
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<fraction>  *■   .    {<digit>}o 

<exponent>  *-  E   {<ad<iition  operator>}o  <whole> 

<ietter>  *-a|b|c|d|e|f|c|ii|i|j|k|l 
m|n|0|p|o|r|s|t|u|v|w|x| 

Y    I    Z 
<digit>  ^0|l|2|3!4|5|6|7|8|9 


Thus  expressions  are  formed  from  symbols  and  numbers  combined  by 
the  arithmetic  operators  +,  -,  *,  /,  and  **,  balanced  parentheses 
and  commas . 

Symbols  are  strings  of  from  1  to  10  letters,  digits,  or  periods, 
which  begin  with  a  letter.   For  example,  the  following  are  valid 
symbols . 

A 

X. 1.2. ABC 

SIN 

FUNCTI0N 
The  following  are  not  valid  symbols, 
1..A 

$A=  SIN(X) 

•TRUE.         lElO 

X  +  1  *A 

Numbers  are  combinations  of  a  decimal  whole  number  with  a 
decimal  fraction  and  exponent.   Either  the  whole  number  or  the 
fraction  must  appear,  but  in  either  case  the  remaining  two  fields 
are  optional.   If  all  fields  are  present  a  number  assumes  the  form 


A..1 

LAB0R PARTY 

ETC... 

DIMENSI0N 

C0S 

TAN 

PR0GRAM 

LS0249 

LA B0R PARTIES 

.N0T. 

-B 

1.* 
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n.m   E  k 
which  is  interpreted  as 

while,  if  some  fields  are  omitted,  a  number  may  have  the  forms 

n  n.m  n   E  k  .m  .m  ?.  k 

which  are  interpreted  as 

n.  n.m  (n. )  10^  O.m       {O.m)    10^ 

where  n  is  a  string  of  one  or  more  digits,  m   is  either  empty  or 
a  string  of  digits,  and  k  is  a.   string  of  one  or  more  digits 
optionally  preceded  by  a  string  of  pluses  and  minuses  determining 
the  sign  of  the  exponent.   The  following  are  examples  of  valid 
numbers  with  the  indicated  interpretations  in  parentheses: 

1        fl.O) 
.5       (0.5) 
1.E--3   (1000.0) 
l.lF+2    (110.0) 
OEIO     (0.0) 

Symbols  are  stored  in  the  resulting  expression  directly  as 
expression  symbols;  i.e.,  the  characters  of  the  symbols  are 
used  left  justified  with  blank  fill.   Numbers  are  stored  as 
expression  numbers;  i.e.,  the  number  is  translated  into  a  single 

precision  floating  point  number,  holding  between  14  and  15  digits 

-''94 
accuracy  with  magnitudes  ranging  from  npproximntely  10     to 

322 
10    for  non-zero  numbers. 

A  symbol  may  ae  used  as  a  function  name  by  following  it  with 

an  argument  list  of  expressions  separated  by  commas  and  enclosed 


• 

(0.0) 

5E-1 

(0.5) 

1.E3 

(1000.0) 

l.lF-2 

(0.0011) 

0 

(0.0) 
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in  parentheses.   If  the  argument  list  is  present  without  a  preceding 
function  name,  the  expression  symbol  3H1.*  is  used  as  the  function 
name  in  the  resulting  expression.   For  example,  the  following 
are  valid  function  references: 

(X  +  1.1)  SIN(3.14*Y.2) 

(0,  (((2**3))) ,K)    TAN(X**(EXP,SIN) ,  +  -+. ,A/B) 

(((C0S(n))))      c0S((((n)))) 

A (3/B )   X(X(X,X),X) 

A  factor  is  formed  by  taking  a  number,  symbol,  or  function 

reference,  and  optionally  preceding  it  by  a  string  of  pluses  and 

minuses.   Each  +  is  ignored,  and  each  -  alternates  the  resulting 

sign.   A  term  is  formed  by  combining  factors  with  the  operators 

**  (exponentiation),  *  (mnltirlication) ,  and  /  (division).   An 

expression  is  formed  by  combining  terms  with  the  operators  +  and  -. 

Translation  is  done  in  a  left  to  right  scan,  with  operators  being 

applied  in  the  following  order: 

1 •    + ,  -      (unary) 
2    ** 

3.  *,  / 

4.  +,  -      (infix) 

Thus  the  following  expression  generates  intermediate  results  as 
indicated: 

A  +  B*-C/D**E*F  -  X**Y**Z**-W/L/M*N 
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A 

A,  B 

A,  B,  -C 

A,  B*(-C) 

A,  B*r-c;,  D 

A,  B*r-c;,  D,  E 

A,  B*r-c;,  D**E 


A, 
A, 
A, 

A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
A+r 
rA+ 


B*r-c;;/rD**E; 
B*r-c;;/rD**E;,  f 

B*r-c;;/rD**E;;*F 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**Ej;*F; 

B*r-c;;/rD**E;;*Fj 

B*r-c;;/rD**E;;*F; 

R*r-c;;/rD**E;;*F; 

B*r-c;;/rD**F;;*F; 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**E;;*F; 

B*r-c;;/rD**E;;*F; 


X 

X,  Y 
X**Y 
X**Y,  Z 

rx**Y;**z 
rx**Y;**z,  -w 
rrx**Y;**z;**r-w; 
rrx**Y;**z;**r-w;,  l 
rrrx**Y;**z;**r-w;;/L 
rrrx**Y;**z;**r-w;;/L,  m 
rrrrx**Y;**z;**r-w;;/L;/M 
rrrrx**Y;**z;**r-iv;;/L;/M,  n 
rrrrrx**Y;**z;**r-wj;/L;/M;*N 


rB*r-c;;/rD**E;;*F;;-rrrrrrx**Y;**z;**r-w;;/L;/M;*N; 
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If  quart   is  a  unit  number,  INLIST  will  read  one  line  at  a  time 
from  the  current  position,  beginning  with  the  next  complete  line, 
and  ending  with  the  line  which  terminates  the  input  expression.   If 
the  expression  is  terminated  by  a  $  within  columns  7  through  72,  the 
remaining  columns  may  be  used  for  a  comment.   If  the  expression  is 
terminated  by  non-blank  characters  in  columns  1  through  5  of  a  line 
following  the  expression,  the  entire  line  may  be  used  for  a  comment, 
as  long  as  at  least  one  of  columns  1  through  5  is  not  blank. 

If  queen   is  a  list,  INLIST  assumes  that  each  of  the  items  on  Quan 
is  a  list  of  Hollerith  constants,  each  such  constant  representing 
ten  consecutive  columns  of  a  line.  The  list  quan   and  each  of  its 
sublists  is  scanned  from  top  to  bottom,  in  effect  ''rewinding"  quan 
with  each  call  to  INLIST.  Thus,  only  the  first  expression  on  quan 
is  translated,  and  everything  beyond  the  terminator  of  this  expression 
may  be  used  for  a  comment. 

The  following  sets  of  line  images  and  calls  to  LSQMNL  will 
be  translated  by  INLIST  into  valid  representations  of  the  indicated 
expressions: 

[[0]]         0$  THIS  IS  A  C0MMENT  (file) 

LSQMNL (LSQMNL (lOH      0$   ))  (list) 

[[x  -  /]]     X  -  Y**Z  (file) 

*   THIS  IS  ALSiS  A  COMMENT 

LSQMNL (LS(^INL(10H      X  -  ,10HY**Z      ),  (list) 
LSQMNL (IH*)) 
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[[S/F  +  Sx/y  -    {X  *  y)^"^]] 

3*Z**.5  +  (file) 

*   5*X   /   Y 
-(X  -  -Y)**+(3*R) 
SEND  0F  EXPRESSI0N 

LS0MNL(LSQMNL(1I'  , 7H3*Z** . 5 , ]H+ , 3H5*X)  ,  (list) 
LS0MNL(1I!  ,4II/Y-(), 
LSQMNL(1H  ,10HX  -  -Y)**+ , 5H(3*R) ) ) 

There  are  other  possible  forms  of  input  for  INLIST.  These  are 
discussed  in  detail  in  §3.2.  One,  which  insures  compatibility  with 
some  of  the  forms  of  expression  output  provided  by  the  routine  LSQPNT 
(§2.3),  will  be  mentioned  here.  If,  instead  of  being  terminated, 
an  expression  is  followed  by  an  =,  the  expression  is  discarded,  and 
the  scan  for  an  expression  begins  anew.  Thus,  an  input  line  of  the 
form 

LA(X.Y,  Z)=0=5$ 
will  be  translated  into  a  valid  representation  of  [[5]].   Though 
the  expressions  before  tlie  last  equals  sign  are  effectively 
comments,  they  must  conform  to  the  syntax  for  input  expressions. 
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Exercises 

2.2.1.  Which  of  the  following  are  valid  symbols  for  INLIST? 
.A  FUNCTIiSN       A+B/C  UIMENSI0N 

A.  1.*  E5  DIMHNSI0NED 

A.1F.  +  5  (A )    1.F.5  A 

2.2.2.  Which  of  the  following  are  valid  numbers  for  INLIST? 

1.0  1234567890.0987654321 

.0        l.lE-5         n-5 
.1        3.14159r.--5    .1U3 

2.2.3.  Consider  the  following  input  stream  to  lie  on  a  file. 
Assume  INLIST  to  be  called  five  times  for  that  file.   What  lists 
would  be  formed?  What  if  the  file  were  a  list? 

LA=IF(A.EO.B)+SIN(5.E-l)/3**2SXYZ$$ 

X ♦¥ (X *  *3)  *  *X .  +X 

$ 

H01V=N0W=  BR0WN=C0W ( BULL)  $  MY 

♦G00D/NESS.0R.A 

$ 

FUNC(ARr,l,ARG2. 

SIN(ARG3)**3*X, 

ARG4)*C((((0))))) 
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2.3.  Output  of  Expressions 

The  functions  LSQPNT  and  LSQ0UT  provide  a  means  of  translating 
expressions  from  their  internal  representation  as  lists  to  a 
FORTRAN  like  notation.   LSQPNT  performs  the  translation  to  the  level 
of  Hollerith  constants,  and  I,SO0UT  formats  these  into  line  images 
on  a  file  or  a  list. 

LSQ0UT(4HUNIT,  quan)   directs  all  further  output  to  the 
logical  unit  quan.      For  example, 

CALL  LSQ0UT(4HUNIT,  6L0UTPUT) 
would  direct  lines  to  the  file  0UTPUT.   LSQ0UT  assumes  this  unit  as 
a  default. 

LSQ0UT(4HLIST,  lis)   directs  all  further  output  to  the  list 
lis.      As  each  line  is  completed  it  is  broken  into  Hollerith  constants 
representing  consecutive  groups  of  ten  columns.  Tliese  constants  are 
then  placed  on  a  list,  with  columns  1  through  10  at  the  top.   This 
newly  created  list  is  placed  on  the  bottom  of  lis. 

LSOPNT(HSj  hoi)    translates  the  expression  lis   into  line 
images.   The  first  line  begins  with 

hoi   = 
in  column  7,  and  the  last  two  lines  are  of  the  form 
$ 

$nND  0F  EXPRCSSI0N 
with  the  $  in  column  2.   Beginning  after  the  eciuals  sign  on  the 
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first  line,  columns  7  through  72  of  the  line  images  are  used  for 
a  FORTRAN  like  representation  of  lis,    and  possibly  of  some  of  its 
subexpressions,  conforming  to  the  syntax  given  for  INLIST  in  §2.2. 

LSQPNT(0,  5HN0SUB)  causes  LSQPNT  to  translate  all  subexpressions 
where  they  occur.   If  this  call  is  made,  then  LSQPNT  will  translate  an 
expression  into  a  single  FORTRAN  like  expression.   Otherwise,  some 
deeply  embedded  subexpressions  may  be  replaced  by  symbols  of  the  form 
N.O,  N.l,  N.2,  etc.,  with  these  subexpressions  translated  after  the 
main  expression. 

LSQPNT(0,  3I1SUR)  undoes  the  effect  of  LSQPNT(0,  5HK0SUB) . 
This  is  the  default  state  of  LSQPNT.   It  makes  output  expressions 
more  readable,  but  is  in  conflict  with  normal  INLIST  usage. 

In  most  cases,  expressions  take  less  list  space  as  line  images 
than  in  their  internal  representation.  The  user  may  take  advantace  of 
this  fact  to  compress  infrequently  used  expressions  as  follows: 

CALL  LSQPNT (O,5HN0SUB) 

LSAVE=LIST(9) 
CALL   LSQ0UT(4I!LIST,LSAVE) 
CALL   LSQPNT(LXPR,4I!LXPP) 
CALL   LSQDES(LXPR) 

LXPR=INLIST(LXPR,LSAVE) 
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The  code  first  conditioned  LSQPNT  to  translate  subexpressions  where 
they  occur.   This  need  be  done  only  once.   Then  tlie  list  LSAVE  was 
created  to  hold  the  line  images  of  the  expression  LXPIL   LSQPNT  was 
called  to  dump  LXPR  to  LSAVE.   Then  LXPR  was  destroyed.   Later  it 
was  recreated  by  a  call  to  INLIST. 

If  no  such  special  use  of  LSQPNT  has  been  made,  an  expression 
may  be  printed  simply  by  a  call  such  as 

CALL  LSQPNT(LXPR,6HMYLIST) 
but  if  such  use  has  been  made  one  should  first  use 

CALL  LSQ0UT(4HnNIT,6L0UTPlJT) 
to  insure  that  output  will  go  to  the  file  0UTPUT. 
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Exercise 

2,3.1.  Run  the  following  program  on  the  data  in  exercise  2.2.3.  Add 
to  the  data  such  expressions  as  may  interest  you.  Compare  the  output 
with  the  input.   (Set  a  short  time  limit  to  make  the  program  stop.) 

PR0GRAM  THST( INPUT, 0UTPUr) 

DIMENSION  SPACE (5000) 

CALL    INITAS(SPACE,5000) 
1  CALL    INLlST(LA,5LINPLrr) 

CALL   LSQPNT(LA,2HLA) 

G0  T0    1 

END 
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2.4.   Destruction  of  Lists  and  Expressions 

The  functions  LSQDES  and  LSQUNM  unbuild  lists  and  expressions. 

LSQDESCHs)  erases  the  list  or  expression  lis.      Lists  that 
are  no  longer  needed  must  be  erased  so  that  their  cells  will  be 
available  for  other  lists.   If  lis   is  a  sublist  of  another  list,  it 
will  not  actually  be  destroyed.   However,  the  user  should  consider 
lis   unavailable  in  all  cases,  since  a  sublist  is  erased  when  the  list 
of  which  it  is  a  sublist  is  destroyed. 

LSQUNM  (Hs,  vavi ,   var2 ,    ••••  var,),   k  =   1,  2,  ...,  59, 
removes  k   items  from  the  bottom  of  list  lis,    the  bottom  item  going 
to  var,,    the  next  to  var,    .,  etc.   For  example,  an  expensive  way  of 
doing  very  little  is: 

LA=LSQMNL(X,Y,Z,W,L,M) 

CALL  LSQUNMCLA,X,Y,Z,W,L,M) 
This  leaves  LA  as  an  empty  list.   LSQUNM  returns  lis   as  its  value. 
Thus  a  way  of  doing  effectively  nothing  is: 

CALL  LSQDES (LSQUNM (LSQMNL (A, B), A, B)) 
which  creates  a  list  with  A  and  B  on  it,  pops  them  off  and  destroys  the 
list. 
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Exercises 

2.4.1.  What  does  the  following  code  do? 

CALL  LSQDIiS(LSQUNM(LSQMNL(A,B,C),B,C,A)) 

2.4.2.  The  program  given  in  exercise  2.3.1  contains  an  error. 
Correct  it.   (Hint:  What  would  happen  if  the  data  produced  3000 
expressions?) 

2.4.3.  Which  lists  are  available  after  the  following  code  is 
executed? 

LA=LIST(9) 

LE=L1STC9) 

LC=LIST(9) 

LU=LIST(9) 

CALL  MANY  C LA , LB , LC , LU , LD , LC , LB) 

CALL  LSnnKS(LR) 

CALL  LSQDES(Ll)) 

CALL  LSQI!NM(LA,LF,LD,LC,LB) 

CALL  LSQDES(LA) 
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2.5.   Other  Output  Facilities 

The  function  PRLSTS  provides  a  mechanism  for  printing  dumps 
of  lists.   This  is  primarily  for  debugging  purposes,  when  lists 
have  been  created  to  represent  expressions,  but  cannot  be  printed 
by  LSQPNT  due  to  some  error. 

PRLSTS(Zts,  int)   prints  the  list  lis   in  the  mode  determined 
by  the  integer  quantity  int,   which  may  be  1,  2,  3,  or  4 ,  for  integer, 
Hollerith,  real,  or  octal  dumps  respectively.   Each  dump  begins  a 
new  page.  The  most  likely  call  is 

CALL  PRLSTS(LXPR,4) 
which  will  dump  LXPR  in  octal. 
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^ESIM  LIST 
?ESIv  Sl='L!S"' 

E'w-  SJ^-IST 
'^E'^I'M  SL^L!5T 

^E3Iv  S-'LlS" 

175i4Gncia'?Gf"33ai3Q3(J 

3"=?5?555?5t5'^^559555 

3l«=5  =  5  =  5^55555555355 

32  =  r  =  i'?5  =  5S555555555 


r-:cic-e  5.5.:.      >rtrut   frtw  a   call  to  PRL5TSaiS.4T   with  LIS 

reareseitijtg  the  exiJTession   [[1  ♦  r  -  r-   •  r^z]] 
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2.6.   Arithmetic  Routines 

The  functions  LSQADD,  LS(?-1EX,  and  LSQRAZ  are  available  to  perform 
arithmetic  opcratior.s  on  pairs  of  expressions.   Each  returns  an 
expression  as  its  value,  and  does  not  change  its  arguments. 

LSQAUD[Z-isi ,  lis2)    returns  the  sum  of  expressions  lisi 
and  lis2 

LSQMEX(Zisi  ,  lis2)    returns  the  product  of  expressions 
lis  I    and  lis2- 

LSQRAZ(lisi,    lis2)    returns  the  result  of  raising  the 
expression  lisi    to  the  expression  lis2   power. 

In  order  to  perform  subtraction  or  division,  addition  of  the 
negative  or  multiplication  by  the  inverse  respectively  must  be  used. 
For  example: 

LA.MI.\S=LSQMNL  (LSQMNL  (- 1 . )  ) 
LD;  =  LSQMEX  CLAMINS, LC) 
LCI  =  LSQRAZ(LC,LA.MIXS) 
LDD=LSQADD(LB,LCM) 
LDQ=LSQ'-1EX(LB,LCI) 
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leaves  LB-LC  in  LDU  and  LB/LC  in  LUQ. 

Neither  LSQADD  nor  LSQMEX  introduce  parentheses.   The  results 
are  obtained  by  combining  the  terirs  of  the  operands  to  form  new 
terms.   For  example,  in  the  following,  LDD  will  be  an  empty  expression, 
since  all  its  terms  will  cancel. 

LAMINS=LSQMNL(LSQMNL(-1 . )) 

LCM=LSQMEX(LAMINS,LC) 

LDD=LSQAUD(LC,LCM) 

LSQRAZ  does  introduce  parentheses  in  many  cases.   A  sum  of  terms 
raised  to  some  non-trivial  power  will  not  be  expanded.   However  a 
single  term  raised  to  a  power  will  have  that  power  distributed  over 
its  factors.   Thus,  in  the  following  example,  LDQ  will  be  the  list 
((1.)),  representing  [[1]],  if  LC  is  a  single  term,  but  not  necessarily 
in  any  other  case. 

lamins=ls0mnl(lsqmnl(-1.)) 
lci = lsqraz (lc, lam ins) 
ldq=lsqml;x(lc,lci) 
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Exercises 

2.6.1.  Write  code  to  p.enerate  the  following  arithmetic  operations 
on  expressions  LA,  LB,  LC,  LD,  etc. 

LA+LB  LA+LB*LC       LA*LB+LC*LD+LR*LF 

LA-LB         LA-LB**LC      LA**2+LB**2+LC**2 
LA*LB         3*LA  LA/(1-LB**2)**.5 

LA/LB  LA**3  LA*LB*LC*LD*LF*LG*LH 

2.6.2.  Devine  the  purpose  of  the  following  function 

FUNCTI0N  LSOSVC(LA,N) 
DIMENSI0N  LA(N) 
LSQSVC=LA(1) 
D0  100  J=2,N 
LB=LSQSVC 

LSQSVC=  LSQAOU ( LSQSVC , LA ( J ) ) 
IF(J.GT.2)CALL  LSQDF.S(LB) 
100   C0NTINUE 
RETURN 
END 

2.6.3.  Write  a  subroutine  LSQMMT  to  multiply  two  2  by  2  matrices. 
-  two  3  by  3  matrices.  -  two  N  by  N  matrices.-  an  N  by  M  matrix  by 
an  M  by  K  matrix. 
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2.7.   Definition  and  Evaluation 

The  functions  LSQVAL  and  LSQDEF  provide  facilities  for  evaluating 
expressions  in  terms  of  user  provided  definitions  of  expression  symbols. 
The  term  "evaluation"  refers  to  the  application  of  all  currently  known 
definitions  to  an  expression.   It  is  proper  to  leave  symbols  undefined. 
Thus  this  process  always  returns  an  expression,  not  a  FORTRAN  constant. 
The  routine  LSQTYP  (see  §2.10)  may  be  used  to  obtain  a  FORTRAN  constant 
from  a  constant  expression. 

As  part  of  standard  program  initialization  (see  LSQINI  and 
INITAS  in  §2.11)  certain  definitions  are  made.   If  these  and 
subsequent  user  definitions  are  suspended,  temporary  definitions 
may  be  introduced  so  that  evaluation  will  effect  substitution.   Such 
a  mechanism  is  provided  by  the  routine  LSQSBS.   The  calls  to  LSQDEF 
and  LSQVAL  used  by  LSQSBS  are  among  those  described  in  §3.6,  where 
the  precise  definition  and  evaluation  methods  used  are  presented. 

LSQVAL(Hs)  returns  the  expression  representing  the  value 
of  expression  lis,   obtained  by  applying  all  currently  known  definitions 
to  the  expression  symbols  of  lis.      For  example,  unless  C0S  is  re- 
defined, the  following  code  will  leave  ((1.))  in  LAV,  representing 
[[1]]. 

LA=LSQMNL(LSQMNL(1.,LIST(9),3HC0S,1.)) 

LAV= LSQVAL (LA) 
The  input  expression  lis    (e.g.  LA  above)  is  unchanged. 
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LSQDEF (sj/mfc ,  lis,   0)  defines  the  symbol  symb   to  be  the 
expression  lis.      The  expression  lis   should  be  considered  erased  by 
this  call.   Where  the  symbol  symb   is  used  as  a  function  name,  its 
arguments  will  be  applied  to  the  expression  symbols  of  lis.      For 
example 

LA=LSQMNL(LSQMNL(1. ,LSQMNL(LSQMNL(1. ,1HX,1.)), 
$   SHSIGMA.l.)) 

CALL  LSQDEF(5HSIGMA,LSQMNL(LSQMNL(1.), 
$   LSQMNL{1. ,3HTAU,1.)) ,0) 
LAV=LSQVAL(LA) 
will  leave  a  representation  of  [[1  +  t(x)]]  in  LAV. 

LSQDEF(si/?r?i>,  lisy,    lisi)   defines  the  symbol  symb   to  be  the 
expression  lisi,    in  which  the  expression  symbols  on  the  list  lis2 
hold  the  places  into  which  any  arguments  applied  to  symb   are  to  be 
substituted.   Both  lis\    and  lis2   should  be  considered  to  be  erased 
by  this  call.   The  list  of  dummy  arguments  lis2   should  be  a  list  of 
expression  symbols  not  used  in  any  other  such  call  to  LSQDEF.   In 
the  course  of  an  expression  evaluation  in  which  symb   is  encountered, 
the  arguments  to  which  it  is  applied  will  be  scanned  from  left  to 
right,  while  lis2    is  scanned  from  top  to  bottom,  temporarily  defining 
each  expression  symbol  found  on  lis2   to  be  the  corresponding  argument. 
The  value  of  lis\    subject  to  these  definitions  is  used  in  place  of 
symb   and  those  of  its  arguments  which  have  been  matched.   Any  unmatched 
actual  arguments  are  applied  to  the  result.   Though  it  is  not  an 
error  to  leave  some  dummy  arguments  unmatched,  the  user  is  advised 
to  read  §3.6  before  so  doing. 
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As  an  example  of  this  call  consider 

CALL  LSQDEF(3HSIN,LSQMNL(LSQMNL(1. , 
$   LSQMNL{LSQMNL(1.) ,LSQMNL(-1. ,LSQMNL( 
$   LSQMNL (1 . , 3HL . . , 1 . ) ) , 3HC0S, 2 . ) ) , 3H1 . * , . 5) ) , 
$   LSQ?4NL(3HL..)) 
which  defines  [[sinil..)]]    to  be  [[/(I  -  cos{l..)^]] . 

LSQSBS(ltsi ,  symb,    lis2,   -1)   returns  the  result  of 
substituting  lis\,    an  expression,  for  symb,   an  expression  symbol, 
in  the  expression  lis2-      All  appearances  of  symb,   with  or  without 
arguments  will  be  replaced.  An  infinite  loop  may  occur  if  lisi 
contains  symb.     The  expression  lis\    should  be  considered  erased 
by  this  call. 
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Exercises 

2.7.1.  Write   a  subroutine  which  when   called  will   define   3IIC0T  as 
tlie   cotangent    function,    .3I1SFX  as   the   secant    function,    and   3IICSC  as 
the   cosecant   function.      Assume   that   3MSIN   is   already  defined  as   tlie 
sine   function,    .31IC0S  as  the  cosine   function,   and   3IITAN  as   the 
tangent    function.      Call   the  subroutine   LSQTRG. 

2.7.2.  Discover  the  purpose  of  tlie   following  program. 

PR0GRAM  TnST(INPUT,0irrPUT) 
DIMr;NS10N  SPACR(5000) 
CALL    INITAS(SPACr:,5000) 

LI)l-TR=LSQMNL(LSQMNL (1 .  ,  3HM  1,1., 3HL22  , 1 .  )  , 
*        LSQMNL(-1.  ,3HL12,1.  ,3IIL21,1.)) 
CALL   LSODnFC3IILll,INLIST(LA,5LINPHT),0) 
CALL   LSQL)EF(3HL12.INLt,ST(LA,5LINPlIT),0) 
CALL   LSOnEf'(3HL21,!NLIST(LA,5LrNPIJT)  ,0) 
CALL   LSQDLr(3IIL22,INLIST(I,A,5F, INPUT), 0) 
LV[)ETR=  LSnVAL  ( LDLTR ) 
CALL   LSQDES(LDF,TR) 
CALL   LSQPNT(LVDETR,6HUETFRM) 
CALL   EXIT 
END 
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2.8.   Truncation 

The  function  LSQTRC  provides  a  mechanism  whereby  expressions 
may  be  formed  in  which  some  of  the  powers  of  a  particular  variable 
have  been  removed.  The  usual  application  is  to  power  series,  say 
[[oq   ■••  cix  *   C2X^  +  ...  ♦  ax     ♦  ...]],  in  which  the  variable  x  is 
sufficiently  small  to  warrant  the  assumption  that  powers  beyond, 
say,  the  kth   may  be  treated  as  zero.  Then  one  would  use  the  polynomial 
[[cg  +  Oix   +  C2X^  •••...+  a,x   J]  in  place  of  the  power  series.   In 
this  section  we  restrict  our  attention  to  this  simple  case  of  polynomials 
in  the  variable  on  which  trucation  is  to  be  performed.   For  the  action 
of  LSQTRC  on  more  complex  expressions  see  §3.7. 

LSQTRC (Zis,  symb,   quan)    returns  an  expression  in  which 
those  terms  of  lis,    a  polynomial  in  the  expression  symbol  symb,   which 
contain  symb   raised  only  to  a  constant  power  between  zero  and  quan, 
inclusive,  are  retained.   For  example 

LA=LSQMNL(LSQMNL(.5),LSQMNL(3.,1HX,21.), 
$   LSQMNL(-6.  .1I1X,44.,1HY,1.)) 

LA1=  LSQTRC ( LA, IHX.O.) 

LA2=  LSQTRC(LA,1HY,0.) 

LA3=  LSQTRC(LA,1IIZ,0.) 

LA4=  LSQTRC  (LA,  llIX,  21.) 

will  leave  representations  of  [  [  .5]  ]  ,  [  [  .5  ♦  3x   ]],  [[.^  +3x   -  6x  tj]]  , 

21 
and  [[.5  +  3x   ]]  in  LAI,  LA2,  LA3,  and  LA4  respectively. 
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LSQTRC(Zis,  sywb ,  quarii ,  qucm2)   returns  an  expression  in 
which  those  terms  of  lis,   a  polynomial  in  the  expression  symbol  symb , 
which  either  contain  symb   raised  to  a  constant  power  between  qucm\ 
and  quan2   inclusive,  or,  if  zero  lies  in  that  range,  do  not  contain 
symb   at  all,  are  retained.  The  order  of  quan\    and  quarii    is  not 
material,  i.e.,  LS0TRC{LA,1HX,3. ,5. )  returns  the  same  expression  as 
LSQTRC(LA,1HX,5.,3.). 
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Exercises 

2.8.1.  The  following  function  will  return  the  coefficient  of  the 
expression  symbol  SYMB  raised  to  the  power  QUAN  in  the  polynomial 
LIS.   Use  this  function  to  write  a  program  wliich  will  differentiate 
polynomials  of  known  degree. 

FIINCTI0N  LSO('C0(SYMB,QUAN,LIS) 

LA=LSQTRC( LIS, SYMB, QUAN, QUAN) 

LB=LSQMNL(LSQMNL(1.)) 

LSOGC0=LSQSBS (LB , SYMB , LA , - 1 ] 

CALL  LSQDESCLA) 

RETURN 

END 


2.8.2.   Write  functions  LSOADT  and  LSOMIIT  which  will  add  and  multiply 
a  pair  of  polynomials  while  truncating  on  some  expression  symbol  to 
a  specified  {lower. 
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2.9.   Differentiation 

The  functions  LDIF  and  LOOKUP  provide  a  mechanism  for  taking 
the  derivatives  of  expressions  with  respect  to  expression  symbols 
using  user  provided  definitions  of  the  partial  derivatives  of 
functions.   By  differentiation  we  mean  the  application  of  the 
following  rules,  for  a  particular  expression  symbol  symb: 

D       Aa^)  -*  D       Aa)ba  la  *  D       Ab)a^loaia) 
symh  symb  symb  ' 


^syrnb^^k^'^^^^^^'^^^^^''  '^2 e^)  ,k) 

D       A  symb')   -*■         1,  if  expression  symbol  symb'   is.  symb 

0,  if  expression  symbol  symb'    is  not  symb 
D       ,  icons)   -*  0,    for  any  constant  expression  cons 


where  the  function  partial   is  determined  by  user  definitions  via 
LOOKUP.   It  is  proper  to  differentiate  an  expression  involving 
functions  for  which  some  or  all  partial  derivatives  are  not  defined. 
In  such  a  case,  the  expression  symbol  7IIPARTIAL  will  be  used  as  the 
function  name  for  partial. 

LDIF(Zis,  symb)    returns  the  derivative  of  expression  lis 
with  respect  to  expression  symbol  symb.      The  resulting  expression 
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is  the  derivative  in  the  sense  of  the  above  rules,  in  effect  a 
partial  derivative  with  respect  to  aymb ,   since  only  explicit 
functional  dependence  is  considered.   For  example,  if  LA  represents 

[[xy   ♦  logi-x)    *  xOj)]] 
then  LDIF(LA,1HX)  will  return  a  representation  of 

[[y   '   IM] 
since  x{y)   does  not  depend  on  x. 

L00KUP(O,  symb,    lis,   int)   defines  the  expression  lis   to 
be  the  partial  derivative  of  the  function  represented  by  the 
expression  symbol  symb,   with  respect  to  the  argument  int,    an  integer 
quantity.   Arguments  are  numbered  from  the  left,  beginning  with  1, 
Where  an  actual  argument  of  symb   is  to  be  used  in  the  partial 
derivative,  the  function  ARG.(fe)  should  be  used  in  lis,   with  k 
as  the  number  of  the  argument.  The  expression  lis   should  be  considered 
erased  by  this  call.   For  example,  the  partial  derivatives  defined 
in  standard  initialization,  see  §2.11,  include  the  following. 

CALL  L00KUP(O,3II1.*,LSQMKL(LSQMNL(1.))  ,1) 
CALL  L00Kl)P(O,3HSIN,LSQMNL(LSQMNL(l.  , 

$   LSQMNL(LSQMNL(1. ,LSQMNL(LSQMNL(1 . ) ) ,4HARG. ,1.)), 

$   3HC0S,l.)).l) 

This  defines  the  obvious  derivatives  of  the  identity  function  and 
the  sine. 
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Ezeraises 

2.9.1.   Apply  the  following  program  to  the  indicated  expressions  as 
data.   Differentiate  those  expressions  with  respect  to  X  by  hand 
and  compare  with  the  output.   (Use  a  short  time  limit.) 

PR0GRAM  TEST( INPUT, 011TPUT) 

DIMENSI0N  SPACF.(5000) 

CALL  INITAS(SPACE,5nOO) 
1    CALL  INLIST(LA,5LINPUT) 

LB=LUIF(LA,1IIX) 

CALL  LSQDES(LA) 

CALL  LSnPNT(LB,5HDERIV) 

CALL  LSQDES(LB) 

G0  T0  1 

END 
data  record: 

1$ 

X$ 

X**2$ 

X**N$ 

1/X**N$ 

SIN(A*X)/C0S(A*X)$ 

F(X,X,X,X.X,X,X,X,X,X)$ 

F(X,Y,Z,W)$ 

L0G(N*X)-N*L0G(X)$ 

($$  BAD  INPUT  TO  TERMINATE  RUN 
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2.9.2.  Write  a  subroutine  I.SOGRD  which  will  calculate  the  gradient 
of  an  expression  with  respect  to  an  array  of  expression  symbols. 

2.9.3.  The  following  function  returns  the  Nth  derivative  of  expression 
LIS  with  respect  to  expression  symbol  SYMR.   Use  it  to  write  a  program 
which  calculates  the  first  ten  terms  of  the  Taylor  scries  of  an  expression 
about  the  origin,  and  form  the  indefinite  integral  of  those  terms. 

FUNCTIPN  I,SQNDR(LIS,SYMB.N) 
ir(N.GT.O)  G0   T0  1 
I,SONDR=LSSCPY(LI.S) 
RETURN 

1  LN=N 
LYS=LrS 

2  LSOM)R=L[)IF(LYS,SYMB) 
IF(LN.NE.N)CALL   LSQDES(LYS) 
I  F(lJ^.EQ.l)  RETURN 
LN=LN-1 

LYS=LSONDR 
GO  T0  2 
END 
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2.10.   Analysis  of  Lists  and  Expressions 

The  functions  SLQRDR,  SF.QLR,  and  SFIQLI.  provide  a  means  of 
scanning  lists  and  expressions  element  by  element.   Lists,  expressions, 
and  their  elements  may  be  further  analyzed  by  use  of  the  functions 
LSQTYP,  LI  STMT,  LSQCMP,  LSTF.QL,  LSQCNM,  and  LQUAL. 

SEQRDR(Zte)  returns  a  quantity  to  be  used  as  a  "sequence 
reader"  for  the  list  or  expression  lis.      A  sequence  reader  contains 
sufficient  information  to  allow  SF.QLR  to  find  the  next  element  to 
the  right  on  the  list,  and  for  SF.QLL  to  find  the  next  element  to  the 
left  on  the  list.  As  initially  provided  by  SF.QRDR,  a  sequence  reader 
is  at  the  head  of  the  list,  pointing  to  the  right  to  the  top  element 
of  the  list,  and  to  the  left  to  the  bottom  element  of  the  list.   The 
functions  SEQLR  and  SEQLL  "advance"  readers,  so  that  each  call  allows 
elements  further  into  the  list  to  be  seen. 

SF.QLR(yar,  intvar)   returns  the  next  element  to  the  right 
(downwards)  on  the  list  or  expression  for  which  the  variable  or 
array  element  var   contains  a  sequence  reader.  The  integer  variable 
or  array  element  intvar     is  set  to  -1,  0,  or  +1  if  the  element 
returned  is  a  quantity  other  than  a  sublist,  a  sublist,  or  the  head 
of  the  list  respectively.   In  addition,  var   is  advanced  to  the 
right,  so  that  it  now  points  one  element  further  to  the  right. 

SEQLL(yar,  intvar)    returns  the  next  element  to  the  left 
(upwards)  on  the  list  or  expression  for  which  the  variable  or  array 
element  var   contains  a  sequence  reader.  The  integer  variable  or 
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array  element  intvar   is  set  to  -1,  0,  or  +1  if  the  element  returned 
is  a  quantity  other  than  a  sublist,  a  sublist,  or  the  head  of  the 
list  respectively.   In  addition,  var   is  advanced  to  the  left,  so  that 
it  now  points  one  element  further  to  the  left. 

SEQLR  and  SEQLL  may  be  used  in  any  combination.   It  is  proper 
to  advance  more  than  once  around  a  list.   It  should  be  borne  in  mind 
that  the  head  of  a  list  is  treated  as  an  element,  which  for  most 
practical  purposes  is  an  undefined  quantity.  Consider,  as  an 
example,  the  following  code: 

SA=SEQRDR(LA) 
100  SAR=Sr:QLR(SA,NR) 

IF(NR.NE.O)  G0  T0  100 

where  LA  is  a  list.   If  LA  contains  any  sublists  then  this  code 
will  leave  SAR  set  to  the  first  sublist,  but  will  loop  forever 
if  LA  does  not  contain  any  sublists.  On  the  other  hand,  the 
following  code  will  terminate  with  a  count  of  the  sublists  of  list 
LA  in  NSLA,  and  a  count  of  all  other  elements  in  MELA.   Note 
that  sublists  of  sublists  are  not  counted,  nor  are  other  elements 
of  sublists. 

SA=SEQRDR(LA) 
NSLA=0 
MELA=0 
100   SAR=SEQLR(SA,NR) 
IF(NR)101,102,103 
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101  mela=m[;la+i 

C0  T0    100 

102  NSLA=NSLA+1 
G0  T0    100 

103 

Thus,  if  LA  were  an  expression,  NSLA  would  contain  the  number  of 
terms  in  LA  and  MLLA  would  be  zero. 

LSQ'HP iquco-i,    var)    returns  -1  if  the  quantity  quan   is  not 
a  list  (hence  not  an  expression),  and  zero  or  preater  if  the  quantity 
quan   is  a  list  or  expression.  The  variable  or  array  element  var   is 
set  to  the  quantity  quan   if  quan   is  not  a  list,  otherwise  it  is  used 
to  provide  additional  information  on  the  nature  of  quan   considered 
as  an  expression.   If  quan   is  to  be  considered  a  list,  but  not  an 
expression  necessarily,  then  yar  will  not  contain  useful  information, 
and  only  the  sign  of  the  result  returned  by  LSQTYP  is  significant. 

If  quan   is  a  constant  expression,  then  LSQTYP  is  zero,  and  var 
is  set  to  the  FORTRAN  constant  equivalent  to  quan    (as  a  floating 
point  number).   As  indicated  in  §2.7,  this  is  the  method  by  which 
numeric  values  for  use  in  FORTRAN  arithmetic  expressions  may  be 
obtained  from  the  results  of  application  of  LSQVAL  for  evaluation. 
For  example,  a  rather  expensive  way  of  calculating  the  square  root 
of  2  is  given  in  the  following  code: 

LA=LSQMNL  ( LSQMNL  ( 1 .  ,  1  i IX  ,  .  5 ) ) 

CALL  LSQUUF(1HX,LSQMNL(LSQMNL(2.)) ,0) 

63 


2.38 

LB=LSQVAL(LA) 
NLA=LSQTYP(LB,XLA) 

This  leaves  XLA  set  to  2"  ,  and  NLA  set  to  zero. 

If  quan   is  an  expression  of  just  one  term,  but  is  not  constant, 
then  LSQTYP  will  return  a  value  ranging  from  10000  through  79999. 
If  quan   is  an  expression  of  more  than  one  term,  say  n   terms,  n 
greater  than  1,  then  LSQTYP  returns  80000+n,  and  var   is  set  to  0. 

If  quan   is  an  expression  of  just  one  term,  but  more  than  one 
factor,  say  n   factors  (exclusive  of  the  coefficient),  n   greater 
than  1,  then  LSQTYP  returns  either  SOOOO+n  or  70000+n;  the  former 
when  the  coefficient  is  1.,  the  latter  otherwise. 

If  quan   is  an  expression  of  just  one  term  which  consists  of 
a  single  factor  (possibly  with  a  non-trivial  coefficient),  where 
the  expression  may  be  considered  to  represent  [ [a  i  ]]  or 
[[a  b{,e\,e2," .  ,e,)    ]]  with  a  as  a  numeric  coefficient,  then 
(taking  k   to  be  zero  when  no  arguments  are  present)  LSQTYP  returns 
lOOOO-*-^,  20000+fe,  30000+k,  or  40000+?;;  the  first  when  a   and  c   are 
both  1.,  the  second  when  a  is  1.  and  c   is  not,  the  third  when  c 
is  1.  and  a   is  not,  and  the  last  when  neither  a   nor  c  is  1. 

In  any  case  where  quan   is  an  expression  of  just  one  term,  var 
is  set  to  the  expression  symbol  used  as  the  variable  or  function 
name  in  the  leftmost  factor. 

For  example,  in  the  following  code  NLA  will  be  10000,  and 
NAMLA  will  be  IIIX,  since  LA  represents  [[x]]. 
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LA=LSQMNL(LSQMNL(1. ,1HX,1.)) 
NLA=LSQTYP(LA,NAMLA) 

LISTMT(Zis)  returns  zero  if  the  list  lis   is  empty  and  thus 
is  the  expression  representing  [[0]].   Otherwise  a  non-zero  value 
(actually  -1)  is  returned. 

LSQCMP(Hsi,  lis2,    3HEXP)  returns  zero  if  expressions 
lis\    and  lis2   «ire  identical,  -1  if  lis\    is  lexicographically  less 
than  lis2,   and  +1  if  lis2    is  lexicographically  less  than  Ztsj.  This 
ordering  is  defined  in  detail  in  §1.4.  There  are  other  calls  to 
LSQCMP  which  may  be  used  on  fragments  of  expresssions  such  as  terms 
and  factors.  These  are  covered  in  §3.8.   For  the  moment  it  should 
be  noted  that  two  expressions  may  be  equivalent  by  application  of 
the  distributive  laws,  but  not  considered  identical  by  LSQCMP,  since 
parentheses  are  significant.  The  major  use  for  this  call  to  LSQCMP 
is  to  provide  a  means  of  sorting  expressions  to  avoid  redundant 
computations. 

LSTEQL(Ztsi,  ti32)    returns  zero  if  the  lists  (or  expressions) 
lis\    and  lis2   are  identical,  and  a  non-zero  quantity  otherwise  (again 
actually  M).   Unlike  LSQCMP,  LSTEQL  may  be  applied  to  arbitrary 
lists,  not  just  expressions,  but  does  not  provide  an  ordering  for  sorts. 

LSQCNM(.s;ymt>i  ,  symb2)    returns  zero  if  the  expression  symbols 
S}jrnh\    and  symb2   are  identical,  -1  if  symbi    is  lexicographically  less 
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than  symh2,    ♦!  if  symb 2   is  lexicographically  less  than  symbi.      For 
eaxample,  LSQCNM(n!A,lHB) ,  LSQCNM(1HA,2I1AB) ,  LSQCNM(2HAA,  3HAAB) ,  and 
LSQCNM(1IIA,211AA)  are  all  -1,  while  LS0CNM(4HAABA,4HAAAZ)  is  +1. 
LSQCNM  is  used  to  determine  ordering  of  expressions  throughout  the 
package,  and  may  be  replaced  by  a  user  routine  of  similar  effect 
if  a  different  ordering  is  desired. 

l.QUM(quani  ,  quan2)    returns  zero  if  quantities  quan\    and 
quan2   are  identical  as  bit  patterns  in  a  machine  word,  and  a  non-zero 
quantity  otherwise.   F.QUAL  may  be  used  on  any  quantities,  including 
expression  symbols,  but  will  not  provide  an  ordering  as  does  LSQCNM. 
It  should  be  noted  that  if  the  arguments  to  EQUAL  are  lists  or 
expressions,  the  value  returned  will  be  zero  only  if  the  arguments 
are  precisely  the  same  list.   Identical,  but  distinct  lists  would 
return  a  non-zero  value.  The  value  returned  by  LQUAL  should  only 
be  used  directly  in  an  arithmetic  IF  statement,  or  simple  logical 
tests,  not  as  a  floating  point  quantity  in  arithmetic  statements. 
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Exercises 

2.10.1.  The  following  function  returns  the  Nth  element  from  the 
top  of  list  LIS.   Write  a  companion  function  LSQNF.R  which  will 
return  the  Nth  element  from  the  bottom  of  list  LIS. 

FUNCTION  LSQNF.T(LIS,N) 
SA=SF.QRnR(LIS) 
00  100  J=1,N 
100  CALL  STRDIR(SFQLR(SA,M),LSONF,T) 
RFTURN 
HND 

2.10.2.  Suppose  LA  is  an  expression  and  MLA=LSOTYP(LA,XLA) .   Describe 
LA  for  the  following  values  of  MI,A  and  XLA. 

10000, .^Hl.*    20005, SIIFUNCT  30000, IIIX      40001, 3MC0S 
0,.%.5         -1,0  80120,0        30001, ^IITAN 

10001, 3IISIN    0,0  50017, 3H1.*    20001  ,4HTANII 

2.10.3.  Revise  the  programs  in  excerises  2.3.1  and  2.9.1  to  cease 
input  and  call  HXIT  on  the  reading  of  a  zero  expression.  -  on  the 
reading  of  an  expression  representing  [  [encf]  ] . 

2.10.4.  Write  a  function  LSQDGR  which  determines  the  degree  of  a 
polynomial  in  expression  symbol  SYMB  by  differentiating  to  an 
expression  which  evaluates  to  zero.  -  by  searching  for  the  highest 
power  of  SWB   in  the  polynomial. 
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2.11.   Initialization 

The  user  has  one  initialization  responsibility  to  SYMB0LANn.   I!e 
must  declare  an  array  to  be  the  list  of  available  space  from  which 
list  elements  are  to  be  drauTi  l)y  a  call  to  the  routine  INITAS.   The 
routines  of  the  package  arc  in  all  other  respects  self-initializing. 
(There  is  an  exception  for  special  print  definitions  noted  in  53.3). 

The  two  most  important  examples  of  this  self-initialization 
arise  in  evaluation  and  differentiation,  where,  on  the  first  use 
of  LSQVAL  and  LDIF,  standard  defintions  are  created  by  calling  the 
routines  LSQINI  and  LSQIDR.   The  user  may  also  call  these  routines 
to  restore  standard  definitions  after  he  has  changed  some. 

INITAS (arra!/,  int)   establishes  the  array  away   of  dimension 
specified  by  the  even  integer  quantity  irt   as  the  available  space  list, 
This  call  must  be  made  before  any  list  processing  is  done,  including 
calls  to  L00KIJP,  ISQPEF,  LSQINI,  or  LSQIHR.   This  call  should  be  made 
only  once,  since,  as  part  of  self-initialization,  some  routines 
establish  lists  to  be  used  in  all  subsequent  calls.   These  lists 
would  be  destroyed  by  a  second  call  to  INITAS.   Very  little  can  be 
done  with  an  available  space  list  of  dimension  less  than  3500,  and 
5000  is  more  practical.   Thus  a  typical  program  should  begin 

PR0GRAM  TYP(INPnT,0llTPirr) 
niMF.NSI0N  SPACh(5000) 
CALL  INITAS (SPACE, 5000) 
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LSOINI  is  called  by  LSQVAL  on  its  first  use  to  define  the 
expression  symbols  3HSIN,  3I1C0S,  3Mr.Xr,  5HL0C,  3IITAN,  fjllARCTAN, 
4IiT,V\lI,  3tll.*,  2HV.,  21in.,  and  3HTF.  .   If  the  user  wishes  to 
redefine  these  expression  symbols,  it  is  pointless  to  do  so  before 
performing  at  least  one  evaluation.   At  any  time  the  user  may  call 
LSQIM  to  reestablish  the  following  definitions: 

A.  Simple  functions  of  one  evaluated  argument 
symbol         use  meaning 

the  sine  of  X  radians 

the  cosine  of  X  radians 

X 
the  exponential  of  X,  i.e.  e' 

the  natural  logarithm  of  X,  i.e.  tog    (X) 

the  tangent  of  X  radians 

6HARCTAN  ARrTAN(X)  the  arctangent  of  X  given  in  radians 

4IITANH    TANH(X)   the  hyperbolic  tangent  of  X  railians 

B.  Special  functions 
symbol         use  meaning 

3111.*     (X)       the  expression  symbol  3111.*  is  the  internal 

representation  of  parentheses.   When  used  with 
one  argument,  the  value  of  that  argument  is  returned. 
(X,Y,...)  Ulien  used  with  more  than  one  argument,  the  evaluates 
of  those  arguments  are  themselves  enclosed  in 
parentheses  to  form  the  value. 


31  IS  IN 

SIN(X) 

3IIC0S 

C0S(X) 

3imXP 

tXP(X) 

31IL0n 

I,0G(X) 

3HTAiN 

TAN(X) 
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2IIV.  V.  (X)  the  expression  symbol  2IIV.  causes  its  argument 
to  be  evaluated  twice,  and  returns  this  double 
evaluate  as  its  value. 

2I1Q.      Q.  (X)     the  expression  symbol  2I!0.  returns  its  argument 
unevaluatod.  Thus  V.(0.(X))  will  return  the 
value  of  X. 

TilllF .  IF.(X,Y,...)  the  expression  symbol  3UIF.  provides  a 

conditional  expression  facility.   It  takes  its 
arguments  unevaluated  and,  scanning  from  left 
to  right,  returns  the  first  even  numbered  argument 
following  an  odd  numbered  argument  which  evaluates 
to  zero.   The  value  returned  is  an  unevaluated 
argument.   If  there  is  an  even  number  of  arguments 
and  no  odd  numbered  argument  evaluates  to  zero, 
zero  is  returned.   If  there  is  an  odd  number  of 
arguments,  the  last  argument  is  not  tested;  rather 
it  is  returned  unevaluated  as  the  value  of  3HIF. 
if  no  other  odd  numbered  argument  evaluates  to  zero. 
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LSOIDR   is   called  by   LDIF  to  define  derivatives   of  the 
expression  symbols    3H1.*,    3IISIN,    3HC0S,    3I1TAN,    3l!r?)T,    SIISEC,    3IICSC, 
6HARCSIN,    61IARCC0S,    6HARrTAN,    6HARCC0T,    6HARCSF.C,    f)HARCCSC,    4HSINH, 
4HC0SH,    4I1TANH,    41IC0TII,    41ISECH,    4IICSCH,    7HARCSINII,    7HARCC0SII, 
7HARCTA.NII,    7HARCC0TI1,    7HARCSECH,    7HARCCSC1I,    3I1EXP.    3I1L0C.      This   call 
is  made   on   the    first   use  of  LDIF.      The  user  may  call    LSQIDR  to 
reestablish   these  definitions   at    any   time. 

symbol  use  meanina  ar.d  derivative  defined 

3H1.*  (X)  parentheses  1. 

311SIN  SIN(X)  sine  C0S(X) 

3I1C0S  C0S(X)  cosine  -SIN(X) 

3HTAN  TAN(X)  tangent  SEC(X)**2 

3IIC0T  C0T(X)  cotangent  -CSC(X)**2 

3HSFC  SEC(X)  secant  TAN(X) *SFX(X) 

3HCSC  CSCCX)  cosecant  -C0T(X)*CSC(X) 

6HARCSIN  ARCSIN(X)  arcsine  1/ (1-X**2)**. 5 

6HARCC0S  ARa:0S(X)  arccosine  -1/(1-X**2)**.5 

611ARCTAN  ARCTAN(X)  arctangent  1/(1  +  X**2) 

6HARCC0T  ARCC0TCX)  arccotangent  -1/(1+X**2) 

6IIARCSEC  ARCSEC(X)  arcsecant  X**-l    *    (X**2-l)**- .S 

6HARCCSC  ARCCSC(X)  arccosecant  -X**-l    *    CX**2-1)**- .5 

4HSINH  SIN!1(X)  hyperbolic  sine  C0SII(X) 

4HC0SH  C0SIICX)  hyperbolic   cosine  ST\'II(X) 

41ITANH  TANIi(X)  hyperbolic   tangent  SECH(X)**2 

41tC0TH  C0TH(X)  hyperbolic   cotangent    -CSCH(X)**2 
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4HSF.CH  sr:CH(X)  hN-perbolic   secant  -SECFi(X)*TANFI(X) 

4HCSCII  CSCIl(X)  hyperbolic  cosecant  -CSCH(X)  *C0TII(X) 

7I!ARCSINH  ARCSINIl(X)    hyperbolic   arcsine  1/(X**2  +  1)  **.  5 

7IIARCC0SH  ARCC0SH(X)   hyperbolic   arccosine  1/ (X**2-l)**.5 

yilARCTANH  ARCTANII(X)    hyperbolic   arctangent  1/(1-X**2) 

71IARCC0TH  ARCC0TI!(X)   hyperbolic   arccotan,?ent  -1/CX**2-1) 

71!ARCSECn  ARCSFXH(X)   hyperbolic   arcsecant  -X**-l    *    (1-X**2)  **- .5 

7IIARCCSai  ARCCSni(H)    hyperbolic   arccosccant  -X**-l    *    (X**2+l)**- .  5 

31IEXP  EXP(X)  exponential  EXP(X) 

3IIL0C  L0C(X)  natural    logarithm  1/X 

It  should  be  noted  that  a  large  number  of  the  trigonometric  and 
hyperbolic  functions  whose  derivatives  are  defined  in  I.SOIDR  do 
not  have   their  values  defined   in   I.SOINI. 
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Exercises 

2.11.1.   The  definition  of  .^Hl.*  established  by  LSOINI  is  not  the 

best  choice  for  all  purposes.   When  (X,Y,Z,...)  is  intended  to 

represent  a  vector  this  definition  is  suitable.   On  the  other  hand, 

for  evaluation  purposes  in  general,  the  following  is  a  better 

definition: 

CALL  LSnDEF(3lll.*,LSQMN[,(LSnMNL(l.  ,  1011$  ,  1  .*$.$.$,  1 .)  )  , 
LSQMNL(10ll$.l.  *$.$.$)) 
Run  the  following  program  with  and  without  this  definition.   Describe 
the  difference  in  the  output. 

PR(3GRAM  TEST ( INPUT, 0l)TPUT) 

DIMLNSI0N  SPACE (5000) 

CALL  INITAS (SPACE, 5000) 

LA=LIST(9) 

CALL  LSnr)ES(LSOVAL(LA)) 

CALL    LSQDES(LA) 
C    INSERT   DEFINITI0N  0F   3111.*   HERi;   T0   REPLACE   STANDARD   DEF 

CALL   LSQDEF(3HC0S,LSOMNL(LSQMNL(1. , 
*        LSQMNL  (LSQMNL  ( 1 . )  ,  LSOMNL  (-  1 .  ,3HSIN,2 .  )  )  ,  3111 .  * , .  5)  )  ,0) 
C    I.E.    C0S   =    (1-SIN**2)**.5 
1  CALL    INLIST(LA,5LINP11T) 

IF(LISTMT(LA).EQ.O)CALL   EXIT 

CALL   LSQPNT(LA,6HS0URCE) 

LR=LSQVAL(LA) 

CALL   LSQr)ES(LA) 

CALL   LS0PNT(LB,4HVAL1) 

LC=LSOVAL(LB) 
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CALL   LSQDI-S(LB) 
CALL   LSQPNT(LC,4I1VAL2) 
CALL   LSQL)i;S(LC) 
G0  T0    1 
END 
data  record: 

C0S**2   +   SIK**2$ 

C0S{X)**2    +   SIN(X)**2$ 

C0S(X)SIN(X)$ 

$   ZER0  i:XPRrSSI0N  T0  F.NI)  INPUT 

2.11.2.   Tlie  definition  of  3III1".  as  a  conditional  expression  makes 
it  possible  for  the  user  to  easily  create  recursive  definitions  of 
expression  symhols.   For  example,  the  following  code  would  define 
the  expression  symbol  9I1I'ACT0RIAL  applied  to  the  argument  N....1  to 

he  the  factorial  of  N 1,  i.e.  N 1*(N 1  -  1)*(N 1  -  2)*.. 

.  *2*  1 

LAM=  LSQMN  L  ( LSQMN  L  ( 1  OH  ,  1 01 IV .  ( I F .  ( 

*  ,101[N 1,1,     .lOlIN 1-1,1, ion, FACT0RTAL. 

*  101I(N 1-1).10H*N 1))      )) 

CALL  LSQl)FF(i>HFA(rr0RIAL,INMST(LA,LAM)  ,LSQMNL(6IIN 1)) 

CALL  LSQDrS(LAM) 

Write  a  program  using    IF.   which  will   evaluate  binomial   coefficients 
by   the   rule: 

N^K   =   N-l4   *   N-l4-l         for  0   <    K  <    N 

c    =    r    =  1 

NO        N   N 
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2.12.   F.xamplcs 

We  now  present  some  worked  examples  of  SYMBOLANG  programming. 

2.12.1.  In  some  applications  it  may  he  a  nuisance  to  have  to 
form  expressions  representing  constants.   .So  we  write  a  function 
LSQC0N(QUAN)  which  will  return  an  expression  equivalent  to  the 
constant  OUAN.   There  is  actually  very  little  to  do.   If  OHAN  is 
zero  an  empty  list  is  the  appropriate  expression.   Otherwise  ((OtlANI) 
is  the  correct  form. 

FUNCTION  LSOr0N(OllAN) 

LSnC0N=LIST(9) 

IF(QUAN.i;0.0)Rl-TI!RN 

CAM  MANY  (LSOr0N,  I.SOMNI,  (OUAN)  ) 

RF.TURN 

FNn 

2.12.2.  Parentheses  arc  a  proat  hlessinp  when  one  wants  to  control 
the  organization  of  an  expression,  hut  often  hinder  such  tasks  as 
comparison  and  evaluation.   The  following  function  returns  a  version 
of  the  expression  US  in  wtiich  parentheses  have  to  a  large  extent 
been  removed.   Indeed  parentheses  will  remain  only  for  expressions 

to  positive  or  negative  fractional  powers  and  to  the  power  -1.  unless 
a  non-constant  power  is  encountered. 

FlINrTI0N  I,SOXPN(I.I.S) 

C0MM0N  AVSL,X(100) 

ASSinN  ion  T0  i,0r 
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LYS=LIS 

CALL  sTRniR(visiT(L0r)  .r.snxPN) 

RFTIIRN 

100  SA=SrORDR(LYS) 
LTU=LIST(9) 

101  SAT=SEQLR(SA,a') 
IF(N.nT.O)CALL  TFRMCLTII) 
SATR=SEORDR(SAT) 
LVU=LSQMNL(LSQMNL  (SF.OLR  (SATR  ,N)  ) ) 

102  LWIJ=LS0MNL(1.) 
NAR=0 

106      CALL   STRDIR(SEOLR(SATR,N),LYS) 

IF(N)103,104,105 
105      CALL   LSODES(LWU) 

LSU=LSOADD(LVU,LTU) 

CALL  LSODFS(LVU) 

CALL  LSOnnS(LTU) 

LTU=LSU 

G?)  T0    101 
104   CALL  MANY(X(100) ,NAR+l,SA,SATR,LTII,LVn,LWU) 

SAG=VISIT(L0C) 

CALL  LSnUNM(XC100),NAR,SA,SATR,LTn,LVU,LWU) 

CALL  MANY (LWU, SAC) 

C?)  T0  106 

103  IF(LSnCNMfLYS,.'5Hl.*).FO.O)f;0  T0  107 
110  CALL  MANY(LWn,LYS) 

CALL  STRniR(SEQLR(SATR,N) ,LYS) 
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IF(N)10,S,109,in5 
109   CALL  MANY(X(inO),SA,SATR,LTU,LVU,LWll) 

CALL  STI'niR(VISITfL0C),LYS) 

CALL  LSQUNM  (X  ( 100)  ,  SA ,  SATR ,  LTU ,  F.VII ,  LWU) 

IF(LSOTYP(LYS,LAS).NF..O)n0  T0  lOR 

CALL    LSQnFS(LYS) 

LYS=LAS 
108      CALL  MANY(LWU,LYS) 

LWlI=LSQMNL(LWn) 
125      LMU=LSOMF.X(LWU,LVll) 

CALL   LSODES(LW()) 

CALL   LSOnES(LVU) 

LVU=LMll 

n0  T0  102 
107   IF(NAR.NT..1)G0  T0  110 

CALL  STRniR(SFnLR(SATR,N) ,LYS) 

IFfN)llI,112,105 
111   LYS=LSOMNL(LSOMNL(LYS)) 

C0  T0  141 
121   CALL  MANY (X(l 00) ,SA, SATR, LTU, LVU.LWU) 

CALL  STRnTR(VISIT(L0C),LYS) 

CALL  LSOUNM  (X  ( 1 00)  ,  SA , SATR ,  LTU ,  LVU  ,  LWII) 
141   CALL  LSOUNM (LWl I, LCU) 

CALL  LSODRSfLWU) 

LWU=LSORAZ(La!,LYS) 

CALL  LSOI)FS(LnU) 
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CALL   LSQnES(LYS) 

^Wll=LsnTYP(LW^J,x^'.'U) 

IF(C^^wu,GE, 50000). («R.(NMJ. no. o))n0  t0  125 

IFfLSOrNM(XWIJ,31ll.*).Nr:.O)n0  T0  12.S 
rF(M0D(MWlI, 10000). KE.l)r,0  T0  125 
SAW=SCORDR(H\'lJ) 
SAW=SEORnR (SFQLR (SAW ,N) ) 
YS=SEOLL(SAW,N) 

rF(N.Gr,.o)n0  t0  125 

IF(YS)116,117,118 

117  CALL  LSOt)ES(LWn) 
f;0  T0  102 

116  M=-YS 

ASSIGN  121  T0  LAG 
G0  T0  130 

118  M=YS 

ASSIGN  122  T0  LAG 

1 30  K=M 
SN=SEQLL(SAW,N) 

CALL  STRniR(SHnLL(SAW,N) ,LARG) 
LCRG=LARG 

LBPG=LSQMNL(LSOMNL(SEQLL(SAW,N))) 
IF(K.I-.O.O)G0  T0  LAG,  (121,122) 

131  L=K 
K=K/2 
L=L-2*K 
IF(L.F.O.O)G0  T0    132 
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LDRG=I,BRG 

LBRG=LSQMEX(LDRG,LCRG) 
CALL    LSQnES(LDRG) 

132  IF(K.F.Q.O)G0  T0   133 
LURG=LCRG 

LCRG=LSQMEX (LDRG , LDRG) 
IF(LDRG.NE.LARG)CALL   LSQnES(LDRG) 
G0  T0    131 

133  IF(LCRG.NE.I-ARG)CALL   LSQDES(LCRG) 
G0  T0   LAG, (121,122) 

122      X=YS-M 

IF(X.NE.O)G0  T0   124 

CALL  LSQDES(LWU) 

LWU=LBRG 

G0  T0  125 
124   LRU=LSQMNL(LSQMNL(1. ,LARG ,3H1 . *,X) ) 

CALL  LSQDES(LWU) 

LWU=LSQMEX(LRU,LBRG) 

CALL  LSQDES(LRU) 

CALL  LSQDES(LBRG) 

G0  T0  125 
121   X=YS+M 

IF(X.NE.O)G0  T0  126 

CALL  LSQDES(LWU) 

LWU=LSQMNL(LSQMNL ( 1 . , LBRG , 3H1 .*,-!.)) 

G0  T0  125 
126      LBRG=LSQMNL(LSQMNL(1 . , LBRG , 3H1 .*,-!.)) 
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G0  T0  124 

END 
This  code  is,  perhaps  unfortunately,  typical  of  the  sort  of  code 
needed  to  derive  one  expression  from  another.   Basically,  one  sets 
up  a  recursive  loop,  first  through  the  terms  of  an  expression 
accumulating  a  sum  and  within  that  loop  another  loop  through  the 
factors  of  each  terms  accumulating  a  product.   In  this  case,  the 
sum  was  accumulated  in  the  expression  LTU,  and  the  product  in  the 
expression  I.VU,  whi  le  the  sequence  reader  SA  was  used  to  scan  for 
each  term  of  the  expression  LYS  being  expanded,  and  the  sequence 
reader  SATR  was  used  to  find  the  factors  of  each  term.  When  a 
function  argument  or  non-constant  power  is  encountered  the  readers 
and  expressions  being  used  to  accumulate  results  are  saved  on  the 
public  list  X(IOO)  and  a  \'ISIT  is  paid  to  the  same  loop  to  expand 
parentheses  there.   The  place  where  the  particular  task  of  this 
loop  is  introduced  is  in  the  test  at  statement  103.   If  we  changed 
that  statement  to  a  CONTINUE  we  would  have  a  proper  routine  for 
producing  a  copy  of  an  expression  which  is  guaranteed  to  have 
correct  ordering. 

2.12.3.   A  task  the  user  is  more  likely  to  wish  to  code  himself 
is  the  evaluation  of  the  determinant  of  a  symbolic  matrix.   The 
following  program  does  this  task  by  a  recursive  expansion  in  terms 
of  minors.   Rather  than  use  scratch  arrays,  a  list  of  struck  columns 
is  used.   As  written,  the  program  is  limited  to  5  by  5  or  smaller 
matrices,  but  the  subroutine  LSQDF.T  which  does  the  work  is  not. 
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program: 


PHOURaM    TtST(  IMPI'T, OUTPUT)  TEST0n!!2 

niMEMSIOM    LSPaY(5uOiJ  >  ,1  k4Y(!?,!3)  TfcSTU(J(^ 

CALL     1MITaS(LSP:aY,'^U^U  )  TtST0(in4 

RfcAU    «:0,M;Ii^.  Ttr'TCiir'tJ 

FORMATdl)  T6?;T00f'o 

CALL    laETlNCLf  AY,f'l;T,l)  TcSTPUl^? 

EKD  TfcSTOord 

rjNCTlO.M    LSuCtT(L^^iY  ,  M:''U')  LbRDn002 

COMMON    AVbL,x(i()n)  LiODPUfi 

DlME^)bIU^i    LKA  Mi^jDI'^/Ml 'h)  ^5000004 

ASSIGN    iPu    TO    LUC  LSGOOgrs 

AS^Ji^niNj    200     Tn     IFlun  LSdO'JOPb 

LMT=LbQMML(LSjfNH-l.}>  LSnuOOC? 

LSTRUCsLlSTC^)  LSOiJOor.b 

IURG=1  Lsounon9 

JQRG=1  LSPDOOIO 

NSIZ=NDIM  LS^OOUll 

CALL    i)TRPI«(  VlSITaOi")  ,LSQ3hT)  LSr.iO0(jl2 

CALL    LSantSCLbTHIJC)  LSnDPu]3 

CALL    L?QDt"^(L-^.n  LSi'jDOUI'* 

RfcTuPN  LsnuPuis 

lF(MS12.bT,:')(;()     TO    ^.u'f'O  USCDCUlb 

MGRGsIOt-r;*!  La'JiJOUl? 

CALL    VISITC  ir  if' jj)  LSOIJPOI.8 

LTU  =  LSOMEX(Li^AY(IUi^ti,J^-^r;),L«AY(MUPGiMOkr,  +  i. ))  LbPunul'i 

LTV  =  LSPMPX(Lf=AYnORa,J"Wf:*1).L«AY(wCP(j,jDR'^))  LSOUD020 

LTw=Lbr3MFx(  L^I  .LTV)  Li>nuno2i 

CALL    L5QDbS(LTV)  LSOG0u?2 

LTXsLSOADLMLTU,LTW)  LSOliOU2«J 

CALL    LSURbS(LTl')  LSOO0U24 

CALL    LSQDfc9(LTW)  LSOD0025 

CALL    TERM(LTX)  USOD00?6 

JORG=JUPb*l  LSODOO?7 

NbIZ  =  MSIZ-^  LSODOO'-^tt 

^!ORG=inRGfl  LSODP029 

CALL    VISTT( IP  r-D)  LSnjPCJO 

CALL    hANY(LiTK''C,  1  Ortnj  L^Jt.'QOU?! 

lOBlisMORG  LiiO0U032 

CALL    bTHDIR(VlSlT(LUi")  ,LTU)  LSnD0o33 

MORG=inKG  LSrDOU>'^'» 

CALL    LSOUfMM(LSTHL'C.  lOKf^O  LSr:O0035 

LTWsLbnr-lEXtLRAvciuRG.j'^wn-D.LTu)  LSOuOu3  6 

CALL    LSQDbS(LTl')  LSODOQ?? 

Ki>IZ  =  ''JSI7  LbPDnO?8 

LSIGN'  =  -1  L3ODn039 

KSI7=KSIZ-1  LS^aou^o 

CALL    i^AIMY(LSTKUC,MO,-if;)  LS'^OOO'H 

CALL    MAix|Y(X(99)  ,KbIZ.L'=:i'^>M.LTw)  Laf;uOO<»2 

CALL    bTRDIR(  VI<:iT(LUC)  ,LTU)  LSOD0043 

CALL    LSUUNM(X(09)  .K3IZ,L'^I'-''»lTk)  LbODPtJ«*4 

CALL    LSQUiMMCLbTwuU.MHKn)  LSODPO'^tJ 

lF(LSinN, GT,  n  )r.o   To   p^'^o  Lsoono^b 

LTMsLbOMEX(Lf"I  ,LTU)  LSOonCZ 

CALL    LSQDcS(LTli)  LbOU0048 

LTijsLTM  LSOun049 

LSIGi^^  =  -LbIG^'  LbODOuSO 
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3ono 

200 
201 

202 

203 


100 


200 


UTMsL 

CALL 

LfUsL 

CALL 

CALL 

LTWsL 

IF(KS 

HORG  = 

CALL 

GU  TO 

NbTZ  = 

CALL 
SArS 

KFLG  = 
CALL 

I F  ( IN  , 

It-"  (NO 

IF(KF 
RO  TO 

KTLG 
GO  TO 
EN" 
SURRO 
D  I  M  fc  Ni 
DO  10 
DU  10 

LRAY 
L 
CALL 

ro  20 
no  ijo 

CALL 
CALL 
RtTuR 
END 


SOrtE 

Lson 

aOAD 

Lsar 

L?QD 
TU 
17. L 

MOHG 

?or 

=  JUR 
WSIZ 
TEH'^ 
tQRD 
0 

SIR 
bT,n 
HG-M 
sMOR 
LG.E 

20n 
=  1 

201 

OTIN 
blON 
U  J  = 
U  K  = 
(  J.K 
tx  =  L 
LSUP 
0  J  = 
0  K  = 
LSO 
LSOD 


X(L'^AV(MOWu.jr'wn-i),LTu) 

tS(LTll) 

Li  (  L  T  «  ,  L  ■>■  ii  ) 

E  S  (  L  T  W  ) 

t-^lLTK') 

t.OJijP    TO    ,i"(if 

♦1 

T(  IPIMJ) 

rj 

U-1 

*1 

(LT'O 
H(LSTi:ur) 

bIR(bPULK(b4,^'),N0.-<G) 

)CALL    TEHM(n,) 

0RG)2ri,^n2,?P3 

G  +  1 

U,J)ur    TO    2T1 


fc    UF 

LhiA 

l.iMH 

1 .  f'i  n 

)  =  IN 
SOUP 
(viT(L 

1,  -JD 
LlES( 
E?<L 


N(LRAV,  ^'l■IM) 


ST(LWAY(  J.K  i  .'JLlMf'UT) 

LRAY,ML!(^) 

,4Hntrw  > 


AY( J.K  )  5 
) 


LSDD0G51 
LSOur'0?'2 
LSOQOO'SS 

LaCUno54 
L5OD0  0?t> 
LSOL)0036 
LSQD0057 
LSOJ0056 
LiiOUOoSV 
LSOJ0060 
LSODOOfrl 
LSQUOO<-2 

LbODOfjei 

LSODOOf''* 
LSOD0065 
LSOD0066 
LSOD00^7 
LSODOOfiS 
USRU00^9 
LSO'jOo7o 
LSPD0U71 
LSQD0072 
LSOD007i 
LSCD0074 
GET  100  02 
GETlPOOi 
GETIOOI'4 
RETIOOPb 
GcTlDorfc 

n  E  T 1  n  0  0  7 

GcTinotJ8 
'.JETIOOOy 
GET  I  0010 
GETIOOll 
3ET10012 
GET10013 
GET10014 
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2.57 
data  record: 


A114 
Al?i 
A 13* 
Alts, 
A2ji 
Air?* 

Ai:3$ 
Ai:4i, 
Ail  $ 
Ai?i 
As53i 
AJ4* 
A41  % 
A4?i 
A434 
A44S, 


The  data  provided  is  a  dimension  between  2  and  5  (in  this  case  4) 
followed  by  the  elements  of  each  row  of  the  matrix. 


2.58 


output: 

DET"    =    A11«A22*A33*A44         -         AiJL*  4^;>*  a  j-i*  A43         -  411*A27 

♦  A.1^*A44    *    A1.1*A?3*A  .5  4*^42    ♦    A1. 1*A?4*A'^2*A43 
All*A2'l*A3:^«A4?    -    A1'''*A21*V33*A<»4    *    4:,?*a21*A>54 

♦  A43    +    a12*a?3*A"»1*aA4    -    A12*A23*A74*.\4-i. 
A12«A24*A31  *A43    *    4i;^<iA24*.433*A<ll.    ♦    Ai3*42l*A3? 
*A44    -    A13*A?l'-A'^4*'-*'*2    -    A13*A?2*A31*'*44    ♦ 
A13»A22*A3'>*A41     +     Al  '!*A24*A3l«A42     -     A  .1.3*  a24*  A3? 

•  A^l    -    Al4*A?i*A3  2*'*''o    ♦    A14*A?1*a33*a42    ♦ 
A14*A2?*A31  "A^S    -    A1'»«A22*A33*A41    -    Ai4*423*A3l 

•  A42    ♦    A14*A23*A-^2*A'il 
% 

SEND  OP"  EXPRESSION 
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3.   BELLS,  WHISTLES,  AND  FRILLS 

The  information  in  this  section  should  not  be  needed  for  most 
SYMBOLANG  applications.   Here  may  be  found  alternative  routines  and 
calls  to  routines  mentioned  previously  which  may  provide  some  coding 
convenience,  or  extra  efficiency,  or  access  to  internal  hooks. 
Furthermore,  some  routines  which  are  not  recommended  for  new  code 
are  included  to  assist  in  reading  existing  programs.   Niceties,  fine 
distinctions,  special  cases,  and  funny  tricks  (if  any)  are  largely 
confined  to  this  section.  The  prograiimer  wishing  to  write  code 
which  involves  detailed  structural  analysis  of  expressions,  should 
find  this  material  of  interest. 

3.1.  Additional  List  and  Expression  Creation  Methods. 

In  addition  to  the  methods  described  in  §2.1,  calls  to  the 
functions  CPYTRM,  INLSTL,  INLSTR,  LIST,  LSQCPY,  LSQCXP,  NEWB0T, 
NEWT0P,  NULSTL,  NULSTR,  NXTLFT,  NXTRGT,  PUTLST,  SUBSET,  SUBST,  and 
SUBSTP  provide  means  of  constructing  lists  and  expressions  from 
scratch  or  from  existing  lists.   Some  of  these  calls  could  also  be 
considered  as  methods  of  destroying  lists,  and  will  be  mentioned 
again  in  13.4. 

CPYTRM(ZtSi,  lis2)    removes  the  top  element  from  list  or 
expression  lis\    and  places  it  on  the  bottom  of  list  or  expression 
lis2-      If  lisx    is  empty,  so  there  is  no  top  element  to  remove,  then 
no  change  is  made  in  either  argument,  and  CPYTRM  returns  0.,  else 
CPYTRM  returns  1.  as  its  value.   If  both  arguments  are  expressions. 
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the  effect  is  to  subtract  the  leftmost  term  from  lisi   and  add  it 
to  li82-      Such  use  is  not  recommended  since  correct  ordering  is 
not  necessarily  preserved,  and  the  destructive  change  to  lis\   may 
affect  expressions  of  which  it  is  a  subexpression.   This  routine 
should  not  be  used  in  new  code. 

INLSTL(H8,  quan)   removes  all  the  elements  from  list  or 
expression  lis   and  inserts  them  to  the  left  of  the  list  element 
specified  by  the  quantity  quan   as  a  sequence  reader  (see  §2.10). 
The  list  element  before  which  the  insertion  is  to  be  made  must 
reside  on  a  list  distinct  from  lis,   and  should  be  the  next  element 
to  which  the  sequence  reader  quan   would  make  a  right  advance.   INLSTL 
returns  the  emptied  list  lis   as  its  value.   No  attempt  is  made  to 
change  the  quantity  quan.      It  is  not  advanced  as  a  reader,  and  thus 
will  still  point  to  the  element  to  the  right  of  the  insertion  on 
the  right  and  to  the  element  to  the  left  of  the  insertion  on  the  left. 
As  in  §2.10,  the  head   of  a  list  is  considered  an  element  of  the 
list  in  this  context.   For  example,  a  list  (1,  2,  3,  4,  5,  6)  may 
be  created  by  the  folowing  code. 

LA=LSQMNL(1,4) 

LB=LSQMNL(2,3) 

LC=LSQMNL(5,6) 

SA=SEQRDR(LA) 

SAC=SEQLR(SA,N) 

CALL  LSQDES(INLSTL(LB,SA)) 
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SAC=SEQLR(SA,N) 

CALL  LSQDES(INLSTL(LC,SA)) 
This  makes  LA  the  lists  (1,4),  (1,2,3,4),  and  (1,2,3,4,5,6)  in  turn, 
while  SA  points  to  1  at  creation,  to  4  on  the  first  advance,  and 
to  the  list  head  on  the  second  advance. 

In  view  of  the  destructive  manipulations  performed  by  INLSTL, 
its  use  on  expressions  should  be  avoided. 

INLSTR(Zts,  quart)    removes  all  the  elements  from  the  list 
or  expression  lis   and  inserts  them  to  the  right  of  the  list  element 
specified  by  the  quantity  quan   as  a  sequence  reader.  The  list 
element  after  which  the  insertion  is  to  be  made  must  reside  on  a 
list  distinct  from  lis,    and  should  be  the  next  element  to  which 
the  sequence  reader  quan   would  make  a  right  advance.   INLSTR  returns 
the  emptied  list  lis   as  its  value.  With  obvious  modifications,  the 
remarks  for  INLSTL  apply  to  this  routine.  The  differences  may  be 
seen  by  considering  the  following  code,  which  creates  a  list 
(1,  4,  2,  5,  6,  3). 

LA=LSQMNL(1,4) 

LB=LSQMNL(2,3) 

LC=LSQMNL(5,6) 

SA=:SEQRDR(LA) 

SAC=SEQLR(SA,N) 

CALL  LSQDES( INLSTR (LB, SA)) 

SAC=SEQLR(SA,N) 

CALL  LSQDES( INLSTR (LC.SA)) 
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This  makes  LA  the  lists  (1,4),  (1,4,2,3),  and  (1,4,2,5,6,3)  in  turn, 
while  SA  points  to  1  at  creation,  to  4  on  the  first  advance,  and 
to  2  on  the  second  advance. 

As  with  INLSTL,  the  use  of  INLSTR  on  expressions  should  be  avoided. 

LIST(O)  returns  a  newly  created  empty  list,  () ,  which  is 
the  correct  representation  of  [[0]].  This  result  agrees  with  that 
of  LIST(9),  except  that  the  list  thus  created  acts  as  if  it  were 
already  a  sublist  of  some  other  list,  and  will  not  be  destroyed 
by  either  a  call  to  LSQDES  nor  by  the  erasure  of  lists  of  which 
it  is  subsequently  made  a  sublist.  The  function  IRAI-ST  described 
in  §3.4  must  be  used  to  erase  a  list  created  in  this  manner.   Such 
a  list  is  fairly  well  protected  against  accidental  erasure.   For 
example,  it  may  still  be  referenced  after  use  as  an  argument  in  a 
call  to  LSQDEF  (see  §2.7)  or  LOOKUP  (see  §2.9),  though  subsequent 
modification  of  the  list  by  the  user  could  have  disastrous  consequences. 
If  an  expression  has  been  created  by  other  means,  and  has  not  been 
manipulated  beyond  the  point  of  creation,  this  call  to  LIST  might 
be  combined  with  a  use  of  INLSTL  to  create  a  protected  version  of 
the  expression.   For  example,  if  LA  is  an  expression  just  formed 
by  INLIST,  the  following  code  will  have  that  effect. 

LB=LIST(0) 

CALL  LSQDES (INLSTL(LA,SEQRDR (LB))) 

This  leaves  LB  as  a  protected  version  of  LA,  but  destroys  LA.   On 
the  other  hand,  a  much  safer  way  of  protecting  expressions  is 
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available,  even  if  one  does  not  wish  to  explicitly  make  them 
sublists  of  a  known,  well  protected  list.   If  LA  is  a  list  or 
expression 

CALL  SETIND(-1,-1,LCNTR(LA)*1,LNKR(LA)+1) 

will  make  the  list  or  expression  look  as  if  it  were  a  sublist  of  one 
more  list  than  it  actually  is.   Subject  to  the  previous  comments 
about  using  such  protection  at  all,  this  last  method  may  be  used 
on  any  list  or  expression  at  any  time  after  its  creation  and  prior 
to  its  destruction.   If  this  method  is  used,  IRALST  is  needed  to 
erase  it  rather  than  LSQDES. 

LSQCPY(yar)  returns  a  newly  created  list  which  contains 
essentially  the  same  elements  as  are  pointed  to  by  the  variable 
or  array  element  var   treated  as  a  sequence  reader  of  some  list 
or  expression.   The  elements  copied  onto  the  newly  created  list 
are  those  which  are  reached  by  right  advances  of  the  sequence 
reader  quan   until  it  reaches  the  head  of  the  list  it  is  scanning. 
The  sequence  reader  quan   is  itself  advanced  in  the  process. 
For  example,  the  following  code  will  leave  LA  as  the  list  (1,  2,  3) 
and  LB  as  the  list  (2,  3) . 

LEX=LSQMNL(1,2,3) 
SA=SEQRDR(LEX) 
SAC=SEQLR(SA,N) 
LB=LSQCPY(SA) 
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U=LSQCPY(SA) 

LA  is  a  full  copy  of  the  list  LEX  since  SA  was  advanced  all  the 
way  to  the  head  of  LEX  in  making  the  partial  copy  for  LB.  This 
routine  may  be  used  to  properly  break  off  terms  from  expressions. 
For  example,  in  the  following  code  the  first  term  of  expression 
LEX  is  placed  on  expression  LA,  and  the  remaining  terms  are  placed 
on  expression  LB.   LEX  is  not  changed. 

SA=SEQRDR(LEX) 

SAT=SEQLR(SA,N) 

LA=LIST(9) 

IF(N.EQ.n)CALL  MANY(LA, LSSCPY (SAT) ) 

LB=LSQCPY(SA) 

For  elements  which  are  sublists,  LSQCPY  will  use  the  routine 
LSQCXP  (see  below)  to  either  provide  the  sulUist  itself  or  a  copy 
via  LS-SCPY.   Thus  LSQCPY  may  not  provide  a  completely  fresh  copy 
of  a  list  or  expression,  and  sublists  of  the  copy  it  returns 
should  be  treated  with  the  same  consideration  afforded  the  sublists 
of  the  original.   Normally  LSQCXP  does  return  the  sublist  itself, 
which  implies  that  LSQCPY  will  return  a  copy  made  only  on  the 
first  level  having  sublists  in  common  with  the  original  list. 

LSQCXP(Zie)  either  returns  the  list  or  expression  lis,   or 
returns  a  full  copy  of  lis   via  LSSCPY.   The  choice  is  determined 
by  the  state  of  the  first  variable  (usually  called  LSQCSX)  in  the 
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C0MM0N  block  labeled  LSQCSX .   If  this  variable  is  zero  then  the  list 
itself  is  returned,  else  LSSCPY(Hs).   This  variable  is  initialized 
by  a  DATA  statement  to  zero.   Thus,  unless  the  user  resets  this 
variable  to  a  non-zero  value,  LSQC.'XP  will  return  the  list  or 
expression  itself.   LSQCXP  is  used  by  the  other  routines  of  the 
package  which  generate  new  expressions  from  old,  such  as  LSQAUD 
and  LSQMF.X,  to  determine  when  common  terms  and  subexpressions  are 
permitted.  Thus,  if  the  user  wishes  to  prevent  the  introduction 
of  common  subexpressions  he  should  use  code  of  the  form 


C0MM0N/ LSQCSX/ LSQCSX 


LSQCSX=1 


This  routine  is  largely  for  internal  use  of  the  package. 

NEWB0T(<7ua?7,  lis)   places  the  quantity  quan   on  the  bottom 
of  the  list  or  expression  lis.     The  value  returned  may  be  used 
as  a  sequence  reader  whose  next  right  advance  would  bring  it  to 
the  element  quan   just  added  to  the  bottom  of  lis.      Until  such  a 
right  advance  is  made,  no  left  advance  should  be  made  since  the 
reader  begins  with  only  a  pointer  to  the  right  (to  quan)   but  with 
no  valid  pointer  to  the  left.  Aside  from  the  value  returned 
MANY (its,  quan)   has  the  same  effect. 
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NEWT0P (<7uan ,  lis)   places  the  quantity  quan   on  the  top 
of  the  list  or  expression  lis.     The  value  returned  may  be  used 
as  a  sequence  reader  whose  next  right  advance  would  bring  it  to 
the  element  quan   just  added  to  the  top  of  lis.     The  same  strictures 
apply  to  this  sequence  reader  as  to  the  value  of  NEWB0T. 

NULSTLC^uan,  lis)   returns  a  newly  created  list  formed 
by  stealing  the  elements  to  the  right  of  the  head  of  list  or 
expression  lis   and  to  the  left  of  the  next  element  to  be  reached 
by  a  right  advance  of  the  quantity  quan   treated  as  a  sequence 
reader,  as  well  as  stealing  that  element  itself,  provided  there 
is  such  an  element.  The  sequence  reader  quan   must  refer  to 
elements  of  the  list  lis.      If  the  element  quan   would  reach  by  a 
right  advance  would  be  the  head  of  the  list  lis,   then  no  change 
is  made  in  lis   and  an  empty  list  is  returned.  Otherwise  this 
constitutes  a  destructive  change  in  lis   and  should  normally  be 
avoided  for  expressions.   As  an  example  of  the  effect  of  NULSTL 
consider  the  following  code  which  leaves  LA  as  the  list  (3,  4) 
and  LB  as  the  list  (1,  2) . 

LA=LSQMNL(1,  2,  3,  4) 
SA=StQRDR(LA) 
.SAC=SEQLR(SA,N) 
LB=NULSTL(SA,LA) 
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NULSTR((7i/an,  lis]    returns  a  newly  created  list  formed 
by  stealing  the  elements  to  the  left  of  the  head  of  the  list 
or  expression  lis   and  to  the  right  of  the  next  element  to  be 
reached  by  a  right  advance  of  the  quantity  quan   treated  as  a 
sequence  reader,  as  well  as  stealing  that  element  itself,  provided 
there  is  such  an  element.  The  sequence  reader  quan   must  refer  to 
elements  of  the  list  lis.      If  the  element  qvan   would  reach  by  a 
right  advance  would  be  the  head  of  the  list  lis,   then  no  change 
is  made  in  lis   and  an  empty  list  is  returned.  Otherwise  this  is 
another  destructive  change  in  lis   and  should  normally  be  avoided 
for  expressions.   As  an  example  of  the  effect  of  NDLSTR  consider 
the  following  code  which  leaves  LA  as  the  list  (1)  and  LB  as  the 
list  (2,  3,  4). 

LA=LSQMNL(1,  2,  3,  4) 
SA=SEQRUR(LA) 
SAC=SEQLR(SA,N) 
LB=NULSTR(SA,LA) 

NXTLFT (cjuon  1  ,  quan2)    inserts  the  quantity  quan\   as  a 
new  element  in  a  list  to  the  left  of  the  element  which  would 
be  reached  by  a  right  advance  of  the  quantity  quan^   treated  as 
a  sequence  reader. 
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aXTRGT iquan I ,  quan2)    inserts   the  quantity  quani   as  a 
new  element   in  a   list   to  the  right   of  the  element  which  would  be 
reached  by  a  right  advance  of  the  quantity  quan2   treated  as   a 
sequence  reader. 

PUTLST(Hs,   quan^  \  ,  quan^2'    •••>  quan^ ,    ,    lis,  quan^  \  , 

'^1 


quari^-,    ,    lis,    ...,    lis,  quan   \,    ...,  quan  ,     ,    lis,    ....    lis,  quart   i, 

...,  quan  ,    ,    lis,    lis)    for  r=   1,    2,    ...,   fej,   k2,    ...,   k     =   1,    2,    ...; 

r 
and  up  to  27  arguments  in  all;  returns  the  expression  lis,   onto  which 

the  terms  (quan   j,  ....  quan   ,  )  have  been  added.   If  these  are 

i 
valid  terms  the  result  is  a  valid  expression.  Thus  each  quan   i  must 

be  a  non-zero  floating  point  constant,  and  so  forth.   For  example, 

a  valid  representation  of  [[x  -  y^]]   may  be  left  in  LA  by  the  following 

code. 

LA=LIST(9) 

CALL  PUTLST(LA,1.  .IIIX.I.  ,LA,-1.  ,  IHY,  3.  ,  LA,  LA) 

The  use  of  PUTLST  should  be  avoided  in  new  code. 

PUTLST(0,  lis,  quan^ I ,    ...  )  with  arguments  as  above 
except  for  the  leading  0,  acts  as  the  above  call  to  PUTLST  except 
that  re-ordering  of  factors  within  the  terms  is  prevented.  This 
call  should  also  be  avoided  in  new  code. 

SUBSBT(i:7t«2n,  lis)    substitutes  the  quantity  quan    for  the 
bottom  (rightmost)  element  of  the  list  lis.      It  is  an  error  to 
use  this  routine  if  the  list  lis   is  empty.  The  value  returned 
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is  the  element  which  quan   replaced.   Since  this  is  a  destructive 
manipulation  of  lis,    lis   should  not  be  an  expression.  The  following 
code  will  leave  LA  as  the  list  (1,  2,  4). 

LA=LSQMNL(1,  2,  3) 
CALL  SUBSET (4, LA) 

SUBSTiqucmi,  quan2)    substitutes  the  quantity  quani    for 
the  list  element  which  would  be  reached  by  a  ri^ht  advance  of  the 
quantity  quan2   treated  as  a  sequence  reader.   It  is  an  error  to 
use  this  routine  if  the  element  to  be  replaced  is  the  head  of  a 
list.  The  value  returned  is  the  element  which  qucmi    replaced.  The 
following  code  will  leave  LA  as  the  list  (1,  4,  3). 

LA=LSQMNL(1,  2,    3) 
SA=SEQRDR(LA) 
SAC=SF,QLR(SA,N) 
CALL  SUBST(4,SA) 

SUBSTP((7uan,  Hi;)    substitutes  the  quantity  quan   for  the 
top  (leftmost)  element  of  the  list  lis.     The  same  strictures  apply 
to  SUBSTP  as  to  SUBSET,  and  the  value  returned  is  again  the  element 
which  quan   replaced.  The  following  code  will  leave  LA  as  the  list 
(4.  2,    3). 

LA=LSQMNL(],  2,  3) 
CALL  SUBSTP (4. LA) 
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Exercises 

3.1.1.   Devine  the  purpose  of  the  following  function.   (Hint  -  LIS 

is  an  expression) . 

FIINCTI0N   LSQRHM(LIS) 

SA=SEORDR(MS) 

SAT=SEQLR(SA,N) 

LSQREM=LSQCPY(SA) 

RETURN 

END 
.^.1.2.   Run  the  following  program.   Discuss  its  sins. 

PR0GRAM  SNAFU ( INPUT, 0UTPirr) 

DIMENSI0N  SPACE (5000) 

GAEL  INITAS (SPACE, 5000) 

LA=INLIST(LA,5LINPUT) 

LB=INI.IST(LB,5LINPl)T) 

LCxLSOALn)(l,A,LB) 

ED=LSQMEX(LA,LB) 
1    CALL  LS0PNT(LA,2H1,A) 

CALL    LSQPNT(LB,2IILB) 

CALL   LSQPNT(LC,2HLC) 

CALL   LSQPNT(LD,21ILI)) 

SA=SEQRnR(LA) 

SAT=SEnLR(SA,N) 

SA=SEQRDR(SAT) 

CALL  NXTRC,T(LB,SA) 

CALL  CPYTPJ^(LA,LB) 

G0  T0    1 
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END 

data  record: 

siN[X)*r0S(X)-siNfX)*r0s(X)**-i$ 

SrN-(X)-C0S(X)+SIN(X)**-l*C0S(X)$ 
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3.2.  Additional  Expression  Input  Capabilities 

The  function  INLIST  described  in  §2.2  is  somewhat  more  liberal 
in  accepting  FORTRAN-like  expressions  for  translation  than  we  have 
indicated  thus  far.  As  was  mentioned  at  the  end  of  §2.2,  the  use 
of  the  "equals"  sign,  =,  is  permitted.  This  is  achieved  by  using 
a  slightly  different  syntax  than  was  shown.  The  definition  of  an 
expression  in  the  input  stream  is  actually 

<expression>  *■   <expression  part>  {  =  <expression  part>}o 
<expression  part>  ■*■   <term>  {<addition  operator><term>}o 

rather  than 

<expression>  *■   <tenn>  {<addition  operator><term>}o 

which  effectively  adds  =  as  a  fifth  class  of  operator  to  be  applied 
after  the  infix  addition  operators  ♦  and  -.  This  syntax  is  applied 
at  all  levels  of  the  scan.  Thus  the  equals  sign  may  be  properly 
embedded  deep  within  an  expression.   For  example,  the  following  is 
valid  input. 

LA=1+SIN(ARG(A=B,C=D)=0)$ 

If  the  call  to  INLIST,  INLIST(uar,  quan)  is  used,  the  expression 
parts  to  the  left  of  an  equals  sign  are  discarded  in  translation 
to  internal  representation.  In  this  case,  the  input  above  would 
be  equivalent  to 


5-15 


1-^SIN(0)$ 

which  makes  the  discarded  expression  parts  little  more  than  comments, 
There  is  another  call  to  INLIST  however,  which  allows  the  user  to 
make  use  of  the  expression  parts  to  the  left  of  equals  signs. 

INLlST(uar,  quan,    3MVAL)  returns  a  list  representing  the 
next  group  of  expressions  found  on  the  logical  unit  or  list  quan. 
The  variable  or  array  element  var   will  also  contain  the  newly 
created  expression.  The  scan  of  the  input  stream  and  the  syntax 
applied  are  the  same  as  for  the  two  argument  call  to  INLIST. 
However,  groups  of  expressions  are  translated,  rather  than  a  single 
expression.  The  result  returned  is  the  value  of  the  first 
expression  in  the  group,  with  evaluation  performed  on  this  first 
expression  after  each  of  the  following  expressions,  if  any,  has 
been  evaluated  and  erased.  The  expressions  considered  are  all 
those  encountered  in  a  scan  of  the  input  stream  up  to  a  double 
termination,  i.e.  two  dollar  signs  appearing  consecutively  within 
columns  7  through  72  of  lines  of  the  input  stream  (with  or  without 
intervening  blanks),  or  two  consecutive  lines  which  have  columns 
1  through  5  non-blank  (and  columns  7  through  72  blank) ,  or  some 
mixture  of  terminations.  The  first  expression  in  a  group  begins 
a  new  line,  but  the  remaining  expressions  in  a  group  begin  directly 
after  the  terminators  of  the  prior  expressions.   For  this  purpose, 
termination  of  an  expression  by  non-blank  columns  1-5  counts  as 
a  dollar  sign  in  column  7  of  the  terminating  line,  and  a  new 
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expression  in  the  same  group  may  begin  with  column  8  of  the 
terminating  line. 

The  translation  of  the  equals  sign  differs  from  the  two 
argument  call  to  INLIST.   In  this  three  argument  call  to  INLIST, 
expression  parts  to  the  left  of  equals  sign  are  not  ignored. 
Instead,  input  of  the  form 

...  ei  =  62  =  ..  .  =  ej,  ... 

is  translated  into 

.  .  .   EQUAL,  (ei  ,  02.  ••-.«;>)   •  •  • 

using  the  expression  symbol  EQUAL,  (i.e.  6HEQLIAL.),  which  is  not 
normally  defined  by  SYMBOLANG  initialization  (see  §2.11).  Thus 
the  input 

LA=l+SIN(ARn(A=B,C=D)=0)$$ 
is  translated  into  the  equivalent  of  the  input 

EQUAL. (LA,1+SIN(EQUAL. (ARn(EQUAL. (A.B), 
*    EQUAL. (C,D)),0))$$ 

Note  that  two  dollar  signs  are  needed  to  terminate  the  scan. 

By  defining  the  expression  symbol  6HEQUAL.  via  LSQIJEF,  the 
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user  may  cause  this  three  argument  call  to  INLIST  to  return  any 
reasonable  translation  for  the  equals  sipn.  One  approach  is  to 
make  EQUAL,  a  means  of  introducing  definitions,  defining  a  name 
on  the  left  hand  side  to  be  the  expression  on  the  right  hand  side. 
Then  the  appearances  of  those  names  defined  in  the  second  and 
later  expressions  of  a  group  would  be  replaced  by  their  definitions 
in  the  first  expression  of  the  group.  This  technique  is  detailed 
in  §3.6. 

There  is  a  function,  LISTEN,  which  consists  of  nothing  more 
than  a  two  argument  call  to  INLIST  for  input  from  the  standard 
input  unit  (5LINPUT),  i.e. 

LIST0N=INLIST(ZT:s,5LrNPirr) 

Thus,  in  order  to  read  expressions  from  a  data  record  in  his  input 
stream,  the  user  may  have  code  such  as 

LA=LIST?IN(LA) 

which  will  read  an  expression  from  current  position  to  a  single 
terminator. 
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3.3.   Additional  Expression  Output  Capabilities 

The  functions  LSQPNT  and  LSQ0UT  described  in  §2.3  have  additional 
calls  available  to  facilitate  meeting  special  output  needs. 

LSQPNT(0,  SIlFLpAT)  causes  all  subsequent  output  of  non- 
zero numbers  by  LSQPNT  to  include  a  decimal  point,  even  if  indistinguishable 
from  integers.   Unless  this  call  is  made,  numbers  within  expressions 
which  have  no  fraction  are  printed  without  a  decimal  point.   For  example 
the  list  ((1.), (2.5,1HX,1.))  which  represents  [[1  +  2.5x]]  would  usually 
be  printed  as 

hoi   =    I      *     2.5*X 
$ 
$nNI)  0F  EXPRESSI0N 

but  after  this  call  would  be  printed  as 

hot   =    I.    *   2.5*X 
$ 
.$END  0F  EXPRESS I0N 

which  is  useful  in  applications  where  SYMBOLANG  is  being  used  to 
prepare  FORTRAN  statements. 

LSQPNT(0,  3HFIX)  reverses  the  effect  of  a  call  to 
LSOPNT(0,  5nFL0AT)  and  allows  decimal  points  to  be  dropped  where 
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the  number  being  output  is  an  integer,  i.e.  has  no  fraction. 
Whether  or  not  this  call  is  made,  the  number  0  always  prints 
without  a  decimal  point. 

LSQPNT(0,  5HN?)PRE)  causes  all  subsequent  output  of 
expressions  to  lack  the  leading 

hoi  = 
and  the  trailing  "$"  and  "$END  0F  F.XPRESSI0N".   Output  will  not 
necessarily  begin  on  a  new  line.   Rather  it  will  follow  on  the 
same  line  as  any  other  output  via  LSQ0UT  which  has  not  been  forced 
out  of  the  internal  buffer  in  LSO0UT  either  by  filling  the  buffer 
or  by  a  call  to  LSO0UT(5HFLUSH) .  The  last  line  of  the  expression 
will  not  be  forced  out  of  the  buffer  in  LSO0nT.  This  form  of 
output  is  of  very  little  value  unless  a  call  to  LSQPNT(0,  5IIN0SUB) 
has  been  made. 

LSQPNT(0,  3IIPRE)  reverses  the  effect  of  a  call  to 
LSQPNT(0,  5HN0PRE),  returning  LSOPNT  to  its  initial  state  in  this 
respect. 

LSQ0UT(51IFLUSH)  flushes  the  internal  buffer  in  LSQ0UT.   If 
the  buffer  happens  to  be  empty,  this  call  has  no  effect.  Thus, 
when  LSQPNT  has  been  used  after  a  call  to  LSOPNT(0,  5HN0PRE)  it 
should  be  followed  by  this  call  to  LSQ0UT. 
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LSQ0UT(O,  hoi)   outputs  the  characters  to  the  left  of  the 
first  blank  in  the  Hollerith  quantity  hoi.      Up  to  ten  characters 
may  be  output.   If  the  leftmost  character  is  a  blank,  this  call  to 
LSQ0UT  has  no  effect. 

LSO0UT(iwt,  holarray)  outputs  the  leftmost  int  characters 
of  the  Hollerith  string  holarrar^ ,  inclusive  of  blanks.  The  integer 
quantity  int   must  be  greater  than  zero. 

LSQ0l)T(6IIMARr,IN,  int)    sets  the  left  margin  used  by  LSQ0I1T 
to  the  integer  quantity  int,   which  may  range  from  1  through  the 
current  value  of  the  right  margin.   The  left  margin  starts  set  to 
column  7,  but  is  adjusted  frequently  by  LSQPNT  to  indent  lines  of 
expression  output. 

LSQ0lJT(4Hr:[)(;F.,  int)    sets  the  right  margin  used  by  L.SO01IT 
to  the  integer  quantity  int,   which  may  range  from  the  current  value 
of  the  left  margin  through  150.  The  right  margin  starts  set  to 
column  72,  but  is  adjusted  frequently  by  LSQPNT.   I.SO^IJT  formats 
lines  in  its  buffer  starting  at  the  left  margin,  and  may  use 
columns  up  to  and  including  the  right  margin,  but  in  order  to  avoid 
splitting  Hollerith  quantities  unnecessarily,  rarely  actually  reaches 
the  right  margin. 

LSQ0UT(5HCFLAn,  int,   hoi)    causes  column  int   to  be  used 
to  mark  continuation  lines  with  the  leftmost  non-blank  character 


104 


3.21 


of  the  Hollerith  quantity  hoi.      The  integer  quantity  int   may  range 
from  1  through  150,  and  need  not  lie  within  the  margins.  Whenever 
LSQ0UT  starts  a  new  line  because  the  current  line  is  full  (rather 
than  because  of  a  call  to  LSO0UT(5ilFLUSH))  the  continuation 
character  is  placed  in  column  int   of  the  new  line.   LSQPNT  makes 
no  attempt  to  override  such  a  call,  so  continuation  marks  may 
be  used  in  outputting  expressions. 

LSQ0UT(5HrFI.An,  0)  cancels  the  selection  of  a  continuation 
mark.  This  is  the  normal  state  of  LSQ0UT. 

LSQ0UT (<?Ka«  1  ,  quarii,    •••,  quan,),   k   =  1,  2 60, 

has  the  same  effect  as  sequential  calls  to  LSO01JT.  The  argument 
list  is  scanned  from  left  to  right,  grouping  arguments  to  match 
the  calls  previously  described.   Thus  one  might  set  up  for  the 
punching  of  FORTRAN  statements  by  the  call 

CALL  LSO0UT(4HUNIT,  5LPUNCH,  6HMARGIN,  7, 
*    4HEDGE,  72,  5HCFLAG,  6,  IH*) 

These  calls  to  LSQ0UT  and  stacks  of  code  within  !„SQPNT  may 
be  used  to  compose  special  output  formats  for  particular  expression 
symbols  used  as  function  names.   For  this  purpose,  the  user  must 
be  aware  of  the  following  C0MM0N  block  definitions  within  LSQPNT 

C0MM0N/ LSQNTR/NTRYPT (10) 
C0MM0N/ LSQARG/ LARG (10) 
C0MM0N/LSQINT/INI 
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Assuming  these  definitions  have  been  made,  then  the  array  NTRYPT 
will  contain  locations  which  the  user  may  VISIT  to  output 
expressions,  terms,  or  constants.   The  array  LARG  will  be  used 
to  hold  the  lists  or  list  fragments  being  processed.  The 
variable  INI  is  a  trap  location  used  in  defining  which  expression 
symbols  are  to  be  handled  by  user  code. 

LSQPNT(INI,  symh ,    loa)   defines  the  location  loa ,   established 
in  an  ASSIGN  statement,  to  be  the  place  for  LSQPNT  to  visit  when 
outputting  the  expression  symbol  syrrb   as  a  function  name.   When 
the  visit  is  made,  LARG (6)  will  contain  a  sequence  reader  which 
points  to  the  right  to  the  first  argument  of  symb. 

The  user  may  use  the  code  to  print  an  expression  within  LSOPNT 
while  within  his  own  code  by  storing  the  expression  in  LAHG(l) 
and  performing  a  VISIT  to  either  NTRYPT(2)  or  NTRYPT(3),  the  former 
to  print  the  expression  without  consideration  as  to  whether  it  is 
a  common  subexpression  or  not,  and  the  latter  to  take  this  factor 
into  consideration.   The  code  to  output  a  constant  may  be  used 
by  storing  the  floating  point  constant  in  LARG (4)  (without  mode 
conversion)  and  performing  a  VISIT  to  NTRYPT(6). 

The  user  returns  control  to  LSQPNT  by  a  call  to  TF.RM(O). 
Thus  for  example,  the  following  code  would  cause  the  expression 
symbol  2H0R  to  print  as  an  infix  operator  .0R.   (Let  the  user 
beware  of  confusion  between  the  expression  symbol  OHA.0R.B  and 
the  result  of  printing  0R(A,B).) 
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(:0MM0N  AVSL.XflOO) 
C0MM?)N/  I.SONTR/NTRYPT  (10) 
C0MM0N/ LSQARH/ LARG ( 10 ) 
C0MM0N/LSQINT/INI 


ASSIGN    100  T0   LAC 

CALL   LSQPNT(INI,2H0R,LAC) 


100  CALL  STRI)IR(SRQLR(LARf;(6),N),LARG(l)) 

101  CALL  MANY(X(1),LARC(6)) 
CALL  VISIT(NTRYPT(.'^)) 
CALL  LSQUNM(X(1),LARC(6)) 

CALL  STR1)IR(SLQLR(LAR(;(6),N),LARC,(1)) 
TF(N.NL.0)CALL  THRMfO) 
CAF-L  LSn0UT(4,4ll.0R.) 
G0  T0  101 

It  should  be  noted  that  the  sequence  reader  in  LARG (6)  had  to 
be  saved  before  using  NTRYPT(3). 

LSOPNT(INI,  8ymh,    0)  revokes  any  prior  selection  of 
special  code  to  output  the  expression  symbol  symb   as  a  function 
name. 
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LSQPNT(NTRYPT(1))  initializes  LSOPNT  sufficiently 
so  that  a  user  may  store  an  expression  in  LARG(l)  and  a  Hollerith 
quantity  in  LARC(2)  and  then  VISIT  NTRYPT(l)  to  completely 
simulate  a  call  to  LSQPNT(LARG(1) ,LARn(2))  even  if  the  code 
is  being  executed  within  such  a  call.   This  allows  definitions 
for  LSQVAL  to  print  expressions,  yet  for  special  printing  definitions 
to  use  evaluation.   (See  §3.6). 

There  is  a  function,  PRIPUT,  which  is  used  in  some  existing 
code  to  output  expressions. 

PRlPUTihol,    Us)    is  equivalent  to  a  call  to  LSQPNT(He,  hot) 
This  routine  should  not  be  used  in  new  code. 
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Exer'(ri.se 

3.3.1.   Write  a  program  which  will  punch  out  the  gradient  of  an 
expression  as  valid  FORTRAN  statements  with  sequential  statement 
numbers,  *  in  column  6  for  continuation,  and  the  original  system 
of  coordinates  changed  to  X(l),  X(2),  etc. 
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3.4.   Niceties  of  List  and  Expression  Destruction 

In  addition  to  the  methods  described  in  §2.4,  and  the  destructive 
effects  of  some  of  the  routines  in  §3.1,  the  functions  DELETE,  IRALST, 
MTLIST,  P0PB0T,  P0PMID,  P0PT0P,  and  RCELL  unbuild  lists  and  expressions. 
Normally,  none  of  these  functions,  as  well  as  LSQUNM,  CPYTRM,  INLSTL, 
INLSTR,  NUISTL,  and  NULSTR,  should  be  used  to  unbuild  an  expression. 
Indeed,  other  than  erasure  by  a  call  to  LSQDES,  no  change  should 
ever  be  made  to  an  existing  expression  without  great  care  being 
taken,  i.e.  almost  never.   The  reason  for  such  caution  is  that 
most  of  the  expression  manipulation  routines,  including  LSOVAL, 
LSQADD,  and  LSQMEX,  freely  incorporate  sublists  of  the  expressions 
they  accept  as  input  into  the  expressions  they  create  as  output. 
Thus  any  change  to  a  sublist  of  one  expression  could  well  cause  an 
undesired  change  in  another  expression.   An  even  more  extreme  case 
is  presented  by  expressions  used  in  calls  to  LSQDEF  and  L00KUP 
which  may  themselves  appear  in  their  entirety  within  created  expressions, 
An  expression  which  the  user  himself  creates  by  use  of  LSQMNL,  LIST, 
and  INLIST  (without  the  evaluative  three  argument  call)  may  of 
course  be  manipulated  freely  prior  to  their  use  as  arguments  of  true 
expression  manipulation  routines.   Otherwise,  if  the  user  wishes  to 
alter  an  expression,  he  must  content  himself  to  work  with  a  fresh 
copy  created  by  LSSCPY  (not  LSQCPY)  in  its  place. 

A  user  with  some  sort  of  death  wish,  who  absolutely  must 
take  apart  originals  and  not  copies,  may  prevent  the  use  of  sublists 
of  existing  expressions  in  created  expressions  by  using  the  C0MM0N 
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block  LSQCSX  to  switch  the  function  LSOCXP  from  returning  its 
argument  to  returning  LSSCPY  of  its  argument.   (See  §3.1). 

Dr.LE'TEiquan)    returns  the  list  element  which  would  be 
reached  by  a  right  advance  of  the  quantity  quart   treated  as  a 
sequence  reader,  and  removes  that  element  from  tlie  list.   It  is 
an  error  to  attempt  to  apply  DELF.TE  to  the  head  of  a  list. 

IRALST(Z-te)  reduces  the  count  of  the  number  of  times 
the  list  or  expression  lis   appears  as  a  sublist,  and  returns  this 
reduced  count  as  its  value.   If  the  reduced  count  is  zero,  its 
is  also  erased.   Unless  the  list  lis   was  created  by  a  call  to 
LIST(O)  this  is  the  wrong  way  to  erase  it.   Calling  IRALST  twice 
for  the  same  list  is  a  disaster. 

MTLIST(Zt8)  removes  all  the  elements  from  the  list  or 
expression  lis,    and  returns  the  now  empty  list  lis   consisting  of 
just  a  list  head  as  its  value. 

P0PB0T{Zis)  returns  the  bottom  (rightmost)  element  of 
the  list  or  expression  lis   as  its  value,  and  removes  that  element 
from  the  list.   It  is  an  error  to  apply  P0PB0T  to  an  empty  list. 

P0PMID(yar)  returns  the  list  element  on  top  of  which 
the  variable  or  array  element  var   is  "sitting"  treated  as  a 
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sequence  reader,  removes  that  element  from  the  list,  and  advances 
vaj*  to  the  left.   Thus  the  next  ri^ht  advance  of  var   will  bring 
it  to  the  same  element  as  it  was  pointing  to  on  the  right  originally. 
It  is  an  error  to  apply  P^PMID  in  such  a  way  as  to  attempt  to 
reirove  the  head  of  a  list.   For  example,  var   should  never  be  a 
newly  created  sequence  reader. 

P0PT0P(li3)  returns  the  top  (leftmost)  element  of  the 
list  or  expression  lis   as  its  value,  and  removes  that  element  from 
the  list.   It  is  an  error  to  apply  P0PT0P  to  an  empty  list. 

RCELL(Z-ts)  where  lis   is  an  empty  list  causes  the  absolute 
and  total  erasure  of  the  list,  no  matter  how  many  lists  contain 
it  as  a  sublist.   Thus 

CALL  MTMST(lis) 
CALL  RCELl(lis) 

will  erase  lis   despite  any  degree  of  protection.   This  is  for 
emergencies  only. 
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3.5.   Further  Arithmetic  Operations  on  Expressions 

The  functions  ADD,  DVS'JM,  NIIMPY,  SUB,  and  SUMPY  were 
used  in  some  existing  code  to  perfonn  arithmetic  operations  on 
expressions.   Their  use  in  new  code  should  be  strenuously  avoided. 

ADDCZtsj,  tisi)   returns  the  result  of  adding  expression 
tie  I    to  expression  tis2-     The  first  expression  is  destroyed  and  the 
second  is  replaced  by  the  sum. 

nVSLIM(Hsi,  synb ,    lisi)   adds  the  result  of  dividing  expression 
lisi   by  expression  symbol  sy^i   to  expression  iis2>  replacing  lis2   with 
this  result  which  is  returned  as  the  value  of  DVSUM. 

NUMPY(Zis,  qium)   replaces  the  expression  lis   with  the  result 
of  multiplying  this  expression  by  the  floating  point  constant  quan, 
and  returns  the  product  as  its  value, 

SUB(ZtSi,  li82)   returns  the  result  of  subtracting  expression 
Ubi    from  expression  Zis2-  The  first  expression  is  destroyed  and  the 
second  is  replaced  by  the  difference. 

SUMPY(Z-t«i,  lisi,    lisi,   symbi,  qium\ ,   ...,  symk,,  auan,),  k   = 
0,  1,  2,  ....  5,  returns  the  result  of  adding  the  product  of  expressions 
tisi   and  ItS2  ^°   expression  lis-^   and  truncating  on  the  expression 
symbols  synb  .   to  the  powers  auan.. 
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3.6.   Definition  and  Evaluation  Revisited 

The  description  of  the  functions  LSQDEF  and  LSQVAL  presented 
in  §2.7  was  limited  to  the  simplest  forms  of  definition  and 
evaluation  of  expressions.   The  actual  mechanism  is  rather  complex. 
The  functions  which  are  available  for  definition  and  evaluation  are 
LSQDEF,  LSQDSF,  LSQGAR,  LSQSBS,  LSQVAL,  and  LSQVVL.   Some  existing 
code  contains  references  to  the  following  older  routines  for 
evaluation,  which  should  not  be  used  in  new  code:   EVALUE,  HITENT, 
INSBST,  INSUBT,  INTENT,  SBST,  and  SUBT. 

One  of  the  major  concepts  used  in  evaluation  by  LSQVAL  is  that 
of  "level".   Arguments  of  functions,  expressions  used  as  exponents, 
and  definitions  of  expression  symbols  are  considered  to  lie  on  a 
level  one  deeper  than  that  of  the  expression  on  which  they  were 
discovered.  The  surface  level,  the  level  at  which  evaluation  bepins, 
is  level  1;  the  arguments  of  functions  in  the  expression  being 
evaluated  lie  on  level  2;  arguments  of  functions  within  those 
arguments,  on  level  3;  the  definition  of  an  expression  symbol  within 
these  last  arguments,  on  level  4;  etc.   Evaluation  is  performed 
recursively,  in  a  left  to  right  scan,  always  evaluating  on  level  n   +  1 
before  evaluating  on  level  n.   Naturally,  when  a  constant  expression 
is  encountered,  it  is  returned  as  its  own  value.   When  an  undefined 
expression  symbol  without  arguments  is  encountered,  it  too  is  returned 
as  its  own  value.   However,  when  an  undefined  expression  symbol  with 
arguments  is  encountered,  its  arguments  are  evaluated  and  the 
expression  symbol  applied  to  the  evaluates.  The  user  may  set 
limits  on  the  depth  to  which  evaluation  will  be  carried.   Beyond 
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such  a  limit,  evaluation  will  consist  of  simple  copying. 

LSQVAL(Hs,  int)   returns  the  result  of  evaluating 
expression  lis   through  level  int,   an  integer  quantity.   If  int   is 
negative,  no  limit  is  in  effect.   If  int   is  zero,  a  copy  of  the 
expression  lis   is  returned.   This  copy  may  have  common  subexpressions 
with  lis.      The  one  argument  call  to  LSQVAL(iis)  is  equivalent 
to  LSQVAL(Hs,  -1),  i.e.  evaluation  on  all  levels  as  in  §2.7. 

If  simple  substitutions  of  expressions  for  expression  symbols, 
without  any  consideration  of  arguments,  are  the  only  definitions 
in  effect,  then  limits  on  the  level  of  evaluation  will  not  result 
in  great  surprises.   However,  if  dummy  arguments  get  involved 
they  may  well  appear  in  the  evaluate  if  the  limit  level  invokes 
the  definition,  since  the  dummy  arguments  could  not  be  evaluated 
when  encountered  beyond  the  limit. 

Since  the  concept  of  level  does  work  well  in  simple  substitutions, 
the  function  LSQSBS  provides  for  its  application. 

LSQSBS(Hsi,  symb,    lis2,   int)   returns  the  result  of 

substituting  expression  lis i    for  the  expression  symbol  symb   in 

the  expression  lis2   through  the  level  specified  by  the  integer 
quantity  int. 

In  expressions  where  it  is  a  priori  known  that  beyond  a  certain  level 
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the  expression  symbol  for  which  a  substitution  is  to  be  made  does 
not  appear,  the  user  can  greatly  reduce  the  cost  of  substitution 
by  using  that  level  as  a  limit.   If  the  limit  level  is  1,  the 
expression  symbol  for  which  substitution  is  being  made  may  appear 
safely  within  the  expression  to  be  substituted  for  the  expression 
symbol  without  causing  an  infinite  loop,  since  the  evaluation  will 
reduce  to  a  copy  as  soon  as  the  substitution  is  made. 

The  definitions  used  by  LSQVAL  also  involve  a  concept  of 
level.   The  same  expression  symbol  may  be  defined  on  many  different 
definition  levels  simultaneously,  but  only  once  on  each  level. 
Where  we  speak  of  the  depth  of  a  level  within  an  expression,  we 
shall  speak  of  the  height  of  a  definition  level.  The  lowest 
definition  level  is  zero.   This  is  the  level  on  which  the  calls 
to  LSQDEF  described  in  §2.7  introduce  definitions.  There  are 
further  calls  to  LSQDEF  which  allow  definitions  on  level  zero  in 
terms  of  FORTRAN  code  to  VISIT,  and  definitions  on  higher  levels. 
At  any  given  time,  the  user  may  make  a  definition  on  level  zero, 
introduce  a  new  highest  level  of  definition,  make  a  simple 
substitution  definition  on  the  highest  level,  discard  all  the 
definitions  on  the  highest  level  (provided  that  level  is  not  zero) 
and  regress  to  the  next  definition  level  down,  suspend  the  use 
of  level  zero  definitions  and  the  right  to  introduce  level  zero 
definitions,  or  return  level  zero  to  full  use.   Definitions  on  a 
given  level  override  definitions  on  any  lower  definition  level. 
Levels  higher  than  zero  are  largely  for  the  internal  use  of  LSQVAL. 
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I,SQ[)PF(0,  5HBEGIN)  introduces  a  new  highest  level  of 
definition. 

LSQDUF (s ymb ,    lis)   defines  the  expression  symbol  symb   to 
be  the  expression  lis.      This  definition  is  made  on  the  highest 
current  definition  level.   It  replaces  any  other  definition  of  symb 
on  the  same  level,  and  overrides  any  other  definition  of  symb   on  a 
lower  level.   If  the  expression  lis   is  simply  a  representation  of 
an  expression  symbol  symb',   then  before  the  definition  is  made  lis 
is  replaced  by  the  definition  of  symb',    if  any,  found  on  the  highest 
non-zero  level  containing  such  a  definition.   An  attempt  to  define 
an  expression  symbol  to  be  itself  will  be  ignored.   An  attempt  to 
define  an  expression  symbol  to  be  an  expression  involving  itself 
will  be  accepted  and  may  cause  an  infinite  loop  (or  loop  through 
the  limit  level).   Whether  a  loop  occurs  or  not  depends  on  how 
the  expression  symbol  appears  in  its  definition.   If  such  an 
appearance  is  reached  in  the  course  of  evaluation,  a  loop  will 
occur. 

The  expression  lis   should  be  considered  to  be  erased  by  this 
call  to  LSQDEF. 

LSQDEP (symb,    0)  "defines"  the  expression  symbol  symb   to 
be  undefined.   This  action  is  taken  on  the  highest  current  level 
of  definition.   It  replaces  any  definition  of  symb   on  the  same 
level,  and  overrides  any  other  definition  of  symb   on  a  lower  level. 
As  long  as  the  level  on  which  this  un-definition  was  made  exists 
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and  no  new  definition  of  symb   is  made  on  a  higher  level,  the  expression 
symbol  symb   acts  as  if  it  were  totally  undefined  despite  any  old  or 
new  definitions  made  on  level  zero. 

LSQDEF(0,  3HEND)  discards  the  current  highest  level  of 
definition  and  all  definitions  made  thereon.   The  next  definition 
level  down  becomes  the  current  highest  level  of  definition.   It  is 
an  error  to  attempt  to  discard  level  zero  in  this  manner. 

Though  it  is  not  necessarily  an  error  to  do  otherwise,  the 
four  two-argument  calls  to  LSQIJEF  described  above  should  normally 
be  used  in  the  following  order: 


CALL  LSODEFCO,  SHBPniN) 

CALL  LSQl)E¥ (symb I,    lisi) 

any  other  substitution  definitions 
CALL  LSQDEF  (si/mfc,  ,  0) 

any  other  "un- definitions" 
CALL  ISQVMilis) 

any  other  evaluations  using  these  definitions 

CALL  LSQDEF(0,  3HENU) 
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That  is  essentially  the  technique  used  within  LSOVAL  when  it 
encounters  a  level  zero  definition  with  dummy  arguments.   The  dummy 
arguments  are  defined  to  be  the  actual  arguments  and  the  definition 
is  evaluated.   Any  dummy  arguments  which  are  not  matched  by  actual 
arguments  are  left  unmatched,  causing  any  other  existing  definition 
of  the  unmatched  dummies  to  take  effect. 

ISQDEF {symb,    lisi,    H32)   defines  the  expression  symbol 

symb   to  be  the  expression  Hsi    in  which  those  expression  symbols 

which  appear  on  the  list  li32   are  dummies  for  any  actual  arguments  to 

which  symb   may  be  applied.   This  definition  is  made  on  definition 

level  zero,  and  replaces  any  prior  definition  of  symb   on  this 

level.   The  list  lis2    is  of  the  form  (symbi,  symb2,    ...,   symb,) 

for  fe  =  0,  1,  2,  ...«.  The  expression  symbol  symb.   is  used  as 

t 

a  dummy  for  the  i        actual  argument  by  calls  to  LSODIiF  from  LSQVAL 
as  above.  Thus  the  user  may  think  of  this  call  as 

3ymb{symb\,   symb2,    ....  symb,)   >  lisi 

Though  use  of  such  a  definition  with  fewer  actual  arguments  than 
dummies  raises  the  problems  described  above,  having  more  actual 
arguments  than  dummy  arguments  yields  reasonable  results.   Any 
residual  actual  arguments  are  evaluated  and  appended  to  the  argument 
lists  of  all  expression  symbols  in  the  evaluate  of  the  definition, 
which  lie  on  level  1  of  the  evaluate. 

The  expression  lisi   and  the  list  lie2   should  be  considered  to 
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be  erased  by  this  call. 

LSQDEFisymb,    lis,   0)  defines  the  expression  symbol  symb 
to  be  the  expression  lis.     There  are  no  dummy  arguments  to  be 
considered.   This  definition  is  made  on  definition  level  zero, 
but  is  otherwise  similar  to  calling  LSQDE¥ {symb ,    lis).      Again  the 
user  must  consider  the  expression  lis   to  be  destroyed  by  this  call. 

LSQDEF (symb ,   0,  0)  revokes  any  prior  definition  of  the 
expression  symbol  symb   made  on  definition  level  zero  and  leaves 
symb   undefined  on  this  level. 

LSQDEF (symb ,  li8\,    lisi,   5nN0CUR)  defines  the  expression 
symbol  symb   to  be  the  expression  lis\   with  dummy  arguments  as 
specified  on  the  list  lis2-      This  call  differs  from  the  similar 
three  argument  call  only  in  that  when  the  definition  expression 
is  evaluated  LSQDEF(s.vmi),  0)  will  be  called  to  make  the  expression 
symbol  symb   appear  to  be  undefined.   This  permits  symb   to  appear 
within  the  expression  lis\   without  causing  infinite  loops. 

LSQDEF (si/mfc,  lis,   0,  5HN0CUR)  acts  as  the  call  above 
except  that  there  are  no  dummy  arguments  involved. 

LSQDEF(s!ymib,  loc ,   int)    defines  the  expression  symbol 
symb   to  he  the  result  of  a  VISIT  to  the  location  loc,    established 
in  an  ASSIGN  statement.   The  integer  quantity  int   determines  the 
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handling  of  arguments  and  whether  previously  computed  evaluates 
of  the  same  function  call  may  be  used  if  available.   If  each  time 
the  code  at  location  loo   is  applied  to  the  same  arguments,  the 
same  result  is  returned  and  no  changes  are  made  to  variables, 
arrays,  lists,  expressions,  etc.  which  might  be  used  elsewhere, 
then  the  integer  quantity  int   should  be  negative.   Otherwise  int 
should  be  positive.   If  the  code  at  too   should  be  presented  with 
unevaluated  arguments,  then  int   should  be  of  magnitude  1.   If  the 
code  should  have  evaluated  arguments,  then  int   should  be  of 
magnitude  2.   Thus  the  four  possible  values  for  the  integer  quantity 
int   are: 

-2  no  side  effects,  evaluate  arguments 

-1  no  side  effects,  do  not  evaluate  arguments 

1  side  effects,  do  not  evaluate  arguments 

2  side  effects,  evaluate  arguments 

The  user  is  again  reminded  that  all  unspecified  values  are  reserved 
for  future  expansion  of  the  system. 

The  evaluation  code  at  location  loc   communicates  with  LSQVAL 
via  the  C0MM0N  block  LSQDMF ,  defined  by 

C0MM0N/LSODMF/ LSQDMF (15) 

where  LSQDMF(13)  will  contain  a  sequence  reader  which  may  be  advanced 
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to  the  right  to  reach  the  first  argument,  to  the  right  again  to 
reach  the  second  argument,  and  so  forth.   LSQDMF(14)  will  contain 
the  expression  symbol  which  invoked  the  VISIT.   LSQDMF(15)  will 
contain  the  number  of  arguments.   An  expression  may  be  evaluated 
by  storing  the  expression  in  LSQDMF(ll)  and  paying  a  VISIT  to 
LSQ[)MF(6).   An  integer  one  less  than  the  level  within  the  expression 
being  evaluated  on  which  the  VISIT  to  loa   was  invoked  may  be  found 
in  LSQDMF(8).   The  limit  level  for  evaluation  will  be  an  integer 
in  LSQDMF(9). 

A  value  lie   should  be  returned  to  LSOVAL  by  a  call  to  TFRM(Hs). 
The  value  must  be  a  valid  expression,  but  need  not  be  a  newly 
created  expression.   LSODES  will  be  applied  to  lis   to  erase  it  if 
it  was  newly  created. 

To  aid  in  writing  code  to  be  visited  in  an  evaluation,  the 
functions  LSQDSF,  LSQGAR,  and  LSQWl,  are  provided. 

LSQ[)SF(/i/n,  symb)   returns  the  result  of  applying  the 
FORTRAN  callable  function  fun,   which  should  be  declared  in  an 
EXTERNAL  statement,  to  the  next  expression,  if  any,  to  be  reached 
by  a  right  advance  of  LSQ1)MF(13).   If  the  expression  found  is 
itself  a  constant  expression,  the  result  of  calling  LSQDSF  will  be 
a  constant  expression.  Otherwise  an  expression  will  be  formed  by 
using  the  expression  symbol  symh   as  a  function  name  and  the 
expression  found  by  the  right  advance  of  LS0DMF(13)  as  its  argument. 
For  example,  calling  the  following  subroutine  would  provide 
definitions  of  the  arcsine  and  arccosine  functions  using  the 
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expression  symbols  611ARCSIN  and  6HARCC0S  respectively. 

SUBR0UTINE  LSQASC 

EXTHRNAL  AS  IN,  AC0S 

ASSIGN  101  T0  L0C 

CALL  LSQDEF(6HARCSIN,  L0C,  -2) 

ASSIGN  102  T0  L0C 

CALL  LSQDHF(6HARCC0S,  L0C,  -2) 

RETURN 

101  CALL  TERM(LSQDSF(ASIN,  6HARCSIN)) 

102  CALL  TERM(LSODSF(AC0S,  6HARCC0S)) 
END 

LSOGAR(O)  returns  the  next  expression  to  be  reached  by  a 
right  advance  of  LS0DMF(13)  if  there  is  any  next  expression.  Otherwise 
the  value  returned  is  zero. 

LSOGAR (s/ymfc)  returns  the  next  expression  to  be  reached  by 
a  right  advance  of  LSQF)MF(13)  if  there  is  any  next  expression.   Otherwise 
an  error  abort  is  taken  with  the  expression  symbol  symb   provided  as 
an  error  message. 

LSQVVLCHa)  returns  a  location  to  VISIT  in  order  to 
evaluate  the  expression  lis.  The  VISIT  will  save  and  restore 
LS0UMF(13)  so  that  further  scanning  of  arguments  will  be  possible. 
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Consider  the  matter  raised  in  §3.2  concerning  the  use  of  the 
expression  symbol  6HEQUAL.  as  the  translation  of  the  equals  sign 
on  input  via  INLIST.  The  following  subroutine  will,  when  called, 
cause  6HEQUAL.  to  be  defined  as  a  means  of  creating  definitions. 

SUBROUTINE  LSOIEQ 

ASSIGN  100  T0  L0C 

CALL  LSQDEF(6HEQUAL. ,  L0C,  1) 

RETURN 

100  LA=  LSQGAR (6HEQUAL . ) 
NLA=LSQTYP(LA,NAM) 

IF  (  (NLA.  LT.  10000).  0R.  (NLA. r.E.  20000)) 
*    CALL  LS0ERR(6HEQUAL.) 
SA=SEORDR(LA) 
STU=SEORnR (SEQLR (SA ,N) ) 
SAC=SEQLR(STU,N) 
LTU=L1ST(9) 

101  SAR=SEQLR(STU,N) 
IF(N.NE.O)G0  T0  102 
MLA=LSQTYP(SAR,MAM) 

IF(MLA.NE.10000)CALL  LSQERR(6HH0UAL. ) 
CALL  MANY (LTU, MAM) 

C0  T0  101 

102  LB= LSQGAR (6HEnUAL.) 

CALL  LSQUEF(NAM,LSQCPY(SEORnR(LB)) ,LTU) 

CALL  TERM (LB) 

END 


124 


3.41 


Note  that  this  limits  F.QUAl,.  to  two  useful  arpuments. 

We  may  also  wish  to  define  the  expression  symbol  5MPRINT  to 
print  the  value  of  its  first  argument  with  its  second  argument  for 
a  name.  The  following  subroutine  will,  when  called,  so  define  SlIPRINT. 

SUBROUTINE  LSQIPN 
C0^^M0N  AVSL,  X(IOO) 
C0MM(ON/LSQNTR/NTRYPT  (10) 
C0MM0N / LSQARG / LA  RG ( 1 0 ) 
C0MM0N/LSQINT/INI 
CALL  LSQPNT(NTRYPT(1)) 
ASSIGN  100  T0  L0C 
CALL  LSQUEF(5!IPRINT,L0C,1) 
RETURN 
100   LA=LSQGAR(5HPR1NT) 

CALL  STRD1R(VISIT(LS0VVL(LA)),LARG(1)) 

LB=LS0GAR(5HPR1NT) 

I F ( LSOTYP ( LB . LARG (2 ) ) . NE . 10000) CALL  LSQERR (5HPRINT) 

CALL  MANY(X(98) .LARGCl)) 

CALL  VISIT(NTRYPT(1)) 

CALL  LSQUNM(X(98),LA) 

CALL  TERM(LA) 

END 

In  both  these  examples,  the  definitions  involve  side  effects 
and  unevaluated  arguments  are  used.  Actually,  when  "unevaluated" 
arguments  are  requested,  arguments  evaluated  with  definition  level 
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zero  suspended  are  provided.  Usually  this  simply  causes  actual 
arguments  to  be  substituted  for  dummy  arguments.  Thus,  if  both 
the  routines  LSQIEQ  and  LSQIPN  have  been  called,  input  lines  of 
the  form 

PR(A,.1)=PRINT(A..1,A..1)$$ 

LA=1+SIN(X)$$ 

PR(LA)$$ 

will  generate  output  of  the  form 

LA  =  1   +   SIN(X) 
$ 
SEND  0F  EXPRESSI0N 

rather  than 

A,  .1  =  1   ♦  SIN(X) 
$ 
SEND  0F  EXPRESSI0N 

since  PRINT  will  be  invoked  by  the  evaluation  of  PR  with  the  actual 
argument  LA  substituted  for  the  dummy  A..1. 

If  the  user  shoiild  do  all  his  expression  printing  by  means  of 
evaluation  of  an  expression  symbol  such  as  5HPRINT,  then  he  may  use 
the  entries  to  LSQVAL  by  means  of  VISITS  within  special  print  definitions. 
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LSQDEF(0,  4HSAVE)  suspends  all  level  zero  definitions 
and  makes  it  an  error  to  introduce  any  level  zero  definitions. 

LSQDr.F(0,  6HUNSAVf;)  restores  definition  level  zero  to 
full  use. 

It  is  always  proper  to  use  these  two  calls  (in  pairs)  within 
FORTRAN  code  used  as  a  level  zero  definition,  since  in  order  to  have 
reached  the  code,  level  zero  must  not  have  been  suspended  already. 
In  all  other  cases  the  user  must  be  careful  not  to  call  LSni)EF(0,'inSAVF) 
twice  in  a  row. 
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EVALUE(Hs,  symbi,  quani ,    ....  symb,,  qucm,)    for  fe  =  1, 
2,  3,  4,  5,  6,  evaluates  the  expression  lis   with  the  expression 
symbols  symb .   set  to  the  floating  point  constants  guan..     The 
evaluate  is  returned  as  the  value  of  the  function  and  also  replaces 
the  expression  lis.      Inasmuch  as  this  is  a  destructive  manipulation 
of  an  expression,  calls  to  EVALUE  should  be  avoided  in  new  code. 

HITENT(Hsi,  symb,    lis2,    lisj)    returns  the  result  of 
adding  a  level  one  substitution  of  the  expression  Hsj  for  the 
expression  symbol  symb   in  the  expression  Zis2  to  the  expression 
lis2   which  is  replaced  by  this  result.   LSQSBS  is  a  better  choice 
for  new  code.   This  substitution  leaves  parentheses  around  the 
uses  of  lis  I. 

INSBST(ZtSi,  symb,    lis2,    lisj)   returns  the  result  of 
adding  an  all  levels  substitution  of  the  expression  lis\    for  the 
expression  symbol  symb   in  the  expression  lis2   to  the  expression 
lis2   which  is  replaced  by  this  result.   I.SQSBS  is  again  a  better 
choice  for  new  code. 

INSUBT(lisi,  symb,    lisi,    lis^,  symbi,  quani symb,, 

quan,)    for  k  =   I,    2,    3,   4 ,   S ,   returns  the  result  of  adding  an  all 
levels  substitution  of  the  expression  lisi    for  the  expression  symbol 
symb   in  the  expression  lis2   to  the  expression  lis^,   and  then  truncating 
on  the  expression  symbols  symb.   to  the  powers  given  by  the  floating 
point  quantities  quart..      The  expression  lis^   is  replaced  by  the 
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result.   LSQSBS  and  LSQTRC  should  be  used  in  new  code. 

INTENTCHsi,  symb,    lis2,    Us^)   acts  as  IIITENT  except  that 
substitution  is  done  on  all  levels. 

SBSTCHsj,  symb,    lis2,    Hs^)   acts  as  INSBST  except  that 
substitution  is  restricted  to  level  one. 

SUKT(lisi,symb ,    lis2,    tis^,   i^ymh\,  quan\,    ....  symb-,,  quan,) 
for  fe  =  1,  2,  3,  4,  5,  acts  as  INSUBT  except  that  substitution  is 
restricted  to  level  one. 


129 


3.46 


Exevcise 

3.6.1.   Assume  that  the  INPUT  and  0UTPUT  files  are  actually  a 
teletype.  Write  a  program  using  the  subroutines  LSQIEQ  and  LSQIPN 
and  such  additional  code  as  you  may  need  to  form  a  desk  calculator 
for  expressions.   You  would  do  well  to  use  the  three  argument  call 
to  INLIST,  so  that  LSQIEQ  will  have  something  to  work  on.   Consider 
including  a  definition  for  the  expression  symbol  ST0P  which  will 
terminate  execution.   Also  consider  ways  to  allow  the  input  expressions 
to  begin  in  column  1. 
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3.7.  Truncation  Recapped 

The  function  LSQTRC,  mentioned  in  §2.8,  truncates  expressions 
using  the  concept  of  levels  described  for  evalauation  in  §3.6. 

LSOTRC(Hs,  symh ,  quani,  qucm2,   tntj,  int2)   returns  an 
expression  derived  from  expression  lis   by  retaining  certain  terms 
of  the  expression  through  level  int\,   an  integer  quantity,  and  all 
terms  on  deeper  levels.  Whether  a  term  is  retained  or  not  depends 
on  whether  it  contains  the  expression  symbol  symh   to  a  non-constant 
power  or  to  a  power  between  the  floating  point  quantities  quan^   and 
quani • 

If  the  integer  quantity  int2    is  0,  2,  4,  6,  8,  10,  12,  or  14, 
then  terms  containing  the  expression  symbol  symb   to  a  constant 
power  lying  between  quani    and  quan2   will  be  retained.   If,  on  the 
other  hand  int2   is  1,  3,  5,  7,  9,  11,  13,  or  15,  then  terms  containing 
symb   to  a  constant  power  lying  outside  of  the  range  between  auani   and 
quan2   inclusive.   If  int2   is  0,  1,  4,  5,  8,  9,  12,  or  15,  the  terms 
containing  symb   to  non-constant  powers  will  not  be  retained;  while 
if  int2    is  2,  3,  6,  7,  10,  11,  14,  or  15,  terms  with  symb   to  non- 
constant  powers  will  not  be  retained.   If  int2    is  0,  1,  2,  3,  8,  9, 
10,  or  11,  then  terms  containing  symb   to  a  constant  power  equal  to 
the  maximum  of  qucmi    and  quan2   will  not  be  retained,  while  for  int2 
equal  to  4,  5,  6,  7,  12,  13,  14,  or  15,  they  will  be  retained.   If 
int2    is  0,  1,  2,  3,  4,  5,  6,  or  7,  then  terms  containing  symb   to  a 
constant  power  equal  to  the  minimum  of  quani   and  quan2   will  not  be 
retained,  while  for  int2   equal  to  8,  9,  10,  11,  12,  13,  14,  or  15, 
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they  will  be  retained.   A  term  which  does  not  contain  the  expression 
symbol  eymb   (other  than  within  parentheses  or  function  arguments)  is 
considered  to  contain  symb   to  the  constant  power  zero.   Appearances 
of  symb   within  parentheses  or  function  arguments  do  not  affect  the 
decision  as  to  whether  a  term  should  be  retained  or  not.   Appearances 
of  symb   with  its  own  arguments  are  considered.  Terms  are  scanned 
from  left  to  right,  and  the  first  appearance  of  a  factor  which  gives 
grounds  for  dropping  the  term  will  cause  the  term  to  be  dropped. 
LSQTRC  deals  in  terms  of  the  minimum  and  the  maximum  of  the 
floating  point  quantities  auarii    and  auan2,    so  that  they  may  be 
interchanged  without  changing  the  result. 

As  in  §2.8,  LSQTRC  may  be  called  with  fewer  than  six  arguments. 
This  will  cause  default  values  to  be  used  for  quani,   quan2,   int\,   and 
inti   when  the  call  fails  to  specify  them.   The  arguments  lis   and  symb 
must  always  be  provided.   The  default  values  are  as  follows: 

2  argument  call  quan^-=0,  quan2=0 ,   inti  =  l,   int2  =  12 

3  argument  call  quan2=0 ,   int\=\ ,   int2=l2 

4  argument  call  inti=l,   int2=12 

5  argument  call  int2=12 

Thus  in  these  defaults,  the  expression  lis   will  be  truncated  to 
constant  powers  within  a  closed  interval. 
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The  user  may  find  it  convenient  to  think  of  the  truncation  mode 
specified  by  int2    in  terms  of  the  interior  or  exterior  of  an  interval 
with  or  without  the  left  and  ripht  endpoints.   The  mode  is  then  the 
sum  of  four  quantities: 

interior  0  for  the  interior  of  the  interval 

1  for  the  exterior  of  the  interval 

constant  0  for  expression  exponents  not  retained 

2  for  expression  exponents  retained 

right  open         0  to  exclude  the  ripht  endpoint  of  the  interval 
4  to  include  the  right  endpoint  of  the  interval 

left  open  0  to  exclude  the  left  endpoint  of  the  interval 

8  to  include  the  left  endpoint  of  the  interval 

Thus  the  default  value  of  int2,    12,  specifies  the  interior  with  I^oth 
endpoints  and  no  non-constant  exponents: 


while  setting  int2   to  7  specifies  the  exterior  with  non-constant 
exponents  and  the  rijjht  endpoint,  but  not  the  left: 
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TRUNC(Zts,  symb,   quan)    is  a  truncation  routine  used  in 
some  existing  code,  which  should  not  be  used  in  new  code.   This 
routine  replaces  the  expression  lis   with  the  result  of  discarding 
all  first  level  terms  which  contain  the  expression  symbol  symb   to 
a  power  greater  than  the  floating  point  quantity  quan,   or  to  a 
non-constant  power. 
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3.8.   Further  Means  of  Analysis  of  Lists  and  Hxpressions 

There  are  many  more  facilities  for  examining  lists  and  expressions 
than  were  covered  in  §2.10.   The  functions  ADVLF;L,  ADVLFR,  AUVI.NL, 
ADVLNR,  ADVLWL,  ADVLWK,  ADVSF.L,  ADVSF.R,  ADVSNL,  ADVSNR,  ADVSWL,  ADVSWR, 
which  we  represent  by  "ADVaBy".  B0T,  BRKAK,  GETC0F.,  INITRD,  IRARDR, 
LCNTR,  L0CT,  L0FRDR,  LPNTR,  LRDPCP,  LRL)R0V ,  LVLRVl ,  LVLRVT,  MAPLFT, 
MADNBT,  MADNTP,  MAURGT,  NAMTST,  POWFR,  PFHD,  SF.QSL,  SEOSR,  S0r,VF.,  T0P, 
TRCAL,  and  TSTC0N  may  be  so  used. 

B0T(Hs)  returns  the  bottom  friphtmost)  element  of  the  list 
or  expression  lis,   which  is  assumed  to  be  non-empty.  The  element 
returned  remains  on  the  list. 

l<fiV{lis)    returns  the  top  (leftmost)  element  of  the  list  or 
expression  lis,   which  is  assumed  to  be  non-empty.  The  element  returned 
remains  on  the  list. 

If  LA  is  a  non-zero  expression,  then  LSOMNL(T0P(LA))  and  LSOMNL( 
B0T(LA))  are  properly  formed  expressions  consisting  of  the  leftmost 
and  rightmost  terms  of  LA  respectively. 

L0CT(Hs)  returns  the  list  or  expression  lis,   provided  lis 
is  actually  a  list.   Otherwise  L0CT  forces  an  error  abort.   This 
routine  is  used  by  subroutines  to  insure  that  arguments  i#iich  are 
supposed  to  be  lists  actually  are  lists. 
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BREAK{ltSi,  symb ,    lis2,    lis^)   adds  to  the  expression  lis2 
those  terms  of  the  expression  lisi   which  contain  the  expression  symhol 
symb   on  the  first  level,  replacing  Hs2  with  this  result.  The  remaining 
terms  of  lisi   are  added  to  the  expression  lis^,   which  is  changed  to 
this  new  value.  This  routine  should  not  be  used  in  new  code.   Calls 
to  LSOTRC  are  preferable. 

GETC0E (symb ,  quart,    lis\,    Z1S2)  adds  to  the  expression  lis2 
the  coefficient  of  the  expression  symbol  symb   raised  to  the  floating 
point  constant  power  quan   in  the  expression  lis\,   replacing  li62 
with  the  result  and  returning  it  as  the  function  value.  This  routine 
should  not  be  used  in  new  code. 

NAMTST ((7i«3n )  returns  zero  only  if  the  quantity  queen   is  a 
list  or  expression,  -1  otherwise. 

P0WER(Hs,  symb)   returns  the  leftmost  power  to  which  the 
expression  symbol  symb   is  raised  in  the  expression  tic,    considering 
only  the  first  level.   Zero  is  returned  if  r^ymh   is  not  found. 

S0LVr(ZtSi,  symb,    lis2)    assumes  that  the  expression  lifti 
is  linear  in  the  expression  symbol  r:ymb   and  solves  tlie  equation 
lts\    =   0  for  symb.      The  solution  is  added  to  the  expression  lir,2 
which  is  replaced  by  the  sum.   The  sum  is  also  returned  as  the 
function  value.   This  routine  should  not  be  used  in  new  code. 
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TRCAL(3ts,  symh\,   symby)    finds  the  lowest  power  to  which 
the  expression  symbol  nymh2    is  raised  in  those  terms  of  the  expression 
lis   wliich  contain  the  expression  symbol  symbi.      Then  the  expression 
lis    is  replaced  by  a  truncated  copy  in  which  those  terms  that  contain 
the  expression  symbol  symh\    to  a  power  greater  than  that  minimal 
power  of  symb2   are  discarded.   Tliis  routine  should  not  be  used  in 
new  code. 

TSTC0N(its,  var)    returns  zero  if  the  expression  lis   is 
not  a  constant  expression.   If  lis   is  a  constant  expression,  the 
value  returned  is  1.0,  and  the  variable  or  array  element  var   is 
set  to  the  equivalent  I'ORTRAN  constant  value.   LSQTYP  should  be 
used  in  new  code. 

SEQSL(yar,  intvar)   returns  the  next  list  element  otlier 
tlian  a  sublist  encountered  in  left  advances  of  the  variable  or 
array  element  var*   treated  as  a  sequence  reader,  and  by  descents 
into  any  sublists  encountered.   Tlie  integer  variable  or  array 
element  intvar   will  be  set  to  1  if  a  list  head  stops  tlie  search, 
and  -1  if  a  list  element  other  than  a  sublist  or  list  head  is  found. 
SF.QLI,  and  SF.QLR  may  be  applied  to  the  sequence  reader  var   after  use 
by  SEQSL,  but  the  user  should  be  aware  that  ascent  back  to  higher 
level  lists  is  not  possible. 
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SEQSRCuar,  intvar)   returns  the  next  list  element  other  than 
a  sublist  encountered  in  right  advances  of  the  variable  or  array 
element  var   treated  as  a  sequence  reader,  and  by  descents  into  any 
sublists  encountered.  The  integer  variable  or  array  element  intvar 
will  be  set  to  1  if  a  list  head  stops  the  search,  and  -1  if  a  list 
element  other  than  a  sublist  or  list  head  is  found.  The  same 
caution  as  to  subsequent  use  of  SCQLL  and  SLQLR  hold  as  for  SEQSL. 

As  an  example,  consider  the  following  code  assuming  LA  to  be 
an  expression. 

SA=SEQRDR(LA) 
SC0=SF,QSR(SA,N) 

If  LA  represents  [[0]]  then  \'  will  be  1;  otherwise  N  will  be  -1, 
SC0  will  be  the  coefficient  of  the  leftmost  term  of  LA,  and  SA 
will  be  a  sequence  reader  for  that  term  (not  for  LA)  which  is 
ready  to  advance  to  the  right  to  the  first  factor. 

MAnLFT(r/Kan)  returns  a  quantity  wliich  may  be  used  as  a 
sequence  reader  for  a  right  advance  to  the  list  clement  to  the 
left  of  the  list  element  whicli  would  be  reached  by  a  right  advance 
of  the  quantity  quan   treated  as  a  seciuence  reader.   Thus  the  list 
element  on  whicli  quan    is  sitting  will  be  reached  I'y  a  right  advance 
of  the  result. 
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MADRGT (c^uan )  returns  a  quantity  wliich  may  be  used  as  a 
sequence  reader  for  a  ripht  advance  to  the  list  element  to  the 
right  of  the  list  element  which  would  be  reached  by  a  rip.ht  advance 
of  the  quantity  quan   treated  as  a  sequence  reader. 

MADNBT(Hs,  int)    returns  a  qtiantity  which  may  be  used  for 
a  right  advance  to  the  element  of  list  or  expression  tin   whicli  is 
int   elements  to  the  left  of  the  head  of  lis,   where  int    is  a  non-zero 
positive  integer  quantity.   If  there  are  fewer  than  int   elements  in 
lis,    then  the  search  continues  around  the  list,  counting  the  head 
as  an  element  each  time  it  is  encountered. 

MAI)NTP(ZtG,  int)    returns  a  quantity  which  may  lie  used  for 
a  right  advance  to  the  element  of  the  list  or  expression  lis   which 
is  int   elements  to  the  right  of  tlie  head  of  lis,   where  int   is  a 
non-zero  positive  integer  quantity. 

The  remaining  routines  discussed  in  this  section  are  principally 
concerned  with  list  readers,  a  more  flexible  tool  of  list  scanning 
than  sequence  readers,  in  that  a  list  reader  may  descend  and  ascend 
within  a  list  and  its  sublists  easily.   A  concept  of  list  level  is 
used  in  conjunction  with  list  readers  that  is  simply  related  to  level 
within  an  expression.   A  level  in  an  expression  lies  on  two  list 
levels:   the  one  on  wiiich  the  terms  of  the  expression  are  encountered, 
and  the  one  on  which  the  term  coefficient  and  factors  lie.   Term 
coefficients  on  expression  level  n    lie  on  list  level  2n-l. 
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LRDR0V(Zis)  returns  a  newly  created  list  reader  for  the 
list  or  expression  lis.      A  list  reader  is  itself  very  similar  to  a 
list,  and  draws  on  the  same  available  space  list  for  a  head  and 
for  elements  to  use  to  hold  enough  information  to  ascend  from 
sublists  into  which  it  descends,  as  do  lists  and  expressions. 

IRARDR((7uan)  erases  the  list  reader  quart,    and  returns  as 
it  value  the  list  level  the  reader  had  reached  when  erased.   List 
readers  should  be  erased  when  not  needed  so  that  their  cells  will 
be  available  for  use  in  lists  and  other  list  readers. 

kDVa&y  (quan ,   var) ,   where  a  is  I,  or  S,  B  is  W,  F,  or  N,  and 
Y  is  I.  or  R,  returns  the  next  list  element  reached  by  an  advance  of 
the  list  reader  specified  by  the  quantity  quan   in  the  direction  y 
(L  for  left,  R  for  ripht)  either  restricted  to  the  current  list 
level  (a  =  L)  or  descending  into  and  rising  from  sublists  when 
encountered  (a  =  S) ,  and  considering  the  searcli  satisfied  either  by 
a  list  head  or  a  list  element  specified  by  B  as  follows:   6  is  W  for 
the  next  list  element  (whether  or  not  a  sublist),  F.  for  the  next 
clement  which  is  not  a  sublist,  and  N  for  the  next  list  element  which 
is  a  sublist.  The  variable  or  array  element  var   is  set  to  -I.O  if 
the  search  is  stopped  by  a  list  head,  0.0  otherwise. 

l^iUrSiquan)   returns  tlie  list  reader  wliich  the  quantity 
quan   represents  and  advances  that  reader  linearly  to  the  head  of 
the  list  or  sublist  on  the  list  level  being  scanned. 
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I.CNTR Cqwcm )  returns  the  list  level  to  which  the  list 
reader  specified  by  the  quantity  quan   has  descended.   LCNTR  may 
also  he  applied  to  a  list,  in  whicli  case  the  number  of  times  that 
list  appears  as  a  sublist  will  be  returned. 

\.0\'Kr>R{quan)    returns  the  list  or  sublist  currently  being 
scanned  by  the  list  reader  specified  by  the  quantity  quan. 

LPNTR(quan)   returns  a  quantity  which  may  be  used  as  a 
sequence  reader  for  a  right  advance  to  the  list  element  to  which 
tlie  list  reader  specified  by  the  quantity  qi<cm   currently  refers, 
i.e.  the  element  to  which  it  last  advanced. 

LRl)\lCP(quan)    returns  a  copy  of  the  list  reader  specified 
by  the  quantity  quan.      The  copy  will  be  in  the  same  state  as  the 
original  with  respect  to  the  elements  it  would  reach  on  the  next  advance. 

LVl.RViiquan)    returns  tlie  list  reader  qi-tan   and,  if  a  sublist 
is  being  scanned,  makes  the  reader  ascend  one  list  level  back  to  the 
point  of  descent. 

LVLRVl (quan)    returns  the  list  reader  quan   and,  if  a  sublist 
is  being  scanned,  makes  the  reader  ascend  as  many  list  levels  as 
necessary  to  return  to  the  point  of  descent  on  level  zero. 
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REEDiquan)    returns  the  list  element  that  the  quantity  quan 
currently  references  when  treated  either  as  a  list  reader  or  as  a 
sequence  reader.   The  element  returned  is  the  one  reached  on  the 
last  advance,  the  one  on  which  quan   sits. 
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The  function  LS'  >MP  described  in  §2.10  has  some  additional 
calls  which  are  useful  in  the  course  of  detailed  analysis  of  an 
expression. 

LSOCMP(ZtSi,  ltc,2,    3I!Tr.R)  returns  zero  if  tlie  terms  lisi 
and  lis2    (i.e.  sublists  of  expressions)  are  equal,  -1  if  lis\    is 
lexicographically  less  than  lis2,    +1  if  lis2    is  lexicographically 
less  than  lis\. 

LSQCMP(yari,  vaY'2 ,    3HFAC)  returns  zero  if  the  factors 
reached  on  the  right  advances  of  the  sequence  readers  in  the  variables 
or  array  elements  vav\    and  var2   are  equal,  -1  if  the  factor  reached 
in  right  advances  of  vavi    is  lexicographically  less  than  that  obtained 
from  var2 ,    +1  otlierwise.   The  next  right  advance  of  the  sequence 
readers  will  bring  them  to  the  powers  of  the  factors  which  are  not 
considered  in  the  comparison. 

LSQCMP(Zisi,  lis2,    3I1W0C)  acts  as  the  call  with  31ITER  in 
place  of  the  third  argument  3Htv'0C,  except  that  both  terms  are 
considered  to  liave  coefficient  1.,  i.e.  the  coefficients  are  ignored 
in  the  comparison. 
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3.9.   Attribute-Value  Lists 

Considerable  use  is  made  within  SYMBpLANG  of  a  SLIP  mechanism 
known  as  the  attribute-value  list,  which  allows  properties  of  lists 
to  be  declared.  The  attributes  used  internally  in  SYMBOLANG  are 
4HNAMR,  711N0RECUR,  6HS1MPLE,  and  5HVALUE.   These  should  be  considered 
reserved  to  the  package.  To  insure  compatibility  with  future  versions 
the  user  should  not  use  as  attributes  any  expression  symbol  of  four 
or  more  characters. 

Tlie  functions  ITSVAL,  MAKFDL,  NfTDLST,  NAMEDL,  and  KEWAL  provide 
access  to  attribute-value  lists. 

ITSVAL (auorz,  lis)   returns  the  value  of  the  attribute  puan 
for  the  list  or  expression  lis.      When  a  list  is  first  created,  all 
its  attributes  have  the  value  zero.   The  quantity  nuan   is  not 
restricted,  liut  use  of  list  in  this  context  is  not  recommended. 

NF.WVAL(^Ka^ii ,  quan2 ,    lis)    gives  the  attribute  quar.\    the 
value  aucm2   on  the  list  lis.      The  old  value  of  the  attribute  is 
returned  as  the  function  \'alue.   for  example  for  a  new  list  LA 

KA=NLWVAL(3HAGE, 36 , LA) 
KB=Ni;iVVAL  (SllAGL ,  37 ,  LA) 

will  leave  KA  set  to  zero,  KB  set  to  36,  and  the  attribute  3IIAnn 
on  the  list  LA  with  the  value  37.  Giving  an  attribute  the  value 
zero  is  as  good  as  removing  it. 
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Attributes  are  not  elements  of  the  list  with  which  they  are 
associated,  anil  cannot  be  detected  by  use  of  sequence  readers  or 
list  readers.   Rather  they  lie  on  a  separate  list  pointed  to  by 
an  otherwise  unused  field  in  the  head  of  the  list  with  which  they 
arc  associated.   The  exact  format  of  this  attribute  value  list 
should  not  be  assumed  by  the  jirogrammer,  so  that  it  may  !)e  freely 
varied  to  improve  the  efficiency  of  future  versions  of  SYMBOLAiNG. 

N;\MRni,(Zts)  returns  the  attribute-value  list  of  the 
list  or  expression  lis,    zero  if  there  is  none. 

MAKEDl,(Zisi  ,  lis2)    returns  the  list  or  expression  lis2, 
and  makes  the  list  lis\    the  attribute-value  list  of  lis2- 

Thus  I.B  may  be  made  a  very  complete  copy  of  list  or  expression 
LA  by  the  following  code. 

LB=LSSCPY(LA) 

LAV=NAMF.DL(LA) 

IF(LAV.NF,.0)  CALL  MAKF.I)L(LSSCPY(LAV)  ,LB) 

NfrDLST(Zis)  returns  the  list  or  expression  lis,    and  empties 
its  attribute-value  list.   Thus,  every  attribute  will  have  the  value 
zero  after  this  call. 
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3.10.   Another  Example 

As  a  final  example  of  SYMBOLANG  programming,  we  present  the 
code  of  a  preprocessor  for  SYMBCLANG  written  in  SYMB0LANG.   The 
bulk  of  the  code  is  in  the  subroutine  DKFS  which  defines  the 
expression  symbols  2HX.  and  2HF.  among  others.   Evaluation  of 
X.(A)  generates  code  for  the  representation  of  expression  A  as 
a  list,  while  evaluation  of  F. (A)  generates  code  to  combine  the 
expression  symbols  in  the  expression  A  using  LSQADD  for  addition, 
LSQMEX  for  multiplication  and  LSQRAZ  for  exponentiation. 

The  user  of  this  preprocessor  accesses  these  and  other  definitions 
by  presenting  FORTRAN-like  statements  which  begin  with  the  symbol  i 
(the  apostrophe  on  most  keypunches).   Lines  beginning  with  ^  in 
column  1  are  treated  as  comments.   Statements  which  do  not  begin 
with  i   are  transmitted  to  the  file  C0MPILE  unchanged.   Those  which 
do  are  evaluated,  and  the  generated  expressions,  if  any,  are 
transmitted  as  subroutine  calls.   INPUT  lines  are  sent  to  the  0UTPUT 
file  with  line  numbers  and  the  generated  code  with  the  leading 
cliaractcrs  I.SQ.   A  fragment  of  such  output  follows  the  program  listing. 

The  user  may  define  his  own  special  functions  by  using  the 
statement 

/PEF.  (symh {symb I , symb 2  , .  .  . )  ,fixpression) 
whicli  defines  sumb   with  dummy  arguments  symbi,   etc.,  to  be  expression. 


146 


3.63 


ICORE 


100 


ICORF 
ICORE 

42/0LlCORE,lfl/l 

0 

0 

=SSTAT 

CM, STAT. RECALL 

STAT 

30 

83,  XI 

ICORF 


PROGRAM  XLSO(  I  NPLIT  .  OUTPU  T  ,  COMP  1  LE  ,  TAPE  1  =  C0MP  I  LE  ,  TAPE5=  I  NP  t  T  , 
TAPE6=0UTPUT) 
CALL  BEGIN 
CALL  MIDDLE 
END 

rUMCTION  LSQERR( ! A) 
COMMON/TRAP/ITRAP.NOU,  IInHAT 
IWWAT=I  A 
GO  TO  ITRAP 
END 
IDENT 
ENTRV 
VED 
DATA 
MXA 
SA6 

MEMORY 

SAl 

SB3 

AX6 

EQ 

END 

SURROUTINfc  BEGIN 
C0MMON/TWAP/ITRAP,NOW,  I  WHAT 
COMMON  AVSL,X(100) . Y'50Un  ) 
ASSIGN  mo  Tn  ITRAP 
L=1C0RE(UUM)-HAD0V( Y)-3 
L=L/2*2 

CALL  INITASCV.L) 
CALL  DEFS 

CALL  LSUPNT(0,5HNOPRE) 
CALL  LS0PNT(fi,5HN0SUn) 
CALL  LSOPNT(0,5HFLOAT) 
CALL  LSQ0UT(5HCFLAG,(S,1H*) 
RETURN 
END 

SUBROUTINE  DEFS 
COMMON  AVSL»X(10n) 
COMMON/TRAP/ TTPAP, NOW, IWHAT 
COMMON/TEMP/LTEMP.L DP S.N TEMP, LINT, NIMT 
COMMON/ LSQNT R/M TRY PT( 10) 
COMMON/LSQARG/LARGdO) 
COMMON/LSQINT/INI 
DIMENSION  ICBU(10).NAM(1),ISBU(2) 

TRAPS  FOR  LSQPNT 

ASSIGN  1010  TO  LOC 

CALL  LSQPNT(  INl,4HFly.iL0C) 

ASSIGN  1020  TO  LOC 

CALL  LSUPNT( I N I , 2HS , , LOC ) 

ASSIGN  1030  TO  LOC 

CALL  LSQPNT( INI ,4HCAT.,L0C) 

DEFINITIONS 


xLSQfiaoa 

XLSOnooS 
XLSQn004 
XLSQ0005 
XLSO0006 
LSO600n2 
LSOE0003 
LSOE00n4 
LSOEOOOS 
LSOE0006 
ICOR00n2 
ICORP003 
ICOR0004 

icowocn5 

ICOR0006 
ICOW0007 

icnRoooe 

ICORP0n9 
ICOROOlO 
ICOROOll 
ICOR0012 
ICPR00t3 
BEGI0002 
B6GI0003 
REGI00CI4 

HEGinun5 

OEGiri006 
BEHinOO? 
RERI0008 
i^EGI00n9 
REGIOOIO 
BEGIOOII 
BEGI0012 
REGI0013 
BEG  I  0  01 4 
B6Gl00t5 
DEFS0Dn2 
nErS0003 
nEFS0004 
DErSOOOS 
DEFS0006 

nfcrsooo7 

DErSOOOS 
ntFS0009 
DEFSOOlO 

DErsnoii 
DEFS0012 

nEFsnoi3 
OErsooi4 
DEFSP015 
DEFS0016 
DEFS0017 

nEFSnoie 
nEES00i9 
DErS0020 
DErSOO?l 
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2000 


1010 

1011 
1020 

1021 

1022 

1023 
1030 

2010 


CONTINUE 
ASSIGN  2010 
CALL  LSODEF 
ASSIGN  2020 
CALL  LSQDET 
ASSIGN  3010 
ASSIGN  3011 
ASSIGN  3012 
CALL  LSQDET 

CALL  Lsonsr 

ASSIGN  3030 

CALL  LSQDEF 

CALL  LSODET 

ASSIGN  3040 

CALL  LSQDEF 

ASSIGN  3050 

CALL  LSQDEF 

ASSIGN  3060 

CALL  LSODEF 

RETURN 

SA"SEQL«(LA 

IF(N.NE,0)C 

IF(LSQTYP(S 

LC=COE 

ENCODE(20,1 

FORMaT( 120) 

CALL  LSQOUT 

CALL  TERM(0 

SAaSEQLRtLA 

IF(N,NE,0)C 

IF(LSQTYP(S 

DECODECIO,! 

FORMATClOAl 

LsO 

DO  1022  J=l 

IF(ICB'J(  J)  . 

CONTINUE 

IF(L,EO,0)C 

LL"»L*3 

ENC0DE(LL,1 

FORMAT( 12,1 

CALL  TERM(L 

CALL  STRDIR 

1F(N,NE,0)C 

CALL  NEWBOT 

CALL  VISIT( 

CALL  LSOUNM 

GO  TO  1030 

LA«LSQGAR(4 

LBaLSQGAR(0 

NLA»LSOTVP( 

IF(NLA,GF,2 

LTVsO 

1F(L8,EQ.Q, 

MLA»NLA-100 

LTVSLIST(9) 


TO  LOG 
CHDEF.'.LOC.l) 

TO  LOG 
(2HX. .L0C,-1) 

TO  LOG 

TO  IPX 

TO  IPT 
(2hf, ,L0C.-2) 
C3HF, .,LOr,-l) 

TO  LOG 
(2HT,,L0C.2) 
(3KT,  .',LOr,l) 

TO  LOG 
(3HRT.',L0C,2) 

TO  LOG 
(2HR, ,L0C.2) 

TO  LOG 
(2hI,,LnC.2) 

RG(6)  ,N) 

ALL  LS0ERP(4Hrix, ) 

A.GOE)  .ME.O  5CALU  LSQERR ( 4Hr I  X . ) 

Oil, irBU)LC 

(0,LANORM(IC8U(l)),0,LANORMnCBU(2))) 

.5 

RG(6) ,N) 

ALL  LSQERR(?HS, > 

A.NJAM)  .NE.lO000)CALL  LSQERR  (2HS,) 

021, NAM) ICrtU 

) 

,10 

NE.IH  )L=L*l 

ALL  LSaFRf'(2HS,  ) 

023,  ISBU)L,(  ICBUd  ),  1=1, L) 
HH,1GA1) 
SOODTfLL, IS^U)  ) 
(SE0LR(LAOG(6),N),LARG(l)) 
ALL  TERM(n.) 
(LARG(6),X(100)  > 
NTRYPT(3)  ) 
(X(IOO)  ,LARG(<S)  ) 

HDEF,  ) 

) 

L  A  ,  N  A  M  ) 

00C0.OR.NLA.LT,10000>GALL  LS0ERR(4Hr)EF,) 

OR.NLA .EQ.in000)RO  TO  2111 
00 


nEFS0022 
nErS0023 
r)Ersoo24 

nEFS0025 

nErsoo?6 

DEFSnO?7 
DEFS0098 
neFSP0?9 
TEFSno30 
DEFS0031 
nErsnu32 
nEFS0033 
nErS0U34 
DEFS0035 
DErS0036 
OEFSn037 
nErS0038 
DtrS0039 
nEFS0040 
DEfS0U41 
nEFS0042 
DEFS0U43 
DEFSn044 
DEFS0045 
nErS0046 
DEFS0047 
OEFS0048 
DErS0049 
DEFS0050 
riEFS0051 
nEFS0052 
nEFS0053 
DEFSn054 
DEFS0055 
DEFS0056 
DErsnos? 
DEFS00S8 
DErS0059 
DEFSn060 
DErSOUfSl 
DEFS00A2 
r>EFS00(S3 
nEFS0064 
nErS00/S5 
OEFS0066 
DEFSOOA? 
HEfSOUAS 
nEFS00(S9 
DEFSOOTD 
"ErSOOTl 
nEFSno72 
OEFS0073 

nErsoo74 
nErS0075 
nEFS0076 
nEFS0077 
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2112 
2111 


2020 


2021 


2023 
2025 
2026 


2024 


2027 


SA«SEQRDR(TOP(LA) ) 

SCOeSEQLR(SA,N) 

DO  2112  J=1,MLA 

KLABtSQTYP(SEOLR(SA,^J),NEM) 

IF(KLA,NF,10000)CALL  LS0ERR(4HDEr,) 

CALL  NEk^ROT(NEM,LTV) 

CALL  LSQDEF"(NAM,LB,LTV) 

NOWsl 

CALL  TERM(LIST(9)  ) 

CODE  FOR  X, (A) 

TO  BE  USED  WITHOUT  EVALUATIOM,  NO  SIDE 

RETURNS  THE  INTERNAL  REP  OF  A  EXCEPT  IF  A=USE,(B)  HR  asQ,( 
IN  THE  FIRST  CASE,  THE  VALUE  OF  B  IS  RETURNED,  IN  THR  SECC 
THE  INTERNAL  REPRESENTATION 

LAbLS0GAP(2HX, ) 
NLAaLSQTYP(LA,NAM) 
IF(NLA,NE,10001)SO  TO  2021 
SAaSEQRDR(TOP(LA) ) 
SC0bSE0LR(SA,N> 

IF(NAM,E0,4HUSE,  )CALL  TERM  (  V  I  S  I  T  (  LSOWL  (  SEQLR(  SA,  N  )  )  )  ) 

IF(NAM,NE,2HQ, )G0  to  2021 

CALL  STRDIR(SEQLR(SA.N),LA) 

SAsSEQRDR(LA) 

SATbSEOLR<SA,N> 

LUMBLSQMMLtLSQMNKl.n 

IF(N,EO,0)QO    to    2023 

CALL    MANV(T0P(LUM),LSQMNL(LSOMNLa.,LSOMNL(LSOMNL(9,  )).AhV 

1,)>,4HLIST,1,) 

CALL    TERM(LUM) 

S8«'SE0RDR(SAT) 

LUMTmi,SQMNL(LSf3MNL(l.#LSQMNL(LSQMNL(SE0LR(SB,N))))) 

SAXsSEQLR(SH,N) 

IF<N)2024, 2026, 2027 

LARsLS0HNL(LSaMNL(l,,SAX.2HX,,l,)) 

CALL  MANY(X(100),LUM,LUMT,LNKR(SA>,LNKR(Se>,LAR) 

CALL  STRDIR(  VIS  IT  (LSOWL  (LAR)),LBR) 

CALL  LSQUNM(X(ioo).LUM,LUMT,SA,SB,LAR) 

CALL  LSODES(LAR) 

CALL  NEWe0T(LBR,30T(LUMT) ) 

GO  TO  2025 

CALL  NEWROT(LSOMNL(LSQMNL(1,,LSOMNL(LSQMNL(1.,SAX,1.) >,2KS 

>),BOT(LUMT)) 

SAXsiSEQLR<SB,N) 

IF(N,EQ,0)GO    TO    2026 

CALL  NEW80T(LS0MNL(LSQMNL(SAX) >,90T(LIMT)  ) 

GO  TO  2025 

CALL  MANY<80T(LUMT),6HLSOMNL,1.  ) 

CALL  NEWR0T(LUMT.BOT(LUM) ) 

SAT8SEQLR(SA.N> 

1F(N.EO,0>GO  TO  2023 

CALL  MANY(BOT(LUM) ,6HLSQMNL,1,  ) 

CALL  TERM(LUM) 


R) 


F, 
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DEFS0078 
DEFS0079 
DEFSOOaO 
DEFS0081 
nErS0082 
DEFSn093 
DeFS00fi4 

DErsnoss 

DEFSP0fl6 
DEFS0087 

nEFsoose 

DEFS0O«9 
DEFSOOPO 

ND.  nEFS0091 
DErS0092 
nEFS0093 
DErS0094 
DEFS0095 
DErS0096 
nErS0097 
DEFS0098 
DEFS0099 
DEFSOinO 
DEFSOlOl 
DEFSnin2 
DEFS0103 
DEFS0104 
DEFS0105 

lx.,nErsoio6 
nEFSP107 
DEFS0108 
DEFSPin9 
DEFSOllO 
DEFSOlll 
nEFS0112 
DEFS0113 
DEFS0114 
DEFS0115 
DEFS0116 
DEFS0117 
DEFSOlie 
OEFS0119 

■  ,1,DEFS0120 
DEFS0121 
DErS0122 
PEFSP193 
DEFS0124 
DEFS0125 
DErS0126 
DEFS01?7 
DEFS0128 
DEFS0129 
OEFS0130 
DEFS0131 

nErsoi32 

0EFSni33 


3.66 


3nio 

3011 
3013 


3014 


3024 


3015 


3012 


LA=LSQGAW(?Hr, ) 

CALL  TERM(VISIT( IPX)) 

PROCESS  r,  (EXPRESSION) 


SAsSEQ 
MT  =  0 
SAR5SE 
I F  (  N  ,  N 
CALL  M 
VsVISI 
CALL  L 
CALL  N 
NTaNT* 
GO  TO 
!r(NT, 
LSOMML 
LNAM«6 
IF(NT, 
NTaNT- 
LAN=LS 
LUN=LS 
CALL  M 
V  =  VISI 
CALL  L 
CALL  L 
CALL  N 
CALL  L 
NT  =  NT- 
LUM8LS 
CALL  M 
V=VISI 
CALL  L 
CALL  L 
CALL  M 
CALL  L 
IF(NT. 
CALL  M 
))> 

VsVJSI 
CALL  L 
LANsLS 
CALL  L 
GO  TO 


RDR(LA) 


QLR 
E.O 
ANY 
T(  1 

SQU 

EWB 

1 

301 

EQ. 

(LS 

HLS 

EQ, 

1 

QM^' 

QMN 

ANIV 

T(L 

SQU 
SOD 

Ewe 

SQO 

1 

QMN 

ANY 

T(L 

SQU 

SOD 

ANY 

son 

EQ. 
ANY 


(SA.iJ) 

)G0  Tn  3014 

(X<10  0).NT,LNKR(SA)) 

PT) 

NM(y(100) ,NT,SA) 

OT(V,X(100) ) 

3 

0)CALL  TERM(LSQMNL(LS0M^IL(1.,LS0^^^^L(LS0^'NL(1 

QMNL(9.)).4Mf-lx.,l,)),4WLlST,l,))) 

QAUD 

1)  CALL  TERM(P0PBOT(X(100) ) ) 


Ld.) 

L(LSOMNL(1..30T<X(inO) ),3HT, 

(XdOO)  .LNAH.NT.LUN.LAN) 

SOVVLaUN)  ) 

NM(X(lon),LMAM,NT,LUN,LAN) 

ES(POPHOT(XflOO)  )  ) 

0  T  (  V  ,  L  A  M  ) 

ES(LUM) 


.1.  )) 


)  ) 


L(LSQ^NL(l..BnT(X(10  0)),3HT,,,l 

(X(lOn)  .LNA'SNT.LUN.LAN) 

SOVVL<LUN) ) 

NM(X(100),LNAM,NT,LUN,LAN) 

ES(POPBOT(X(100) ) ) 

(LAN, V.LNAM.l, ) 

ES(LUN) 

0)CALL  TERM(LSQMNL(LAN) ) 

(X(100),LNAH,NT,LSQMNL(LSOMNL(l..LSCMML(L'iN),:^HT,.,l 


T(LSQVVL(8OT(X(100) ) ) ) 

SQUNM(X(100).LNAM,NT,LUN) 

0MNL(1.,V) 

SQOES(LUN) 

3015 


PROCESS  F,(TFRM) 

SAR  =  SEQRiiR(SAR) 

SATC»SEOLR(SA«,N) 

CALL  NEWR0T(LS0MNL(LSQMNL(1,  ,LSCMNL<LSnMNL(l.,LSOM^L(LS0f'NL( 

SATC)).6HLSQMNL,l.)).6MLSQMNL,l,)),X(in0)) 

NTal 

IF(SATC,\'E,l)Gn  TO  3016 

SARAcSAR 

SARR=SEQLR(SARA,N) 


nEFSni34 
nErsni-^5 
nErsoi36 
nEFsni:^? 

rjErsni39 
nErsPl40 
DErsni41 
ntrsni42 
nEPSni43 
nErsni44 
nErsni45 
nErsni46 

DEFSni47 

nErsoi48 

nEFSPl49 

nErsni50 
nErsoi5i 
nEFSni52 
nErSPl53 
nErsPi54 
nErsni55 
nErsniS6 
nErS0157 
nErSni58 
nErsniS9 
nEFSPl(«,0 
lEFSPlftl 
nEFSniis2 
nErS0163 
')EFsni'>4 
nErsni65 
DErsnift6 
nErsoi67 
nEfSnice 
nbFsni69 
nErsni70 
nErsni71 
,nEFSni72 
nEFS0173 
Def Sni74 
nErsni75 
nErsni76 
nEFSni77 
nErsrii78 
DEf Sni79 
^ErSClwO 
IiEFSniRl 
nErsnifl2 
nErsnia3 
PErsni'^4 
HErSf^lf^^ 
nEFsni86 
nErSPlfl7 
PEf SPiftS 

nEFSPlR9 
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ir(N,GT,G)GO  TO  3017  nErsni<30 

CALL  LSQD6S(POPHOT<X(1(10)  )  )  nfcrS^l'51 

NTbO  Dersri92 

3016  SARAaSAR  D6rsni03 
IFfLSQGNr (SAP,NAM,NAR> .GT,0)GO  TO  3017  ncrsniq4 
NT»NT*1  nErSni05 
ir{NAM,E'3,3Hl,»-)R0  TO  3019  nerSni<)6 
LUN«lSQMNL(1,  )  ')ErS0197 

3018  CALL  NEWfJOT(SEQLR(SARA,N),LUN)  flErSMoe 
ir(N,EQ,0)GO  TO  3018  nErSfH09 
CALL    NEW30T(LS0MNL(MANY<LUN,1  ,  )  ),X(irj0)  )  D6rS02no 

3022  SPOW=SEQLR(SAR.N)  nErSn2ni 
lF(tM.GE,n)GO  TO  3020  nEfSn2n2 
IFtSPOW.ME.DGO  TO  3021  nErSn2n3 
GO  TO  3016  OErSn204 

3019  IFdvlAR.ivlE.DCALL  LSQPRR(2Hr,)  nErSn2n5 
LUNcLSQMvJLCLSQMNLd.  .SFQLR(SARA,M),3Hr.  ,  ,1.  )  )  nErsn2n6 
CALL  MANY(X(10n),LNKR(SAR),NT,HJN)  DErSn2n7 
VsVISIT(LSOVVL(LUN)  )  DErsn2n8 
CALL  LSQUNM(X<ton).SAR,NT,LUN)  nErSn2n9 
CALL  NEWHOT(V,X(100)  i  DErsC2l0 
CALL  LSODES(LUM)  n£rsn2ll 
GO    TO    3022  nErsn2l2 

30  21      LUN8LS0;1NL(LSQHML(1,.LS0MNL(LS0KNL(1  .  ,LS0MNL(LSQHNt(SPOU|)  >,    n6FS0  2t3 

*  6HLSQMNL1I,  )  )  .6HLSQMNL.1.  )  )  nErsn2l.4 

3023  LRUN=LSQMNL(LSr)MNL(l..LSGMNL(LSCM\L(l,.POPROT(X(10n)),3HT...   nErSn2l5 

♦  l,)),LSQMNL(LS0MNL(l.,LLM.3HT,,,l,)),6ML?0PA7,i,))  nErSn2l6 
CALL  MAi\IY(X(10n)  ,LNKR(SAR)  .NT.LRUN)  DETSn?!? 
V5VISIT(LS0VVL(LRUN) )  DErSn2ie 
CALL  LSQUNM(X(10n) .SAR.NT.LRUN)  D6rSn2l9 
CALL  NEWRQTC V,y(100) )  nErsn220 
CALL  LSQDES(LHUN)  nErS02?l 
GO  TO  3016  nErsn2?2 

3020  LUNeLSOMNL(LSQMNL(l. .SP0W,3Hr, , ,1, n  0ErSn2?3 
CALL  MAMY(X(1U0) ,LNKR(SAR).NT,LLN)  nErSn2?4 
VsVISIT(LSQVVL<LiJN)  )  DErsn225 
CALL  LSQUNM(X(10n).SAR,NT,LUNA)  0Ersn226 
LUN=LSQCPY(SFQRDR(V) )  DErSn227 
CALL  LSO[)ES(LUNA)  n£rsn2?e 
CALL  LSQDESCV)  nPrsn2?9 
GO  TO  3023  nEFSf)230 

3017  LNAM=6MLSQMEX  n£rSn231 
GO  TO  5024  OErS0232 

C  PErsn233 

C         T,  nErSfl234 

C  ntrsn235 

3030  LA  =  LSQGAR(2HT,  )  aErS'1236 
SAsSEQRDRdTEMP)  n6rsr237 

3032      SAN=SEQLR(SA,N)  nEFSn23« 

1F(N,GT,0)GO  TO  3031  nErS0239 

IF(LSQCMP(SEQLR(SA.N),LA,3HEXP1 ,ME.0)GO  TO  3132  nEFS0240 

CALL  TERM(LSOMN)LaSQHNL(l,.SAN,l,  )  )  )  nEFS0241 

3031  NLA=LSOTYP(LA,MAM)  nEFS0242 
IF(NLA,NE, 10000)50  TO  3033  DEFSn243 
SA«SEQRDR(LDES)  nErbn244 

3037      SAN»SEOLR<SA,N)  nErS0245 
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3036 
3033 
3034 


3040 
3047 


3045 


3046 


3044 

3043 
3042 


3041 


3050 


3060 


3062 


ir(N,GT,0)GO  TO  3036 

ir(LSOCNM(SAN,\iAM).En.O)CALL  TERM(LA) 

GO  TO  3037 

CALL  NEWF<OT(NAM,LDFS) 

CALL  TERM(LA) 

NTEMP=NTEMP*1 

ENCODE (6. 30  34, NAM )NTFHP 

rORMAT(3HL?0,03) 

CALL  MANJY(LTEMP,NAM,LA) 

LA  =  LSQMNL(LSOMML(l.#NAM,l,  >  ) 

GO  TO  3036 

RT. 

IRHaO 

LA"LS0GA«(3HRT.  ) 

NLA»LSOTYP(LA,MAM) 

IF(,NLA,NE,10000)GO  TO  3043 

SAaSEQ'^UR(LDES) 

MDES=LIST(9) 

SAMsSEOLUSA.N) 

irdM.GT.DGO    TO    3044 

ir(LSQCN'1(rgAM,9AM)  .En.O)    GO    TO    3046 

CALL    NE>^BOT(SAM,MDES) 

GO    TO    3045 

NTEMP=NTEMP*l 

ENC3DE(6,30  34,NAM)ixlTEMP 

CALL    MANY<LTEMP,KiAM,I.A5 

CALL    NEWyOT(NA*^,MDES) 

GO    TO    3045 

CALL    LSOI)ES(LDES> 

LDEScMDES 

SAaSEQRL)M(LTEMP) 

SANJ»SEQLW(SA,N) 

IF(N,GT,0)GO  TO  3041 

IF(lSQCMP(SEQLR*(SA.N)  ,LA.3MEXP)  ,NE.0)Gn  TO  3042 

CALL  TERM(LA) 

IF( IRM,NE,0)CAll  TF«m(lA) 

NTEMP=NTEMP*1 

ENC0DE(6,3n34,NJAM)  NTE^P 

CALL  MANY(LT5MP,NAM,LA) 

GO  TO  3036 

R. 

IRMal 

LA  =  LSQGAR(2HR,  ) 

GO  TO  3047 

I. 

LA  =  LSQf^AR(?HI  .  ) 

NLASLSQTYP(LA, TYP) 

IF(NLA,tO,C)CALL  TFRM(LS(3MNL(LSQMNl(1,.LA,4HFIX,,1.))) 

SA=SF0RDR(LIMT) 

SAM=SEQLR(SA,N) 


nEFSn246 
nfcrsn?47 
nErsn248 
nEFSn24y 
ntf"S0250 
nErS025l 
nEFS0252 
nEFSn2«;3 
nEFSn254 
HErSO^^S 
nEFSn2^6 
nEFSn2S7 
DEFSP258 
nersn259 
n6rS0260 
DEFSn261 
n£FS02^2 
nErsn2iS3 
nEFSn264 
nEFSn?65 

nEFSfl266 
nEf S02A7 

nEFsn26e 
nEFsn269 

OEFSn270 

nersn27l 
nEFsn272 
DEFSn273 
nEFS0274 
riEFS0275 
nEFS0276 
nEFSP277 
nErsn278 
DErSP279 

nEFsn2«u 
nErS0281 
nEFSn2«^2 
n&rsn2fi3 
nEFsn?ft4 

nEFS'^2^5 
nEFSn2«6 
nErsn2R7 
nEFSOi''<e 

nErS02q9 
nEt"S02O0 
DErsC291 
0EFSn292 
nerS0293 

nEFSn294 
nErsn295 
nErsn296 
nfcFsn297 
nErsn298 
nfcFSri299 
nEFSn3no 
nEFSnini 
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5.69 


3061 


1 

100 


21 

210 


6 
5 

50 
4 


41 


IF(N,GT,0)GO  TO  3061 

ir(LSQCMP(SEOLf (SA.N),LA,3H6XP) ,ME.0)GO  TO  306? 

CALU  TE'^M(LS0MML(LSQ'^''JL(1..SAN,1,))) 

NTEMP  =  .MTEMP*1 

ENCODE (6i 3034, NAM )NTFMP 

CALL  MAMY(LIMT.NAM,LA) 

CALL  TEKM(LSOMML(LSOMNL(1  ,,NAM,i, ) ) ) 

END 

SUBROUTINE  MIDDLE 

DIMENSION  lNBUr(O0),TSM(5).0SB(8),M(l) 

COMMON/TRAP/ ITRAP, NOW, IWHAT 

COMMON/TEMP/LTEMP.LDFS.NTEMP.LINTiNJNT 

COMMON/PAGE/LIME,  IPAr5E 

DATA  lERH/n/, inrr/0/ 

ASSIGN  999  TO  ITRAP 

READ(5,100) INBLT 

FORNAT(90R1 ) 

IF(ENDFILE  5)2.3 

ir(  lEHR.NE.O)  GO  TO  21 

CALL  EXIT 

ENC0DE(60,2in,OSR) lERR 

F0RMAT(4H  *i»*,T6,*  L"?!;!  ERRORS*) 

CALL  AgORT(OSB.6L*nXlT. ) 

CALL  LPRINT(  INRUF.fi) 

IF(  lOFF.NE.Q)  f,0  TO  5 

IF(  INBIJF(I)  .NE.IR?!)  r.Q    TO  4 

ENCODEdO.lDO.H)  nNBlJF<  I  )  ,  Is?»ll) 

IF(M,EQ,4HLAST)     GO     Tn    7 

1F(M,E0,3H0FF)     GO    TO    6 

GO    TO    1 

10FF  =  1 

GO    TO    i 

IF(  INBUF(I)  ,NE.lR;«f ) 

ENCODE(10,100,M) ( IN 

!F(M,NE,2H0N)    GO    TO 

IOFFbQ 

GO    TO    1 

WRITE(1,100)     INBUF 

GO    TO    1 

CALL    SETRAY( !SR,5,lR 

IF( < IN8UF(l) .NF.IR     ) 

GO    TO    ^0 

1F( INBUF(6) .NE.IR     )     GO    T"    50 

M  =  5 

DO    41    J=l,5 

LU»6-J 

1F<INBIJF(LU>.E0.1R 


r,0    TO    50 
IjFd  ),  Is2, 
50 


) 


11) 


DEI 
nE( 
DE( 

HE! 

nfcf 

DEt 
OEf 


AND,  (  INBUF  (1  ).LT.1.R0,OR.  INBUF(  D  ,GT 


1R9) )m: 

M] 
M  ] 

M : 

M  1 


)    GO    TO    Al 

!F(  INBUF  (LU)  ,LT,1RC,0R.  INBUF(Ll))  ,GT.1R9) 

1  SB<M)s INBUF (LU) 

M  =  M-1 

CONTINUE 

1STATN=10H 

IF(M.LT,5)ENCOnE(5.100,ISTATN)!SB 

!SEND=0 

DO    77    J=7,72 

IF(  lN8UF( J) ,EQ,1R     )     GO    TO    77 


GOTO    50 
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1F(  IMBOF  (  J)  .NE.lR^)  r,o  TO  50 

M  =  J 

GO  TO  78 

77  CONTINUE 
GO  TO  50 

78  INRUr(M)=lP 
LIS=LIST(9> 

79  ENCODE( 7 2,101.0 SB) (IM8UF( I ), 1=7,72) 
101       rORMAT(6X,72R11 

LYS=LIST(9) 
DO  80  J=1,R 

80  CALL  NFWR0T(0Sf1(  J),LyS) 
CALL  NEWBOT(LYS,LIS) 
INBUr(l)=lRC 
IN«Ur(2)=lR 
INBUF(3)=1RL 
1NBUF(4)=1RS 
INRUF(5)=1RQ 
WRITE(l.lOO) INRUF 
READ(5,1U0>  iNBijF 
IF(ENDFILE  5)61., fl2 

81  1SEND=1 
GO  TO  93 

82  IF( IMBUF(6) .EO.IR  )Gn  TO  83 
IF( INBUF(I) ,ME,1R  )G0  TO  83 
CALL  LRRINT(  It>|F^UF,0) 

GO  TO  79 

83  NOWso 
NTEMPsO 
LTFMP=LIST(9) 
LDFS=LIST(9) 
NINT=U 

LINT=LIST(9) 
LAM=INLIST(LAM,LrS,0) 
CALL  LSQUES(LIS) 

912       IF(LISTNT(LINT) .FQ.0)GO  TO  911 

SAN=POPTOP(LINT) 

LYS=LIST(9) 

CALL  L5QOUT(4HLI';T,LVS.6HMARGIN,1,0,  I  STATN  ,  6HMARr,  I  v,  7  ,  n  ,  c  AM  ) 

CALL  LS00UT(3,3H  =  ) 

ISTATN=10H 

SANsPOPTnP(L!NT) 

CALL  LSQt'NTCSAN,!  ) 

CALL    LS00UT(5MFLi)SH) 

CALL    LSQl)ES(SAN) 

CALL    LOoT(LYS) 

GO  TO  912 
911       IF(LISTMT(LTFMP) .EO.0)PO  TO  910 

SAN=POPTOP(LTEMP) 

LYS=LIST(9) 

CALL    LSQOUT(4hLlST,LVS,6MMARGIN,l,0,ISTATN,6HMARr,  I\,7,12, 
*  12HCALL    STRDIRO 

ISTATN=10H 

SAMaPOPTOP(uTEMP) 

CALL  LSQPNT(SAM,l) 

CALL  LSQOOT(1,1H,,0,«?AN,1,1H),5HFLUSN> 


15^ 
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CALL  LOUT(LYS) 
GO  TO  911 

910     iF(Now.Eu. n)no  to  920 

CALL  LSQDE^(LAM) 

930  IF(L1STMT(LDES) .EQ.OJGO  TO  931 
LYS=LIST(9) 

CALL  LSOOUT(4HLIST.LYS.6MMARF;IN,1,0,ISTATN.6HMaRGI\,7,12, 
•     12HCALL  LSODESf  ,0,POP8nTaDE??).l,lW)  ,5HFlUSH) 
CALL  LOUT<LYS) 
GO  TO  930 

931  CALL  LSQDESCLTEMP) 
CALL  LSQDES(LINT) 
CALL  LSQDES(LDFS) 

903       IF(  ISE'^iD.NE.  0)7,3 
920       ASSIGN  930  TO  LAC 

IF(LSOTYP(LAM,NJAM)-200nO)9?2,923,9?3 

922  IF(.NAM.NE.3HD0,  )G0  TO  023 
ASSIGN  924  TO  LAC 
LUMSLAH 

SAT=TOP(LAM) 
SAR=SEQRDR(SAT) 
SAC  =  SEQLR(:^AR,N) 

924  CALL  STRDI»(SE0LR(?;AP,M),LA.'1» 
IF(M.EQ,0)fiO  TO  923 

CALL  LSODES(LUH) 
GO  TO  930 

923  LYS=LIST{9) 

CALL  LSQOUT(4HL1ST.LYS.6MMARGIN,1,0,ISTATN,6MMARGP',7) 

ISTATNslOH 

IF(LSQTYP(LAM,MAM)-2nonO)925,9?6,9?6 

925  IF(NAM,NE.4HCAT. )G0  TO  926 
927       CALL  LSOPNT(LAM,l) 

CALL  LS00UT(5HFLUSM) 
CALL  LSODES(LAM) 
CALL  LOUT(LYS) 
GO  TO  LAC, (930.924) 

926  CALL  LSQOUT(5,'5HCALL  ) 
GO  TO  927 

999       IERR=IERR*1 

ISP(1)=10H  ♦•*  ERROR 

ISP(2)=10H*** 

ISB(3)=IWHAT 

CALL  LPRINT( ISR,3) 

1F(  IWHAT.NF.6HNUCELL)  r.Q    TO  903 

CALL  REMARK(28L  ***  MEMORY  OVERFLOW  ♦*• 

GO  TO  21 

END 

FUNCTION  LPR1NT{L8UF, IC) 

DI^^ENSION  LBUF(90) 

COMMON/PAGE/LINE, IPAGE 

DATA  LINE/55/, IPAGE/0/,LPP/55/ 

DATA  ISLIN/0/ 

LlNEsLlNEfl 

IFCLINE.LE.LPP)  GO  TO  200 

IPAGF=IPAGE-^1 

CALL  DATE(XD) 


m: 
m: 
M : 

Ml 

LPRIi 

LPPI 

LPPI 

LPPI 

LPPI 

LPRI 

LPRII 

LPRIl 

LPRI 
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3.72 
CALL  HOUR(XT) 
WRITF(6,101)XT.XD, IPAGF 

rOWMAT(*l  LSO  1.0  SYMBOLANG  PRFPPOCESSAR 

■*.A10.*  DATF  -*.A10i*  PAGF  -*.I6,/) 

LINE=3 

IFMC.ME.O)    GO    TO    30  0 
ISLIM=ISLIN*1 
WRITE(6,201)  ISLIN.LSlJF 
FORMATdX,  I8.1X,90R1> 
RETURN 

WRITE(6,301)(LPUr(J),Jal,IC) 
F0RMAT(6X,«LSQ    *.8A10) 
RETURN 
END 

FUNCTION  LOUT(LYS) 
DIMENSION  I0BU(8) 
IF(LISTMT(LYS) .EO.0)GO  TO  1 
SAM=POPTOP<LVS) 
L  =  l 

CALL  SETRAY( 10RU,8.lnH  ) 

1F(LISTMT(SAM) .EO.0)GO  TO  3 
CALL  STRDIR(P0PT0P(SAM),  lOBUa)  ) 
L  =  L*1 
GO  TO  2 

CALL  LPRINT( I0»U,8) 
WRITE(1,4)I0BU 
F0RMAT(8Atn  ) 
CALL  LSQDES(SAM) 
GO  TO  5 

CALL  LSQDES(LY5;) 
RETURN 
END 


LPRIOOll 

LPRinOl.2 

TIMF  LPRI0013 

LPRinoi4 

LPPinOl5 

LPRIPU16 

LPRinoi? 

LPRinois 

LPRinOl9 

LPRin020 

LPRI0021 

LPRI0022 

LPPlnO?3 

LPR10024 

LOUT0002 

LOUT00n3 

LOUT0004 

LOUTn005 

LOUTn006 

LOUT0007 

LOUTH  0  08 

LOUT0009 

LOUTOOIC 

LOUTOOtl 

LOUT0012 

LOUT0013 

L  0  U  T  n  0 1 4 

LOUT0015 

LOUTnOt6 

LOUT0017 

LOUT0018 

LOUTP019 
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3.73 
LSQ    1,0    SYMBOLAMG    PREPRC5CESS0R 


TIME  •  10,35.51 


DATE 


02 


SOURCE  CARDS  ARE  SCANNED  FOR  AN  t    ON  THE  FIRST  CARD  OF  A 
STATEMENT.   IF  NOT  FOUND  THE  STATEMENT  IS  TRANSFEREO  TO  THE 
COMPILE  FILE. 

IF  AN  t    IS  FOUND  RETWEEN  COLUMNS  7  AND  MI  THE  SCURCE 
STATEMENT  IS  OUTPUT  AS  A  COMMENT  AND  THE  EVALUATE  OF  THE 
SOURCE  STATEMENT  IS  OUTPUT  AFTER   ^CALL*. 

THE  NEXT  CARD  DEFINES  THE  EFFECT  OF  = 
^DEF.(A...lsA...2iSTRDIR(A...?,R,(A...l>)) 

NOW  FOR  A  PROGRAM 

PROGRAM  TESTCINPUT, OUTPUT) 

DIMENSION  SPACE(5000) 

CALL  INITAS(SPACE.5000) 

ADD  A  COLUMN  OF  EXPRESSIONS 

^LA=X. (0) 
CALL  STRDIR(LIST(9).LA) 

LB=INLIST(LR,5L1NPUT) 

IF(LISTMTaR)  .EO.0)GO    TO    2 

^LSQPNT(LR.S. (LB) ) 
CALL    LSOP\IT(L0.    2HLB) 

>fLA  =  F.  (LA*LB) 
CALL  STRDIR(LA,LS0001) 
CALL  STRDIR(LSQADD(LB,LA),LA) 
CALL  LSODES(LSQOOl) 
CALL  LSODES(LR) 

GO  TO  1 

^LSQPNTfLAjS. (LA) ) 
CALL  LSQPNTaA,  2HLA) 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


OW  LETS 
LAaF.(L 
STRDIR( 
STRD!R( 
STRDIR( 
STRDiRt 
STRDIR( 
STRDIR( 
STRD!R( 
STRDIR( 
STRD!R( 
STRDIR( 
LSODES( 
LSODFS( 
LSODFS( 
LSODES( 
LSODESt 
LSODES< 
LSODES( 
LSODES( 
LSQDES( 


TAKE  T 
A*LA**2 
LSQMNL( 
LSQMNLC 
LSQMNL( 
LSORAZC 
LSQRA7( 
LSOADn( 
LSQRA7( 
LSQADDC 
LA.LSOO 
LSQADO( 
LSQOIO) 
LSOn07) 
LSQOOA) 
LSQ005) 
LSQOOO 
LSQ003) 
LSOOO?) 
LSQftOl) 
LSQOll) 


HE  SUM  OF  SOME  POWERS 
♦LA**3*LA**4) 
LS0MNL(2. )),LSO001) 
LS0MNL(3. )),LSO002) 
LS0MNL(4. ) ) ,LSO0n3) 
LA,LSO003) ,LSC004) 
LA,LSO002) ,LSC005) 
LSO004,LSO005>.LSQ006) 
LA.LSOOOl) ,LSG007) 
LSO006,LSQ007), LSQOIO) 
11) 
LSQOIO, LA) ,LA) 
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5.1 


INDEX 


3H1. 


see  parentheses 


}    (advance  list  reader)    3.51,  3.56 


accuracy  of  numbers  in  expressions   2.7 

ADD   (add  expressions  -  obsolete)    3.29 

AUVLKL   "^ 

ADVLER 

ADVLNL 

ADVLNR 

AUVLWL 

ADVLWR 

ADVSEL 

ADVSER 

ADVSNL 

ADVSNR 

ADVSWI, 

ADVSWR 

AI,0AD  (for  use  in  argument  list  scanning)   1.24 

analysis  of  expressions   2.35-2.40,  3.51-3.59 

arithmetic  on  expressions   2.21-2.22,  3.29,  3.62 

ASSIGN  statement    1.20 

attribute-value  lists   3.60-3.61 


B0T   (bottom  of  list)    3.51 

BREAK   (decompose  expression  -  obsolete)    3.51,  3.52 


Chippewa  Operating  System  1.3 
code  generation,  lORTRAN  3.18 
C0MM0N   1.16 

/ESQARC/ 

/LSOCSX/ 

/LSQDMF/ 

/ESQ I NT/ 

/LSQNTR/ 

/SYMSYS/ 


3.21,  3.62 


3.21 

3.6-3.7,  3.26-3.27 

3.37-3.38,  3.39 

3.21 

3.21 

1.16 


i6o 


5.2 


constraints  on  SYMBOLANG  programs   1.15-1.16,  2.42 

r0NT   (contents  of  machine  address)    1.25 

Control  Data  6600   1.1 

core  requirements   see  constraints 

CPYTRM   (copy  term  -  obsolete)    3.1-.3.2,  3.26 

creation  of  expressions   2.1-2.5,  3.1-.3.11 

definition   see  evaluation 

OHLETE   (delete  list  element)    3.26,  3.27 

destruction  of  expressions   2.17,  3.26-3.28 

determinant   2.54-2.58 

differentiation   2.31-2.32,  2.45-2.46 

DVSIJM   (divide  expressions  -  obsolete)   3.29 

entry  points  in  SYMBOLANG   1.17-1.18 

EQUAL   (compare  quantities)    1.25,  2.35,  2.40 

6I1EQUAL.    3.16,  3.40 

evaluation  of  expressions   2.24-2.26,  2.43-2.44,  3.30-3.45 

EVALUE   (evaluate  expression  -  obsolete)    3.30,  3.44 

expression   1.5-1.10 

conditional  1.22,  2.44,  2.48 

constant  1.6,  2.49 

dumps  of 2.19 

syntax  of  (input)   2.5-2.6,  3.14 

syntax  of  (internal)   1.5 

factor   1.5,  3.59 

GETC0E   (get  coefficient  -  obsolete)    3.51,  3.52 

head  of  a  list   1.4 

IIITENT   (substitution  in  expression  -  obsolete)    3.30,  3.44 

IBM  7094    1.3 

IFTHEN   (FORTRAN  callable  conditional  expression)    1.22-1.23 

3HIF.    2.44,  2.48 
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5.3 


INHALT  (contents  of  machine  address)    1.25 

INITAS  (initialization  routine)    1.16,  2.24,  2.42 

initialization   1.16,  2.42-2.46 

INITRD  (reinitialize  list  reader)   3.51,  3.56 

INLIST  (input  expression)    1.10,  2.5-2.11,  3.4,  3.14-3.17,  3.26,  3.40 

INLSTL  (insert  into  list)    3.1,  3.2-3.3,  3.4,  3.26 

INLSTR  (insert  into  list)    3.1,  3.3-3.4,  3.26 

input  of  expressions   see  INLIST 

I.NSBST  (substitute  in  expression  -  obsolete)   3.30,  3.44 

INSUBT  (substitute  in  expression  -  obsolete)   3.30,  3.44-3.45 

INTARG  (for  use  in  argument  list  scanning)   1.24 

INTE.NT  (substitute  in  expression  -  obsolete)    3.50,  3.45 

INTGER  (provide  integer  name)   1.19,  1.25 

IRALST  (erase  list  -  use  with  caution)    3.4,  3.5,  3,26,  3.27 

IRARDR  (erase  list  reader)   3.51,  3.56 

ITSVAL  (obtain  value  of  attribute)    3.60 

job  structure   1.15 

LCNTR  (level  counter)   3.51,  3.57 

LUIF   (differentiate  expression)   2.31-2.32,  2.45 

level 

definition  3.32 

expression  3.30 

list  3.55 

list   1.3-1.5 

constraints  on  the  elements  of  a  1.5 

LIST   (create  list)    2.1,  3.1,  3.4-3.5,  3.26,  3.27 

LISTNfT   (test  for  empty  list)   2.35,  2.39 

LIST0N   (input  expression)    3.17 

L0CARG   (for  use  in  argument  scanning)   1.24 

L0CT   (test  list  for  validity)    3.51 

L0FRnR   (list  being  scanned  by  list  reader)    3.51,  3.57 

L00KUP   (define  derivative)    2.31,  2.32,  2.42,  3.4,  3.26 
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I,PNTK   (list  pointer  of  list  reader)    3.51,  3.57 

LRDRCP  (copy  list  reader)    3.51,  3.57 

LRI)R0V  (establish  list  reader)    3.51,  3.56 

LSQAUU  (add  expressions)    2.3,  2.21-2.22,  3.7,  3.26 

LSQCMP  (compare  expressions)    2.35,  2.39,  3.59 

LSOCNM  (compare  expression  symbols)    2.35,  2.39-2.40 

LSQCPY  (copy  expression)    3.1,  3.5-3.6,  3.26 

LSQCXP  (copy  expression)    3.1,  3.6-3.7,  3.27 

LSQDIiF  (define  value  of  expression  symbol)    2.24,  2.25-2.26,  2.42, 

3.4,  3.16-3.17,  3.26,  3.30,  3.32,  3.33-3.38,  3.43 

LSQDES  (erase  expression)    2.17,  3.4,  3.5,  3.26,  3.38 

LSODSF  (simple  function  definition)    3.30,  3.38-3.39 

LSQGAR  (get  argument)    3.30,  3.38,  3.39 

LSQIUR  (initial  derivative  definitions)    2.42,  2.45-2.46 

LSQINI  (initial  evaluation  definitions)   2.24,  2.42,  2.43-2.44, 

2.46,  2.47 

LSQM1-;X  (multiply  expressions)   2.21-2.22,  3.7,  3.26 

LSQMNL  (create  list)    2.1,  3.26,  3.51 

LSQ0UT  (general  output  routine)    2.13,  3.19-3.22 

LSQPNT  (output  expression)    2.13-2.15,  2.19,  3.18-3.19,  3.20,  3.21 

3.22-3.24 

LSQRAZ  (raise  expression  to  power)   2.21-2.22 

LSCJSBS  (substitute  in  expression)    2.24,  2.26,  3.31,  3.44,  3.45 

LSQTRC  (truncate  expression)   2.28-2.29,  3.45,  3.47-3.49,  3.52 

LSQTYP  (type  of  expression)   2.24,  2.35,  2.37-2.39,  3.53 

LSQUNM  (unbuild  list)    1.20,  2.17,  3.26 

LSQVAL  (evaluate  expression)   2.24-2.25,  2.37,  2.43,  3.26,  3.30 

3.31,  3.32,  3.35 

LSQVVL  (evaluate  expression)    3.30,  3.38.  3.39 

LSSCPY  (copy  list  or  expression)    2.3,  3.6,  3.26,  3.27 

LSTHQL  (compare  lists)   2.35,  2.39 

LVLRVl  (raise  list  reader  level)   3.51,  3.57 

LVLRVT  (raise  list  reader  level)   3.51,  3.57 

machine  address   1.24-1.25 

MADLFT  (locate  left  element)   3.51,  3.54 
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MADNBT  (locate  element  relative  to  list  bottom)   3.51,  3.55 

MADNTP  (locate  element  relative  to  list  top)   3.51,  3.55 

MAD0V  (machine  address  of)   1.25 

MADRGT  (locate  riglit  element)   3.51,  3.55 

MAKEDL  (assign  attribute-value  list)   3.60,  3.61 

MANY   (add  elements  to  list)    1.20,  2.2,  3.7 

M0N0   (provide  integer  name)    1.20 

MTDLST   (remove  attribute-values)    3.60,  3.61 

MTLIST   (empty  list)    3.26,  3.27,  3.28 

NAMEDL   (find  attribute-value  list)    3.60,  3.61 

NAMTST   (test  for  being  a  list)   3.51,  3.52 

NEWB0T   (add  to  bottom  of  list)    3.1,  3.7,  3.8 

NHWT0P   (add  to  top  of  list)    3.1,  3.8 

NEWVAL   (new  value  for  attribute)    3.60-3.61 

notation   1.3-1.11 

NULSTL   (break  off  elements  as  a  list)    3.1,  3.8,  3.26 

NULSTR   (break  off  elements  as  a  list)   3.1,  3.9,  3.26 

NUMPY  (multiply  expression  by  number  -  obsolete)    3.29 

NXTLFT   (insert  element  to  left)    3.1,  3.9 

NXTRGT  (insert  element  to  right)   3.1,  3.10 

operators,  hierarchy  of   2.8 

ordering   1.10,  1.12-1.14,  2.39-2.40,  3.59 

output  of  expressions   2.13-2.15,  2.19-2,20,  3.18-3.24 

overwritng,  execution  time   1.22 

0VWRT   (execution  time  overwriting  utility)    1.22 

parentheses   (represented  by  3H1.*)    1.7,  2.8,  2.43,  2.45,  2.47-2.48 

removal  of  2.49-2.54 

7IIPARTIAL   2.31 

P0PB0T   (pop  off  bottom  element)   3.26,  3.27 

P0PMID   (remove  an  element)    3.26,  3.27-3.28 

P0PT0P   (pop  off  top  clement)    3.26,  3.28 

power   1.5 

P0WER   (extract  power  from  expression)    3.51,  3.52 
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preprocessor   1.2,  3.62-3.73 

primitives,  SLIP   1.25-1.26 

PRIPUT   (output  expression  -  obsolete)    3.24 

PRLSTS   (dump  list  or  expression)   2.19 

PlITLST  (build  expression  -  obsolete)   3.1,  3.10 

QUKST,  preprocessor   1.2 
2HQ.    2,44 

RCF.LL   (return  cell)    3.26,  3.28 

reader,  list   3.5,S-3.58 

reader,  sequence   2.35-2.37,  3.5-3.6,  3.53-3.54 

REAL   (provide  floating  point  name)    1.20 

recursion   1.20-1.21 

REED  (read  reader)   3.51,  3.58 

RUN  2.3   (FORTRAN  compiler)    1.1,  1.15-1.18 

SAME   (provide  floating  point  name)   1.20 

SBARGS   (count  arguments  in  call)   1.23-1.24 

SBST  (substitute  in  expression  -  obsolete)   3.30,  3.45 

SCOPE  3.1.2  operating  system   1.1,  1.15-1.18 

SELECT  (satisfy  externals)   1.15 

SEQLL   (advance  sequence  reader  to  left)    2.35-2.37,  3.53,  3.54 

SEQLR  (advance  sequence  reader  to  right)   2.35-2.37,  3.53,  3.54 

SEQRDR   (sequence  reader)   2.35-2.37 

SEQSL   (advance  seq.  reader  structurally  left)   3.51,  3.53,  3.54 

SEQSR  (advance  seq.  reader  structurally  right)   3.51,  3.54 

simplification   see  ordering,  SMPL 

SLIP  list  processing  package   1.1,  1.3-1.5,  1.18 

SLIP  program   1 . 15 

SMPL  (simplify  expression)   1.12 

S0LVE   (solve  linear  equation)    3.51,  3.52 

STRDIR   (store  directly)   1.19,  1.22,  1.25 

STRINU   (store  indirectly)    1.25 

SUB   (subtract  expressions  -  obsolete)   3.29 

SUBSET  (substitute  at  list  bottom)   3.1,  3.10-3.11 
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SUBST   (substitute  in  list)    3.1,  3.11 

substitution   see  evaluation 

SUBSTP   (substitute  at  list  top)   3.1,  3.11 

SUET   (substitute  in  expression  -  obsolete)   3.30,  3.45 

SUMPY   (multiply  expressions  -  obsolete)   3.29 

SYMB66   ii 

term   1.5,  3.59 

TFRM  (terminate  VISIT)    1.20,  3.22,  3.38 

tests   see  analysis 

THEN  IF   see  IFTIIF.N 

T0P   (top  of  list)    3.51 

TRCAL  (truncate  expression  -  obsolete)   3.51,  3.53 

TRUNC   (truncate  expression  -  obsolete)   3.50 

truncation   2.28-2.29,  3.47-3.50 

TSTC0N   (test  if  expression  is  constant  -  obsolete)   3.51,  3.53 

unit,  input/output   1.16,  2.13 

VALARG   (for  use  in  scanning  argument  lists)   1.24 

VISIT   (recursive  use  of  FORTRAN  code)    1.20,  3.22,  3.36,  3.38, 

3.39,  3.42 
2HV.    2.44 
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This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Commission,  nor  any 
person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
express  or  implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  in  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
in  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  information,  apparatus, 
method,  or  process  disclosed  in  this 
report. 

As  used  in  the  above,  "person  acting  on  behalf 
of  the  Commission"  includes  any  employee  or 
contractor  of  the  Commission,  or  employee  of 
such  contractor,  to  the  extent  that  such  em- 
ployee or  contractor  of  the  Commission,  or 
employee  of  such  contractor  prepares,  dis- 
seminates, or  provides  access  to,  any  infor- 
mation pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with 
such  contractor. 
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